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Abstract 

A massive black hole resides in the center of most, perhaps all galaxies. The one in the 
center of our home galaxy, the Milky Way, provides a uniquely accessible laboratory for 
studying in detail the connections and interactions between a massive black hole and the 
stellar system in which it grows; for investigating the effects of extreme density, velocity 
and tidal fields on stars; and for using stars to probe the central dark mass and to probe 
post-Newtonian gravity in the weak- and strong-field limits. Recent results, open ques- 
tions and future prospects are reviewed in the wider context of the theoretical framework 
and physical processes that underlie them. 
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1 Introduction 



The massive black hole (MBH) in the center of the Milky Way is the nearest example of a universal 
phenomenon: central galactic MBHs. It was first detected as an unusual non-thermal radio source, 
Sagittarius A* (Sgr A*) (Balick & Brown 1974). Over the following decades, observations across the 
electromagnetic spectrum, together with theoretical arguments, established with ever-growing confi- 
dence that Sgr A* is at the dynamical center of the Galaxy and that it is associated with a very massive 
and compact dark mass concentration. This has ultimately led to the nearly inescapable conclusion that 
the dark mass is a black hole. At present, observations of Sgr A* and its immediate environment offer 
the strongest empirical evidence for the existence of MBHs. 

This review deals with one of the many aspects of the MBH phenomenon: stellar processes near 
the MBH in the Galactic center (GC). To motivate this choice of subject matter and to establish the 
astrophysical context and the scientific questions of interest, it is necessary first to address the questions: 
why study MBHs? why focus on stars near a MBH? and why in the GC? 

1.1 Astrophysical context 

Observations indicate that most nearby galaxies contain a massive compact dark object in their center, 
whose mass lies in the range 1O 6 M <wi< fewx 10 9 M Q (Kormendy & Richstone 1995; Magorrian et 
al. 1998; Gebhardt et al. 2003; see recent review by Ferrarese & Ford 2005). It is widely believed that 
these dark objects are MBHs, and that they exist in the centers of most, if not all galaxies. Their number 
density and mass scale are broadly consistent with the hypothesis that they are now-dead quasars, which 
were visible for a relatively short time in their past as extremely luminous Active Galactic Nuclei (AGN), 
powered by the gravitational energy released by the accretion of gas and stars 1 (Soltan 1982; Yu & 
Tremaine 2002). Some present-day galaxies have AGN, although none as bright as quasars. However, 
most present-day galactic nuclei are inactive, which implies that accretion has either almost ceased or 
switched to a non-luminous mode. Their inactivity is not due to the lack of gas supply; most galaxies 
have more than enough to continue powering an AGN. The "dimness problem" is one of the key issues 
of accretion theory, which deals with the physics of flows into compact objects (e.g. Narayan 2002). 

The MBH mass m is known to correlate very well with c sp h, the stellar velocity dispersion of the 
spheroidal component of the host galaxy (the spheroid is the bulge in disk galaxies, or the entire galaxy 

'it is possible that very low-mass MBHs like the one in the GC have acquired most of their mass by mergers 
with other black holes. 
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in ellipticals. Ferrarese & Merritt 2000; Gebhardt et al. 2000; Tremaine et al. 2002), 



m ~ 1.3 x 10' 
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Since the m/a^ correlation holds on large scales where the MBH is dynamically unimportant, it sug- 
gests a causal connection between the formation of the MBH and its host galaxy. There are also claims 
of correlations between cr sph and the circular velocity of spiral galaxies on very large scales, which im- 
plies a connection between the MBH mass and the dark matter halo of the galaxy (Ferrarese 2002; Baes 
et al. 2003). Whether it is the galaxy that somehow influences the growth of the MBH, or the MBH that 
somehow regulates the formation of the galaxy, is still unclear. The nature of the role that MBHs play 
in the formation and evolution of galaxies is one of the major unsolved problems of galaxy formation. 

The Galactic MBH is quite normal. Like most MBHs, it is inactive, and it follows the m/a sp h 
relation (Valluri et al. 2005). With m~ (3-4)xlO 6 M ,it is one of the least massive MBHs discovered, 
in keeping with the relatively small bulge of the Milky Way. What makes it special is its proximity. 
At ~ 8 kpc from the Sun, the Galactic black hole is ~ 100 times closer than the MBH in Andromeda, 
the nearest large galaxy, and ~2000 times closer than galaxies in Virgo, the nearest cluster of galaxies. 
For this reason it is possible to observe today the stars and gas in the immediate vicinity of the Galactic 
MBH at a level of detail that will not be possible for any other galaxy in the foreseeable future. 

In spite of its relative proximity, observations of the GC are challenging due to strong, spatially 
variable extinction by interstellar dust, which is opaque to optical-UV wavelengths. Emission pro- 
cesses related to the very low-level accretion activity in the GC can be observed in the radio and mm 
wavelengths on the long wavelength range of the spectrum, or in X-rays on the short wavelength range. 
However stars, whose spectra are approximately black-body and extend over a limited wavelength range, 
emit most of their luminosity in the optical-UV range, some in the IR and negligibly in longer or shorter 
wavelengths. Observations of stars in the GC must be conducted in the infrared, through the atmospheric 
transmission "windows" (when observing from Earth), primarily in the ^-band 2 (2.2 pirn). While the 
visual extinction along the line of sight to the GC is Ay ~ 30 mag , which corresponds to a transmis- 
sion ratio of only 1 : 10 12 photons, the /T-band extinction is only Ak ~ 3 mag (Rieke et al. 1989), or a 
transmission ratio of 1 : 15 photons. However, IR astronomy poses a much more difficult technological 
challenge than optical astronomy because of the strong ambient background 3 . Rapid advances in IR 
detectors in the last decade have finally made it possible to take full advantage of the proximity of the 
Galactic MBH. Individual stars can now be tracked as they orbit the MBH in the very crowded stellar 
field around it. Their mass, age and evolutionary stage can be determined from their spectral features, 
and their radial (line-of-sight) velocities can be measured from their Doppler-shifted spectra. 

Because of the huge mass ratio between a star and the MBH, stars orbiting near it are effectively test 
particles. This is to be contrasted with the gas in that region, which can be subject to non-gravitational 
forces due to thermal, magnetic or radiation pressure. These can complicate the interpretation of dy- 
namical data and limit its usefulness. The term "near" is taken here to mean close enough to the MBH so 
that the gravitational potential is completely dominated by it, but far enough so that the stars can survive 

2 The A'-band is centered on 2.2 [im and is 0.6 /im wide. The absolute (intrinsic) /f-band magnitude, Mk, is 
related to the monochromatic luminosity at 2.2 /im, Lk, by Mk=-2. 51og 10 Lx+ 84.245, for Lk in ergs -1 fim~ l . The 
apparent (observed) magnitude is related to the absolute one by K=Mk+DM+Ak, where the distance modulus 
is DMs51og IO (/?o/10pc)= 14.5 mag for a distance to the GC of /?o = 8kpc (Eisenhauer et al. 2003) and where 
Ak^3 mag is the /if -band dust extinction coefficient in the direction of the GC (Rieke et al. 1989). For the Sun, 
M m = 3 .4 1 and L KQ = 2 . 1 54 x 1 32 erg s" 1 ^m" 1 . 

3 History and economics also played a role. Optical astronomy had a long head-start over IR astronomy, which 
had to wait for the advent of modern electronics and cryogenics. In addition, the market for IR devices is much 
smaller than that for optical ones. As a result, the development of efficient IR detectors lagged decades behind 
optical detectors. 



T. Alexander / Physics Reports 000 (2005) 000-000 



5 



(that is, beyond the MBH event horizon, or beyond the radius where stars are torn apart by the tidal 
field). In this range, stars directly probe the gravitational field of the MBH as long as their trajectory 
does not take them so close to the MBH that dissipative processes, such as tidal heating, affect their 
motion. The event horizon of the MBH in the GC is much smaller than the tidal radius for most stars, 
and so effects due to General Relativity (GR) lead to deviations of only a few percents from Newtonian 
motion. To first order, the stellar orbits can be treated as Keplerian, which substantially simplifies the 
analysis. However, with accurate enough astrometric and spectroscopic observations it may be possible 
to detect post-Newtonian effects in the orbits and to probe GR. 

The stars near the MBH are of interest in themselves, in particular when they can no longer be 
treated as point particles due to strong interactions with other stars or with the MBH. The observed 
stellar density near the MBH is extremely high, in accordance with theoretical predictions of the stellar 
distribution near a MBH. The presence of so many stars so close to the MBH leads to a variety of 
interactions between the MBH and the stars, which may be relevant for feeding the MBH and can 
potentially produce an observable signal of the MBH existence. The high density, large orbital velocities 
and strong tidal field, well above those in any other Galactic environment, effectively make the central 
~0.1pc of Galaxy a "stellar collider". In this extreme environment, otherwise rare dynamical processes 
that can affect the inner structure of stars occur relatively frequently. The formation and evolution of 
stars in the extreme environment so close to the MBH challenge theories of stellar dynamics and star 
formation. 

The Galactic MBH and the stellar environment around it are probably representative of many sim- 
ilar systems in the universe. It is therefore likely that insights gained from the study of this uniquely 
accessible system can be applied to MBHs in general. In particular, insights about the role that stars 
near MBHs play in the feeding, growth and evolution of the MBH and the galactic nucleus. 

1.2 Science questions 

The study of the stars near the MBH in the GC offers an opportunity to address many specific questions 
in the general themes mentioned above. These broadly fall in four categories. 

The nature of the dark mass The identification of a BH is mainly done by eliminating all other pos- 
sibilities, and so an inquiry into the nature of the dark mass must begin by measuring its primary 
parameters: mass, size, position and velocity (§4). The measured compactness of the mass distri- 
bution can discriminate between a BH and other alternatives (§4.2). The central dark mass may be 
composed of several components that have to be disentangled, for example a binary MBH (§4.3), 
or a single MBH surrounded by a dense cluster of stellar mass compact remnants and a dark 
matter cusp of elementary particles (§4.2.2, §4.2.3). The position and velocity of the dark mass, 
which lies at the dynamical center of the Galaxy, are also primary parameters for reconstructing 
the Galactic structure and rotation and for calibrating distance indicators (§4.1.3). 

Post-Newtonian physics Stellar orbits very near the MBH should display post-Newtonian deviations 
such as periapse angle shifts, Lense-Thirring precession or gravitational redshift, and thus test 
various aspects of GR (§5.1). Compact remnants and low-mass main sequence (MS) stars in- 
spiraling into the MBH could emit detectable gravitational wave (GW) radiation (§6.2.2). The 
deflection of light by the gravitational potential of the MBH may be detected by gravitational 
lensing effects (§5.2). 

Formation and growth of MBHs Observations of stars can provide clues on the way the MBH forms 
and grows. Gas can be supplied by stellar winds from nearby giants or young massive stars, or 
by stars passing very near the MBH (§2.2, §6.2.4). The efficiency of tidal disruption and capture 
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can be probed by searching for stars that are in the process of being captured (§6.2.3) or have 
narrowly avoided such fate (§6.2.4). Stars can interact with a cold remnant accretion disk and 
reveal its existence (§6.2.5). The MBH may also grow by merging with other BHs, whether 
stellar mass BHs (SBHs) (§3.3) or intermediate mass black holes (IBHs) (§4.3, §7.2.3). 

Stars under extreme conditions the unique and extreme conditions near the MBH (§3) may be the 
reason behind the unusual stellar populations observed there. Observations may reveal evidence 
for processes such as extreme mass segregation (§3.3), collisional destruction (§3.4.1), tidal spin- 
up (§3.4.2) or 3-body exchanges (§7.2.4) and provide new insights about unusual modes of star 
formation near a MBH (§7.2.1). 

1.3 Scope and connections to related topics 

The GC is a complex interacting system. Its different components, processes and regions are interre- 
lated. Many of the topics discussed below can be expanded to apply to MBHs in general. Keeping 
this review focused on stellar processes near the Galactic MBH inevitably involves setting somewhat 
arbitrary limits to its scope. 

This review deals with stellar processes in the inner 1-2 parsecs of the GC, where the MBH domi- 
nates the dynamics, and in particular in the innermost region where stars can still exist. Stellar processes 
near a MBH may result in observable effects and thus provide methods for detecting MBHs. However, 
some of the effects discussed below are rather subtle and their usefulness in investigating stellar dynam- 
ics near the MBH is limited to the GC, as it is unlikely that they could be observed with comparable 
precision near MBHs in other galaxies. 

In spite of the fact that stars are born of gas and dust and to gas and dust (and compact remnants) 
return, gas processes, such as accretion and outflows, will be discussed only to the extent that they are 
directly related to stellar processes near the MBH. Observations of the Milky Way and a few nearby 
galaxies show that distinct stellar populations exist near their MBH (Lauer et al. 1998). It is likely 
that the environment near a MBH affects the mechanisms of star formation. However, star formation is 
poorly understood even in typical galactic environments, let alone near a MBH, and will be discussed 
here only briefly. 

The discussion will focus on low mass MBHs (m< 10 8 M Q ), like the Galactic MBH, where the 
tidal disruption radius (for solar type stars) lies outside the event horizon. This allows the possibility 
of tidal disruption and limits the role of GR in the dynamics of stars near the MBH. This situation is 
not generalizable to massive MBHs (m> 10 8 M Q ), which are qualitatively different in these respects. 
At the other end of the mass spectrum lie IBHs (10 2 <m < 10 4 M Q ), whose existence is still a matter of 
speculation (e.g. Kaaret et al 2001; Hopman, Portegies-Zwart & Alexander 2004; Portegies Zwart et al. 
2004; see review by Miller & Colbert 2004). IBHs will also not be discussed here, except where related 
directly to the GC. Likewise, although major (comparable mass) MBH mergers and binary MBHs may 
play an important role in the evolution of MBHs in general (e.g. Merritt & Milosavljevic 2004), it seems 
unlikely that this is relevant for the low mass Galactic MBH and the steep stellar cusp around it, and so 
this topic will be mentioned only briefly here. 

Additional information and discussion of some of the topics outside the scope of this review can be 
found in other reviews about the GC: Genzel & Townes (1987; early general review on inner 10 pc); 
Genzel, Hohlenbach & Townes (1994; comprehensive review with focus on stellar and gas processes); 
Morris & Serabyn (1996; focus on large spatial scales, gas and star formation); Melia & Falcke (2001; 
focus on accretion and the radio source Sgr A*); Yusef-Zadeh, Melia & Wardle (2000; focus on interac- 
tions between stars, gas, dust, clouds and supernovae remnants) and Alexander (2003; focus on stellar 
phenomena). More observational and theoretical background can be found in Falcke & Hehl (2003). 
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The acronyms and notations used in this review follow where possible the usual conventions in the 
literature. For convenience, the ones frequently used here are listed in table (1.1). 

2 Observational overview: Stars in the Galactic center 

The observational picture of the GC, and in particular that of the central few parsecs, is a rapidly evolving 
one. A too detailed focus on the most recent observations runs the risk of premature obsolescence, as 
ongoing observational campaigns continuously yield new data and insights. This overview will therefore 
be limited to painting a broad-brush picture, and defer detailed discussion to the relevant topical sections 
below. 

It is useful to begin by a brief summary of some stellar properties that are of particular relevance in 
this context. The initial stellar mass sets the stellar lifespan, and together the two determine which dy- 
namical processes are relevant (for example, a slow dynamical process that operates on massive objects 
may be irrelevant because the massive stars are too short-lived, while the long-lived stars are not massive 
enough). The initial mass also sets the evolution of the stellar IR luminosity, which determines if, and 
at which evolutionary stage the star can be observed. Figure (2.1) shows the stellar lifespan, the MS 
K-bwA luminosity and the MS spectral type designation as function of initial stellar mass, compared to 
the estimated dynamical two-body relaxation time in the GC and compared to current and future K-bnn& 
detection limits (for evolved stars, see Fig. 3.5 in Genzel, Hollenbach & Townes 1994). 

2.1 The central 100 parsecs 

The stars well outside the dynamical sphere of influence of the Galactic MBH on the < lOOpc scale are 
relevant for understanding the stellar environment in the inner parsec because they define the "boundary 
conditions" for the inner GC; because, given enough time or under certain dynamical conditions (§3.3) 
stars can migrate from large distances to the center, and because the stellar population far from the MBH 
can serve as a control sample for separating the effects of the MBH from intrinsic stellar phenomena. 

The GC lies at the center of the Galactic bulge (typical length-scale ~2 kpc), which is composed of 
an old (~ 10 Gyr), passively evolving stellar population. However, the stellar population in the central 
~ 100 pc, whose density rises toward the center as r~ 2 , appears to have been forming continuously over 
the lifetime of the Galaxy. The star formation is probably fed by the large reservoir of molecular gas 
(the central molecular zone) that is concentrated in the inner ~200 pc (Serabyn & Morris 1996). The 
stellar population in the GC is thus a mixture of the inner extension of the old bulge population and 
intermediate-age and young stars from more recent star formation epochs (Philipp et al. 1999; Mezger 
et al. 1999; Figer et al. 2004). 

About half of the young stars in the region are found today in three particularly massive, young 
(<5 Myr) clusters: The Quintuplet (~ 30 pc from the center in projection, M~ 10 4 M Q , p~fewx 
lO 3 M pc~ 3 ), The Arches (~30 pc from the center in projection, M> 1O 4 M , 7?~0.2 pc, p~3x 
10 5 M Q pc~ 3 ), and the central cluster around the MBH (Figer 2003). The three clusters contain in total 
hundreds of MS O-stars, tens of Wolf-Rayet (WR) stars and a few luminous blue variable stars (§2.2), 
which are ~ 10% of all the massive stars (initial mass>20M Q ) in the entire Galaxy. They produce 
~ 10% of the Galactic ionizing luminosity, and are responsible for the formation of 0.01 M & yr _1 in stars 
(a rate per volume of ~ 10~ 7 MQyr _1 pc~ 3 , 250 times higher than the Galactic mean). This is however 
only ~ 1% of the total Galactic star formation rate. The disparity reflects the marked bias in the initial 
mass function of these clusters toward the formation of massive hot and luminous stars (Figer et al . 
1999; Figer 2005). The overall similarity between the young stellar populations of the Arches and the 
Quintuplet and that of the central cluster (§2.2) is noteworthy in view of the fact that the two off-center 
clusters show no evidence of harboring a central BH. 
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Table 1.1: Acronyms and notation frequently used in this review (by topic, alphabetically) 



Term 


Defined 


Meaning 


Term 


Defined 


Meaning 


Acronyms 




Stellar 


dynamics 




AGN 




Active Galactic Nucleus 


J 


p. 18 


Specific angular momentum 


BH 




Black hole 


n± 


p. 16 


Stellar number density 


DF 




Distribution function 


''coll 


p. 15 


Collisional radius 


GC 




Galactic center 


tc 


p. 16 


Collision time 


GR 




General Relativity 


t e 


p. 16 


Evaporation time 


GW 




Gravitational waves 


tj 


p. 59 


Angular mom. relaxation time 


IBH 




Intermediate mass black hole 


t r 


p. 16 


Relaxation time 


IMF 




Initial mass function 


t s 


p. 16 


Mass segregation time 


IR 




Infrared 


a 


p. 19 


Logarithmic slope of n* 


MBH 




Massive black hole 


e 


p. 17 


Specific orbital energy 


MS 




Main sequence 


P 


p. 20 


Mass density 


NS 




Neutron star 


a 


p. 15 


ID stellar velocity dispersion 


SBH 




Stellar black hole 


v> 


p. 17 


Relative gravitational potential 


WD 




White dwarf 








WR 




Wolf-Rayet star 


Stellar astrometry and orbits 








a 


p. 18 


Keplerian semi-major axis 


Stellar properties 


e 


p. 18 


Keplerian eccentricity 




p. 27 


Stellar binding energy 


i 


p. 36 


Inclination angle 


K 


p. 4 


K-band magnitude 


P 


p. 33 


Projected distance from MBH 


U 


p. 25 


Total stellar luminosity 


P 


p. 18 


Keplerian orbital period 


M* 


p. 14 


Stellar mass 


r P 


p. 27 


Orbital periapse 




p. 17 


Number of stars 


r a 


p. 35 


Orbital apoapse 


R< 


p. 14 


Stellar radius 


Z 


p. 33 


Line of sight distance from MBH 


U 


p. 9 


Stellar lifespan 




p. 46 


Velocity in terms of speed of light c 


r* 


p. 25 


Effective temperature 


T 


p. 46 


Relativistic parameter 


V e 


p. 15 


Stellar escape velocity 










p. 61 


Stellar circular velocity 


Star-MBH interactions 




p. 27 


Stellar dynamical timescale 




p. 57 


Tidal penetration parameter 








n 


p. 15 


Tidal disruption radius 


MBH properties 




r, 


p. 75 


Tidal scattering rate 


in 


p. 14 


MBH mass 


r f 


p. 59 


Tidal disruption rate 


n, 


p. 15 


MBH radius of influence 


t p 


p. 27 


Periapse passage timescale 


is 


p. 15 


Schwarzschild radius 








s 


p. 49 


Spin parameter 


Gravitational lensing 








A 


p. 51 


Lensing magnification 


GC properties 




D LS 


p. 50 


Lens-source distance 


A K 


p. 4 


K-band dust extinction 


D L 


p. 50 


Observer-lens distance 


DM 


p. 4 


Distance modulus 


Dos 


p. 50 


Observer-source distance 


Ro 


p. 4 


Sun-GC distance 


X 


p. 51 


Image angular position in terms of 9e 


t H 


p. 15 


Hubble time (~age of GC) 


y 


p. 51 


Source angular position in terms of 9e 


@o 


p. 39 


Local Galactic rotation speed 


r L 


p. 52 


Gravitational lensing rate 








Oe 


p. 50 


Einstein angle 
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Figure 2.1: Main sequence lifetime and i^-band magnitude for zero-age solar metallicity stars (based 
on the stellar tracks of Schaller et al. 1992). An approximate association of mass with spectral type 
is indicated in the upper horizontal axis (Cox 2000). The apparent T^-band magnitudes are given for a 
plausible range of distance and extinction coefficient values: the top curve for the fiducial distance to 
the GC of /?o=8 kpc (Reid 1993; Eisenhauer et al. 2003) and K-band extinction of Ak =3 mag (Rieke et 
al. 1989), and the bottom curve for the more recent determination of Ro = l .62 kpc and =2.8 mag. A 
rough estimate of the two-body relaxation time in the GC, ( r ~3x 10 9 yr (Alexander 1999), the current 
detection limit in the i^-band with the VLT (K=IS mag with completeness of 63%, Genzel et al. 2003b) 
and the anticipated limit for the VLTI interferometer with adaptive optics (Paresce et al. 2003) are also 
shown. 



X-Ray observations of the inner 20 pc also reveal ~ 2000 X-ray sources, which follow the distri- 
bution of the stars in the IR. The hard spectral index of over half of these sources suggests that they 
are magnetically accreting white dwarfs (WDs) and X-ray pulsars in high mass X-ray binaries (where a 
neutron star (NS) accretes the wind of a massive companion). 

It is useful to construct a stellar population model to represent the continuously star forming popula- 
tion of the central GC. When the star formation rate is constant, the present-day mass function (number 
of stars per stellar mass) at time t is simply related to the initial mass function (IMF) dW*/dM*| by 



oc 



min[*,k(M*)], 



(2.1) 



where t* is the stellar lifespan. The stellar contents of such a model are listed in table (2.1). This model 
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Table 2.1: A mean stellar population model for the inner GC 



Component 


Mass fraction b 


Number fraction c 


Notes 


Gas consumed 


1.8 


— 


Including reprocessed gas 


Gas reprocessed 


0.8 


— 


Assumed ejected from GC 


Live stars (M* > 0. 1 M G ) 


0.74 


0.85 


Mean mass 0.43 M Q 


Live stars (M* > O.8M ) 


0.22 


0.13 


Mean mass O.84M 


Compact remnants d 


0.27 


0.15 


Mean mass O.94M 


WDs e (O.6M ) 


0.03 


0.03 


Progenitor mass 0.8-1. 5M 


WDs e (O.7M ) 


0.08 


0.06 


Progenitor mass 1.5-2. 5 M 


WDs e (1.1M ) 


0.12 


0.05 


Progenitor mass 2.5-8. OM 


NSs e (1.4M ) 


0.03 


0.01 


Progenitor mass 8.0-30M Q 


SBHs / (1OM ) 


0.01 


0.0005 


Progenitor mass > 3OM 



" Continuous star formation over 10 Gyr at a constant rate with a Miller-Scalo IMF (Miller 

& Scalo 1979) and solar metallicity (Alexander & Sternberg 1999). Calculated by a stellar 

population synthesis code (Sternberg, Hoffman & Pauldrach 2003) using the Geneva stellar 

evolution tracks (Schaller et al. 1992). 

b Relative to total remaining mass in stars and remnants. 

c Relative to total remaining number of stars and remnants. 

d Average supernovae rate of 1.2 x 10~ 5 yr _1 per 1O 7 M of gas consumed. 

e Meylan & Mayor (1991). 

1 Timmes, Woosley & Weaver (1996). 



assumes a 10 Gyr old, continuously star forming population with a Miller-Scalo 4 IMF (1979) and solar 
metallicity. The ^-band luminosity function (number of stars per magnitude) of this model is well 
approximated by a power-law down to faint magnitudes (K<2§ mag), A\ogN- k /AK = s, with .s~0.35 
(Alexander & Sternberg 1999). This value broadly agrees with the observed T^-band luminosity function 
in the inner 0.25 pc (Davidge et al. 1997) and that of the Galactic bulge as observed through Baade's 
window, ~0.5 kpc from the GC (Tiede, Frogel & Terndrup 1995; Zoccali et al. 2003). 

2.2 The central parsec 

The central parsec lies inside the dynamical sphere of influence of the MBH (§3.1.2). The properties 
of the stellar population there are observed to change with decreasing distance from the MBH, once at 
~0.4 pc and then again at ~0.04 pc (Fig. 2.2). 

Both the stellar spatial distribution and the stellar population outside ~ 0.4 pc are a continuation 
of the large scale n+ oc r~ 2 mixed distribution of old and young stars in the inner 100 pc (Genzel et al. 
2003b). The old population is dynamically relaxed and follows the Galactic rotation pattern (McGinn 
et al. 1989; Haller et al. 1996; Genzel et al. 1996). 

Inside ~0.4 pc the stellar spatial distribution flattens to ocr -1 4 (Genzel et al. 2003b). On ap- 
proximately the same length-scale (whether or not this is a coincidence is unclear), there appears a 
dominant population of massive young stars (Forrest et al. 1987; Allen, Hyland & Hillier 1990; Krabbe 
et al. 1995; Paumard 2004). This population consists of dynamically unrelaxed blue supergiants with 
distinctive helium emission lines in their infrared spectra (Tamblyn 1996; Genzel et al. 1996; Genzel 

4 A useful analytic approximation of the normalized Miller-Scalo IMF (1979) is In (dW*/dM*| ) = 
-ln(2MjM )ln(4OOM*/M )/4-7.11 for O.1<M*/M < 125. 
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Figure 2.2: A false color composite image of the central parsec of the GC taken by the Very Large 
Telescope in the H (1.65/im), K s (2.16/Ltm) and U (3.76/im) infrared bands (Genzel et al. 2003b). The 
position of the MBH is indicated by the two small arrows. A mixture of stellar types is clearly apparent: 
Old, intermediate mass cool red giants and young massive hot blue stars. (Reprinted with permission 
from the Astrophysical Journal) 



et al. 2000). The lack of hydrogen lines and evidence of outflow in the spectral line shapes ("P Cygni 
profiles") identify them as WR stars or Luminous Blue Variables (LBVs) (Najarro et al 1997). WRs 
are stars with masses of order fewx 10M Q , lifespans of fewx 10 6 yr undergoing rapid mass loss by a 
stellar wind, which removes the hydrogen-rich envelope and uncovers the helium-rich core. LBVs are 
WR progenitors that are still in the process of losing their envelopes. There are tentative detections of 
X-ray emission from the shocked winds of a couple of these stars (Baganoff et al. 2003). 

The emission line stars co-rotate tangentially in an opposite direction to the Galactic rotation, in 
two partially overlapping disk-like structures (Levin & Beloborodov 2003; Genzel et al. 2003b; Fig. 
2.3). The two disks are strongly inclined relative to each other. In projection on the plane of the sky, the 
inner disk appears to rotate clock-wise, and the outer disk counter-clockwise. The inner extent of the 
disk population is ~0.04 pc. Two prominent concentrations of stars appear among the disk stars. The 
IRS 16 complex is located on the line-of-nodes (the intersection between the disk plane and the plane of 
the sky) of the clockwise disk, and so this apparent stellar over-density can plausibly be explained as 
a projection effect (Genzel et al. 2005, in prep.). The IRS 13 complex is part of the counter-clockwise 
disk, which is nearly face-on. IRS 13 appears to be a genuine over-concentration of stars with common 
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Figure 2.3: Evidence for disk-like stellar structures around the Galactic MBH. Top left: Distribution in 
projected velocity space of 13 bright emission line stars in the inner ~5" from a compilation by Genzel 
et al. (2000). Top right: The same, after transforming to the disk coordinates. 10 of the 13 stars lie in 
a thin disk with opening angle < 10°. (Levin & Beloborodov 2003. Reprinted with permission from 
the Astrophysical Journal). Bottom: Reconstructed orientations of the the inner clockwise ring and the 
outer counter-clockwise disk (Genzel et al. 2003b. Reprinted with permission from the Astrophysical 
Journal), shown from two vantage points. The coordinates x and y are in the plane of the sky, and z is 
along the line of sight, with the observer at z=-oo. 



velocities (Maillard et al. 2004; Schodel et al. 2005). 



Inside ~ 0.04 pc the population changes again. There are no bright giants, red or blue, and only 
faint blue stars are observed there. This population is know as the "S-stars" or "S-cluster", after their 
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identifying labels 5 . Deep near-IR photometric (Krabbe et al. 1995; Genzel et al. 2003b) and spectro- 
scopic (Genzel et al. 1997; Eckart, Ott & Genzel 1999; Figer et al. 2000; Gezari et al. 2002; Ghez et al. 
2003b) observations of that region were all consistent with the identification of these stars as massive 
MS stars. It is now spectroscopically confirmed (Ghez et al 2003b; Eisenhauer et al. 2005) that the 
S-stars are B0-B9 MS stars (corresponding to masses in the range ~3M Q - ~ 15M Q ). The rotational 
velocities (stellar spin) are similar to those of Solar neighborhood stars of the same spectral type. Thus, 
there is no indication of anything unusual about the S-stars, apart for their location very near the MBH. 
Astrometric observations (Ghez et al. 2005; Schodel et al. 2003; Eisenhauer et al. 2005) indicate that, 
unlike the emission line stars farther out, the S-stars orbits appear random (isotropic) and uncorrected 
with the orientation of either disk. There are some marginal statistical indications of a bias toward radial 
orbits (Genzel et al. 2000; Schodel et al. 2003, but see Ghez et al. 2005). 

The brightest star in the S-cluster, S2/S0-2 (K= 13.9 mag), is a transitional O8-B0 star with a mass 
of ~ 15M Q , effective temperature of ~ 30,000 K, intrinsic bolometric luminosity of L~ 1O 3 L and 
a main sequence lifespan of ~ 10 7 yr. All the other stars in the cluster are less massive, longer-lived, 
cooler and fainter. A B9 MS star has a mass of ~3M , a radius of ~2.1/? and a main sequence 
lifespan of ~4x 10 8 yr. A typical star would correspond to a B2 main sequence star, with a mass of 
~ 10M Q , a radius of ~4.5/? Q and main sequence lifespan of ~2x 10 7 yr (Cox 2000; Schaller et al. 
1992). It is not clear what is the relation, if any, between the population of the He stars on the 0.04-0.4 
pc scale and the S-cluster stars inside 0.04 pc. The nature and origin of the young stars near the MBH 
are the subject of ongoing investigations (§7). 

3 Stellar dynamics at extreme densities 

The extreme conditions in the central ~0.1 pc around the MBH are unique in the Galaxy. The strong 
tidal field (total disruption of a solar type star at a distance ten times the event horizon), high velocities 
(up to ~0.2c at disruption) and high stellar density (up to ~ 10 8 M Q pc~ 3 , 10 9 higher than in the vicinity 
of the Sun) effectively make the environment near the MBH a "stellar collider" (a total stellar collision 
rate of O(10 -5 ) yr 1 in the central 0.02 pc). Under such extreme conditions, strong stellar interactions 
may probe and modify the stars' internal structure, and otherwise rare dynamical processes may become 
frequent enough to have observable consequences. Before discussing the evidence for the existence of 
a high density cusp in the GC (§3.2) and its consequences for stellar collisions (§3.4), it is useful to 
review briefly the physical processes and scales that govern the evolution of a stellar system around a 
MBH. 

3.1 Physical processes and scales 
3.1.1 Dynamical processes 

The dynamical processes in a gravitating stellar system can be summarized by classifying stellar inter- 
actions by their distance scale (e.g. Binney & Tremaine 1987; Heggie & Hut 2003). In this context the 
term interaction refers both to remote gravitational interactions where the stars exchange momentum 
and energy, as well as actual physical contact between stars. 

5 Two labeling conventions are used for the IR sources near the MBH, the MPE labels (of the form S#) and the 
UCLA labels (of the form S#-#, see Ghez et al. 1998). The dual labels of some stars of special interest very near 
the MBH are listed here for convenience. 

MPE: SI S2 S4 S8 S9 Sll S12 S13 S14 

UCLA: S0-1 SO-2 S0-3 SO-4 S0-5 SO-9 SO-19 SO-20 SO-16 
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On the largest scale, the stellar orbit is determined by the superposed interactions with the MBH 
and all the other stars, which are well approximated by a smooth gravitational potential. On a shorter 
scale, the stellar motions are randomized by two-body scattering when two stars approach each other so 
that their mutual interaction dominates that of the smoothed potential. Over time, scattering leads to the 
relaxation of the system and the redistribution of energy and angular momentum. Massive perturbers, 
such as stellar clusters, IBHs or molecular clouds can accelerate the relaxation (Zhao, Haehnelt & Rees 
2002). In addition, very near the MBH, where the orbits are nearly Keplerian, angular momentum re- 
laxation (but not energy relaxation) can be accelerated by the process of "resonant relaxation" (Rauch & 
Tremaine 1996). In the course of relaxation, the stars, whose mass range spans 2-3 orders of magnitude, 
are driven toward energy equipartition. However, equipartition cannot be achieved in a self gravitating 
system, and in particular in the steep potential well of a MBH. A two-body interaction will tend to slow 
down the heavier star and speed up the lighter one, but since the typical orbital radius depend only on the 
velocity (specific energy), the heavy star will sink to the center, while the lighter star will drift outward. 
Over time, this process leads to mass segregation — the more massive stars are concentrated near the 
MBH and the lighter stars are pushed out of the inner region. 

The system continuously loses stars by two-body scattering, either abruptly, by a single strong en- 
counter that ejects a star, or gradually, by diffusion to higher energies (evaporation). The lost star takes 
away positive energy from the system. The system then becomes more bound and compact, the collision 
rate increases, yet more stars are lost, and a runaway process is launched. This so called "gravother- 
mal catastrophe" or "core collapse" reflects the fact that a self-gravitating system has a negative heat 
capacity — it become hotter when energy is taken out. Core collapse, if unchecked, will lead to the 
formation of an extremely dense stellar core surrounded by a diffuse extended halo. 

Once the density becomes high enough, very short-range inelastic collisions are no longer extremely 
rare, and the fact that the stars are not point masses, but extended objects with internal degrees of 
freedom, becomes significant. In such collisions energy is extracted from the orbit and invested in the 
work required to raise stellar tides, or strip stellar mass. The outcome of inelastic collisions depends 
critically on whether a MBH is present or not. If the collision is slow, as it is in the core of a globular 
cluster where there is no MBH, then the typical initial orbit is just barely unbound (parabolic). In this 
case, the tidal interaction may extract enough orbital energy for "tidal capture", and the formation of a 
binary. However, when the stars orbit a central MBH, the collisions are fast (The Keplerian velocity near 
the MBH exceeds the escape velocity from the star) and the initial orbits are very unbound (hyperbolic). 
Even very close fly-bys cannot take enough energy from the orbit to bind the two stars, and so they 
continue on their way separately after having extracted energy and angular momentum from the orbit. 
The stars can radiate away the excess heat on a time scale shorter than the mean time between collisions, 
but it is harder to get rid of the excess angular momentum. It is therefore likely that high rotation is the 
longest-lasting dynamical after effect of a close hyperbolic encounter, and that stars in a high density 
cusp are spun-up stochastically by repeated collisions (§3.4.2). Finally, at zero range, almost head- 
on stellar collisions can lead to the stripping of stellar envelopes (§3.2), the destruction of stars, or to 
mergers that result in the creation of "exotic stars". These are stars that cannot be formed in the course 
of normal stellar evolution, such as a Thorne-Zytkow object, which is an accreting neutron star (NS) 
embedded in a giant envelope (Thorne & Zytkow 1975). 

3.1.2 Characteristic scales 

Several important length-scales and timescales govern the dynamics of the stellar system around the 
MBH. These are listed here, with estimates of their values in the GC. A star of mass M*= 1M Q and 
radius R* = IR Q and a MBH mass of m = 3.5 x 1O 6 M (Schodel et al. 2003) are assumed. Physical 
lengths are expressed also as angular sizes for a distance to the Galactic Center of Rq = 8kpc, which 
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corresponds to 1"=0.039 pc (Eisenhauer et al. 2003). 

Event horizon The size of the event horizon of a non-rotating black hole, the Schwarzschild radius, is 

2Gm 1 9 _ 7 
rs = — 9 — = 10 cm ~ 3 x 10 pc ~ 9/xarcsec . (3.1) 

c l 

Tidal radius The tidal radius, r t , is the maximal distance from the MBH where the tidal forces of 
the MBH can overwhelm the stellar self-gravity and tear the star apart. Its exact value depends 
somewhat on the stellar structure and the energy of the orbit, but up to a factor of order unity it is 



( ^-J =10 13 cm~3x 10" 6 pc ~90/^arcsec. (3.2) 

Collision radius The collision radius, r co n, is the minimal distance from the MBH where large angle 
deflections by close gravitational encounters are possible. Closer to the MBH the relative velocity 
between the interacting stars, v, exceeds the escape velocity from the stellar surface, V e , and the 
distance of closest approach required for a large angle deflection becomes smaller than the stellar 
radius (i.e. when V 2 ~ GM+/R* < Gm/r ~ v 2 ). Stars undergoing such a physical collision are 
destroyed since the orbital energy exceeds the stellar binding energy. At distances smaller than 

(171 \ 
— J 7?*~0.08pc ~ 2arcsec, (3.3) 

two-body relaxation gradually becomes inefficient. 

MBH radius of influence The radius of influence, r/ !; is the region where the MBH dominates the 
dynamics. It is customary to define rh = Gmj a 1 , where a is the ID stellar velocity dispersion 
away from the MBH. However, in most cases this is not a well defined quantity because a depends 
on radius (a is constant only in an isothermal density distribution, which is not the case in the 
GC). Kinematically, the MBH potential dominates that of the stars out to a distance where the 
enclosed stellar mass is M*(< r) ~ 0(m). For the isothermal distribution M*(< r/,) = 2m. The 
measured enclosed mass in the GC indicates that M*(<2pc)~ra and M*(<4pc)~2m (Schodel 
et al. 2003), and so the MBH radius of influence in the GC can be estimated at 



ri, 



10 19 cm ~ 3 pc ~ 90 arcsec . (3.4) 



Age of GC It is assumed that the age of the GC, being the center of the old spheroid component of the 
Galaxy, the bulge, is of the same order as the age of the universe (roughly, the Hubble time), 

f ff ~10Gyr. (3.5) 
Dynamical timescale The dynamical, or orbital time tj, is the time it takes a star to cross the system 



r r 



.3 



t d ~ - ~ 2vrW 2 x 10 5 yr(at3pc) ~ 300yr(at0.03pc) , (3.6) 

v V GM(<r) 



where r is the typical size of the system and M(< r) is the total enclosed mass (MBH and stars). 
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Two-body relaxation timescale The two-body relaxation time, t r , is related to a, the mean stellar mass 
(M^) and the stellar number density by (e.g. Binney & Tremaine 1987) 

tr ^1 O(10 9 yr), (3.7) 

where a Maxwellian velocity distribution is assumed and log A ~ 0(10) is the Coulomb logarithm. 
A is the ratio between the largest and smallest impact parameters possible in the system for elastic 
collisions 6 . Since near the MBH cr 3 oc r~ 3 / 2 , the relaxation time is not a strong function of r in 
a stellar cusp where roughly n* oc r~ 3 / 2 , as is the case in the GC (§3.2). The old stars near the 
MBH are expected to be dynamically relaxed since since t r < tu- The relaxation time can also be 
approximately expressed in terms of iV*, the number of stars enclosed within radius r , 

ar) ~ V(MjJ 27rlog(0.4A^' M 

where P is the circular orbital period. 

Segregation timescale The mass segregation timescale for stars of mass M+ is of the same order as the 
relaxation timescale, and scales wifhM* as (§3.3) 

Evaporation time The evaporation time scale is 

t e ~300t r . (3.10) 

Evaporation is unimportant in the GC since t e 3> ?// (and also because the GC is not an isolated 
system, but is embedded in the bulge, so evaporated stars can be replaced). 

Collision timescale The rate, per star, of grazing collisions between two stars of mass and radius M", 
Rl and Af* R%, each, is 



t; 1 =4^n*a (RZ+Rl 



i | G(M«+M>1) 
2a 2 {R°+Rl) 



10 _9 yr (at0.02pc), (3.11) 



where it is assumed that the stars follow a mass independent Maxwellian velocity distribution 
with ID velocity dispersion a (this is a good approximation near the MBH, see §3.1.3). The two 
terms contributing to the total rate are the geometric cross-section and the gravitational focusing 
term due to the mutual attraction of the colliding stars, which dominates when the typical stellar 
velocities are much smaller than the escape velocity from the stars. For M"=M^, this corresponds 
to a 2 < GM+/2R+ = V e 2 /4. The collision rate diverges for a/V e — >0 or a/V e -^oo and is lowest 
when a/V e = 1 /2. For the high velocities near the MBH, gravitational focusing is only important 
for collisions with compact remnants where R*<^R & . 



6 The relaxation time can be derived approximately from the rate for large deflection angle collisions, t r 1 
n+cS, where the cross-section is given by S ~ ir (G (M*) /a 2 ) . The actual rate is larger by a factor log A due I 
the additional contribution of the many weaker collisions at large impact parameters. 
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3.1.3 Keplerian power-law stellar cusps 

Theoretical studies of the co-evolution of a MBH and the stellar system around it indicate that a wide 
range of dynamical scenarios all lead to the formation of a stellar density cusp, a region of diverging 
density around the MBH (e.g. Bahcall & Wolf 1976; 1977; Young 1980 Lee & Goodman 1989; Quinlan, 
Hernquist & Sigurdsson 1995). The dynamical evolution of the stars in the radius of influence depends 
on the ratio between the relaxation time and the age of the system. If t r >tu, two-body relaxation can 
be neglected, the MBH grows adiabatically and the resulting stellar distribution depends on the initial 
conditions. For example, an initially isothermal distribution will give rise to a oc r~ 3 / 2 cusp (Young 
1980). Other initial conditions will result in a variety of stellar distributions, some falling as steeply 
as 7i + ocr~ 5 / 2 (Lee & Goodman 1989; Quinlan, Hernquist & Sigurdsson 1995). However, when ?,-<?#, 
two-body relaxation erases the initial conditions, and the final configuration will depend only on the 
boundary conditions and the stellar mass function. 

The formation history of the Galactic MBH is not known. Estimates of the stellar density and 
velocity indicate that t r <tn near the MBH (§3.1.2), and so the old stars there are expected to have 
relaxed dynamically by two-body interactions. This is consistent with their observed isotropic velocity 
field and spatial distribution (Genzel et al 1996; Ghez et al 1998; Genzel et al 2000; §2). The DF 
of a spherically symmetric system with an isotropic velocity field is a function of e only. Bahcall & 
Wolf (1977) show that the distribution function (DF) / (stars per drdv interval) of such a system that 
is relaxed and consists of stellar masses in the range M\ < < M 2 and a moderately varying mass 
function dJV*/dM*, can be approximated near the MBH by a power-law with a mass-dependent index 

/(e;2lf*)oce», Pm=^, (3-12) 

where e=-v 2 /2 + ^(r) is the specific energy of the star, tp=-(f)>0 is the relative gravitational potential, 
<fi is the gravitational potential, and f=0 for e < 0. The power-law index varies in the range 0<pm< 1 /4 
across the mass spectrum. The corresponding stellar number density (stars per dr) is (see Eq. 3.17) 

n+(r;M+) oc r ~ 3/2 ~ PM . (3.13) 

Thus the most massive component of a relaxed multi-mass system is expected to have an r~ 7 /4 density 
profile, and the lightest one an r~ 3 / 2 density profile. In the simple case of a single mass population 
{Mi =M 2 ), Pm = 1/4 and n± oc r -7 / 4 (Bahcall & Wolf 1976; see also Binney & Tremaine 1987, p. 547 
for a simple derivation by dimensional analysis). Note that these results apply also to the Galactic dark 
matter halo elementary particles in the GC (Gnedin & Primack 2004), which are expected to interact 
gravitationally with the cusp stars and settle into a r~ 3 / 2 cusp near the MBH (§4.2.3). 

Dynamical simulations (e.g. Murphy Cohn & Durisen 1991; Baumgardt, Makino & Ebisuzaki 
2004) confirm that when the population consists of a single stellar mass, or when the most massive stars 
come to dominate the central stellar mass density due to mass segregation, their DF has indeed a power- 
law index /?m~0.25. The simulations also indicate that the central distribution of the less massive stars 
is flatter, with /?Af<0. Note that Eq. (3.12) was derived under various assumptions about the boundary 
conditions far from the MBH and about the stellar mass function, and its general applicability has yet to 
be investigated in full by numeric simulations. Nevertheless, it will be assumed here to hold generally, 
since it is broadly consistent with the numerical results and is analytically convenient. Note that the 
applicability of the Bahcall-Wolf solution to the GC is based on the assumption that the GC has evolved 
in isolation for at least a relaxation time before the present epoch. In particular, it is assumed that it has 
not experienced a strong perturbation such as associated, for example, with a merger of two comparable 
mass galaxies and their MBHs (a "major merger"). This assumption is consistent with the low mass of 
the MBH, and with the steep stellar cusp observed around it (§3.2), since major mergers are thought to 
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result in flattened stellar cores due to ejection of stars by the two merging MBHs (Milosavljevic et al 
2002). 

The Bahcall-Wolf solution applies to point particles. This assumption no longer holds very near the 
MBH, where the collision rate is high because of the very high stellar density and velocity. Stars on 
tight orbits around the MBH cannot survive for long, and so eventually most of the population there will 
consist of stars that are on very wide, marginally bound (parabolic) orbits, which spend only a small 
fraction of their time in the collisionally dominated region. These marginally bound stars have a flatter 
spatial distribution, of the form (e.g. Binney & Tremaine 1987, p. 551) 



^ocr _1/2 . (3.14) 

Near the MBH, at r<^ri, the potential is Keplerian to a good approximation, ifi(r) = Gm / V, and 
the orbits can be characterized by their semi-major axis a (or their period P) and their eccentricity e 
(0 < e < 1 for bound orbits), 



Gm 7 J a 3 

a = — , e 2 =l- — , P = 2ir\—, (m»M 7t ), (3.15) 
2e Gma V Gm 

where / is the magnitude of the specific angular momentum, such that the orbital energy is given by 

1 7 Gm 1 9 J 2 Gm ,„ „ 
£= v 2 + _ = v 2_ — + —, (3.16) 

2 r 2 2r L r 

where v r is the radial velocity. The stellar densities per dr, de and da that correspond to a DF of the 
form f(e)=Ae p (Eq. 3.12) are, respectively, (e.g. Schodel et al 2003, Eqs. A1-A8) 

"*« = (2-) 3/2 ^^T(Gm) 3 /^ A ,- 3 /^, 
1X5/2+/?) 

n*(e) = V2,Tr\Gm) 3 Ae p - 5/2 , 

n*(a) = 2 2 - p n\Gmfl 2 + p Aa l l 2 - p . (3.17) 
In the general spatially spherically symmetric case (anisotropic velocity field), the DF f(e,J) de- 
pends both on e and /. The distribution of specific energy and angular momentum, n(e,J) (stars per 
dedJ) is 

n*(e,7) = STT 2 Jf(e,J)P r (e,J), (3.18) 

where P r is the radial period. For Keplerian orbits P, = P and the distribution of specific energy and 
eccentricity (stars per dede) is (cf Cohn & Kulsrud 1978) 



«*(s,e) = 



2V2TT 3 (GM) 3 f(e,J)e~ 5/2 



(3.19) 



It then follows that for Keplerian orbits with isotropic velocities the normalized distribution of eccen- 
tricities n(e) (stars per de interval) is simply (cf Binney & Tremaine 1987) 



n*(e) = 2e. (3.20) 

Thus, isotropically distributed stars on Keplerian orbits tend to have high eccentricities, with a fraction 
l-e 2 having eccentricities >e (for example, 10% have e>0.95). 
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The velocity distribution is tied to the spatial density distribution through the Jeans Equation (the 
continuity equation of the stellar orbits in phase space in terms of the stellar density and velocity disper- 
sion) 



v 2 = d\nn± dlna 2 
a 2 dlnr d\nr ' 

where v c = \/ GM(< r)/r is the circular velocity and a steady state, isotropic, non-rotating system is as- 
sumed. The assumption of steady-state is justified because the dynamical timescale is much shorter than 
the relaxation timescale. The assumptions of approximate isotropy and non-rotation are observationally 
justified for the old stellar population and for the stars in the central 1" (§2.2). Since very near the MBH 
the stellar mass is negligible, M(<r)~m, the velocity dispersion is Keplerian, a 2 oc r -1 , and the density 
is a power-law, oc r~ a (here a =3/2+/?), it follows from the Jeans equation that 



„».( ' (3.22) 
\5/2 + pmJ r 

This result justifies the approximation that the velocity dispersion in a relaxed stellar system around a 
MBH is mass-independent, since cr^ changes by less than 10% over the entire mass range, in marked 
contrast to the a\j oc M~ x dependence of equipartition. 

Theoretical arguments indicate that the velocity field of a relaxed system very near a MBH should 
be approximately Maxwellian and isotropic locally (Quinlan, Hernquist & Sigurdsson 1995). This is 
indeed observed to be the case in the GC. The measured radial velocity distribution is well fitted by 
a Gaussian distribution (Genzel et al. 1996) and the measured velocities in the plane of the sky and 
along the line-of-sight are consistent with an isotropic velocity dispersion (Ghez et al 1998; Genzel et 
al 2000). The Jeans equation (Eq. 3.21) indicates that the steeper the cusp (larger p M ), the larger v c /a, 
and so the fraction /„ of unbound stars (those with v> v e = \/2v c ) is smaller 7 , 

/2 1"°° 
- / u 2 e- u2 ' 2 Au. (3.23) 

The fraction of unbound stars varies greatly across the range of interest because the integration over 
the exponential tail of the Maxwellian distribution depends sensitively on the lower limit. For example, 
in a flat core (n* ~ const.), /„(-3/2) ~ 0.6, whereas in a steep cusp (e.g. n^ocr -2 5 ), /„(1) — 0.07. 
Because unbound stars have wide orbits and spend most of their time far away from the MBH, the 
stellar population in a shallow cusp is well mixed and representative of the average population over a 
large volume. In contrast, the stellar population in a steep cusp is localized and can therefore develop 
and maintain properties that differ from those of the general population. 



3.2 The stellar cusp in the Galactic center 

The basic premise of the attempt to quantify the global properties of the stellar distribution in the inner 
GC is that there is an underlying quasi-steady state distribution that includes most of the mass and most 
of the stars, and which is traced by the old faint stars. In contrast, the young stars and bright giants, 
which dominate the light, are less reliable tracers because they could be far from steady-state, or subject 

7 The stars are bound to the GC as a whole (MBH and stars), but when they approach the MBH they move with 
v > v e in the two-body MBH-star system. 
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to statistical fluctuations due to their small numbers, or susceptible to environmental effects such as 
stellar collisions or photoionization (Sellgren et al. 1990). 

The proximity of the GC allows, in principle, a direct test of the theoretical predictions about the 
distribution of stars around a MBH. However, even in the GC the derivation of the 3D mass and number 
density distributions from the observations is hindered by several substantial obstacles. 

1. The star counts are inevitably incomplete. Only stars brighter than the detection threshold can 
be observed. To observe fainter stars, the instrumental photometric sensitivity must be improved, 
and with it the angular resolution in order to resolve individual stars in the crowded field near 
the MBH and avoid the source confusion problem. The dynamic range of the detectors must 
also be increased, otherwise the bright stars in the frame will saturate the detectors and out-shine 
their faint neighbors. These instrumental requirements can rarely all be satisfied simultaneously. 
The degree of incompleteness can be estimated and compensated for in the statistical sense, by 
simulating the detection procedure with mock data. 

2. Obscuration ("extinction") by a non-uniform distribution of interstellar dust can bias the derived 
stellar distribution. To compensate for such effects, it is necessary to map the dust across the 
field, for example by comparing the observed stellar colors with their assumed intrinsic ones. 
The analysis of the extinction gradient across the central parsec suggests that any bias due to dust 
extinction is unlikely to be large (Alexander 1999; but see strong gradients in extinction maps by 
Scoville et al. 2003). 

3. It is crucial to separate the old and young stars in order make a meaningful comparison to stellar 
dynamical models, since the two populations are dynamically distinct (§2.2). The phase space 
distribution of the young stars reflects the initial conditions at birth and their still unknown for- 
mation mechanism (§7), whereas the distribution of the old stars has had time to adjust to the 
presence of the MBH and is probably dynamically relaxed. The stellar type can be determined, 
at least for the brighter stars, by their spectra or infrared colors. 

4. The observed 2D star counts have to be translated to 3D densities. This can be done by first mak- 
ing assumptions about the symmetries of the DF (e.g. spherical symmetry and isotropic velocity; 
Genzel et al 1996), deprojecting the observed surface density and line of sight velocity dispersion 
(Eq. 4.2) and relating the deprojected quantities to the enclosed mass via the Jeans equation (Eq. 
3.21). Alternatively, the star counts can be fitted to a projected 3D model distribution with a small 
number of free parameters (e.g. a broken power-law cusp), for example by maximal likelihood 
criteria (Alexander 1999; Genzel et al 2003b; see also Merritt & Tremblay 1994 for a more so- 
phisticated approach, yet to be applied to GC data). This method does not take into account the 
velocity information, which in any case is not available for all the stars. 

5. Only the luminous stars are directly observed. However, most of the mass is expected to be 
in the very faint low mass stars, compact remnants or other forms of dark mass, if such exist 
(§4.2). In order to obtain the mass density p(r), it is necessary either to measure the diffuse mass 
dynamically, which is very difficult to do in the center, where the potential is dominated by the 
MBH (§4.2.3), or else assume a stellar number-to-mass ratio, N+/M, in order to translate the 
stellar number density n*(r) to mass density. Such a procedure is very uncertain, since N+/M 
probably depends strongly on the distance from the MBH due to mass segregation. 

The light distribution in the GC approximately follows a ~ r~ 2 power-law from as far away as 
~ 100 pc, where it cannot be directly related to the central super-massive BH, to ~ 1 pc from Sgr A* 
(Serabyn & Morris 1996; §2). Early attempts to determine the stellar distribution inside the central 
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Figure 3.1: Measurements of the surface density distribution of the stellar counts around the Galactic 
MBH (Genzel et al 2003b). The observed surface density density corresponds to a space number density 
of n*ocr~ a with a = l.OzLO.l inside 0.4 pc, and a = 2.0±0.1 outside. (Reprinted with permission from 
the Astrophysical Journal). 

parsec concentrated on identifying the core radius r c , where the surface density falls to half its central 
value, by measuring the surface density of the brightest stars (the only ones observable at the time) or 
of the total surface brightness. These studies yielded conflicting results, from r c ~0.05 pc (Allen 1983; 
Scoville 2004) to ~ 1 pc (Rieke & Rieke 1988). It appears that the differences between the various 
estimates depend on whether the light distribution or the number counts are used and whether or not the 
bright stars are included in the fits. A consensus has yet to emerge on this question. 

Generally, star counts follow the mass more faithfully than the total light distribution, which is 
completely dominated by the light of giants, even though they are only a small fraction of the population 
both by number and mass. Analyses of the stellar distribution based on star counts (Alexander 1999; 
Genzel et al 2003b) indicate that a stellar cusp rises monotonously all the way into the central ~0.1" 
(0.004 pc) (Fig. 3.1). The number density was determined by fitting the star-counts to a 3-parameter 
broken power-law (break radius and the inner and outer exponents) using maximal likelihood and taking 
into account completeness corrections. The number density was then converted to a mass density by 
assuming a constant N+/M ratio and normalizing to the enclosed stellar mass M+(< r)=M(< r)—m far 
from the MBH. The total enclosed mass M(< r) was measured dynamically on the < 10 pc scale, and 
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the MBH mass m was measured dynamically close to the center (§4.1). The resulting mass model is 
(Genzel et al 2003b; assuming /?o = 8 kpc) 

p(r)=po(^) , po = 1.2xlO 6 M pc- 3 , r =0.4pc, a=| ^±0.'! rSro ' ° M) 
This translates to an enclosed stellar mass in the inner cusp (r < ro) of 

M.«r) = 6xlO 5 (r/r o ) L6 M , (3.25) 

and extremely high densities near the MBH, p(0.04pc) = 3 x 10 7 M pc -3 and p(0.004pc) = 8 x 10 8 M Q pc" 3 
Note however that the total mass enclosed in the high density regions is very small, only ~ 8000 M Q 
inside 0.04 pc and ~200M Q inside 0.004 pc. These estimates of extreme densities are consistent with 
the observed gradual depletion of bright giants with decreasing distance to the center, which is natu- 
rally explained by a combination of projection effects and collisional destruction in a steep density cusp 
(§3.4). 

The observed power-law is consistent with the Bahcall-Wolf solution for the less massive stars in 
a relaxed population. This may hint that the mass density inside ro is dominated by dark objects more 
massive than the < 10M Q stars which constitute the majority of the observed population. If this is 
indeed the case, as is further explored in §3.3, then N+/M must rise toward the center, contrary to what 
was assumed, and so the mass model (Eq. 3.24) should be interpreted as representing the run of the 
density in the luminous stars only. The total mass density distribution may have a steeper slope than the 
observed number density, and the mass density close to the center may be yet higher than the estimates 
above 8 . 



3.3 Mass segregation 

It is quite likely that the stars observed near the MBH constitute only a small fraction of the total 
extended mass there. Mass segregation is expected to drive a flow of NSs and SBHs to the center 
since these compact remnants are more massive than the mean stellar mass in the unsegregated parent 
population (table 2.1) and since, unlike the short-lived massive MS stars, the remnants are effectively 
eternal and have time to segregate significantly (Moms 1993; Miralda-Escude 2000). SBHs, which are 
the most massive stellar remnants, are expected to be the most centrally concentrated. 

The farther from the center a compact object is born, the longer it takes it to sink in. The process 
can be described in terms of dynamical friction, the drag experienced by a massive object M, as it 
moves through a background of lighter masses (the drag is due to the work done by the massive object 
in focusing the lighter masses across and behind it). When the velocity DF is Maxwellian, and in the 
limit M. ^> {M^), the deceleration due to dynamical friction, v'df, depends only on the total mass density 
of the background stellar population and not on the individual masses (Chandrasekhar 1943; Binney & 
Tremaine 1987, Eqs. 7-17, 7-23), 



Vdf = - 



4vrlogAG 2 M.p(r) 



2X , ' 

erf(X)-^exp(-X 2 ) 

7T 



(3.26) 



where X = v/\^2a (Note that the dynamical friction timescale for v^a, t&= |v/v'df| ~f r (M*) /M., is the 
mass segregation timescale, Eq. 3.9). To estimate the size of the "collection basin" from which compact 
objects can reach the center over the lifetime of the Galaxy tji ~ 10 Gyr, consider the orbital evolution of 



8 For a given, dynamically determined enclosed mass Mo within ro, a power-law density p(r) 
[47tMq /(3-a)rg] (r/ro)~ a increases with increasing a inside r Q = roexp[l /(a-3)] and decreases outside r a . 
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a compact remnant of mass M. that is initially on a circular orbit with velocity v c at some radius r from 
the center. The dynamical friction drag exerts a torque L=v^r that causes the orbit to decay and the mass 
to spiral inward. It is convenient to approximate here the central density distribution as an isothermal 
cusp, p(r) = vl/4irGr 2 , whose circular velocity v c = Via is independent of the radius, so that L = v c f = 
v'dfr . This assumption does not introduce a large error since the stellar density distribution outside the 
central regions is approximately isothermal (§3.2), and most of the time spent by the inspiraling object 
is at large radii. The solution of the torque equation relates the time available for reaching the center, 
fdf , to the maximal initial distance max r dt -, 

f log A GM^\ 1/2 

maxr df = I y-jy— - — I ~5pc, (3.27) 

for tdf = tH, an SBH of mass M, = 10M Q , and for log A ~ 10 and v c ~ 150kms _1 on the < 10 pc scale 
(Kent 1992). This is in very good agreement with estimates based on more realistic density profiles for 
the GC (Morris 1993; Miralda-Escude 2000). 

It is difficult to estimate reliably the number of SBHs, N., that have sunk to the center in the course 
of the Galaxy's lifetime, or the characteristic spatial extent r, of the central SBH cluster. Such an 
estimate depends on various uncertain quantities and properties: the IMF, the star formation history of 
the GC, the relation between the initial stellar mass and the final SBH mass (the SBH mass function) and 
the dynamical conditions in the GC. In addition, the magnitude of the SBH natal "kick" velocity could 
affect the segregation efficiency, since faster SBHs sink more slowly (Eq. 3.27). Estimates lie in the 
range N,=lAx 10 4 within r.=0.7 pc (for M.=7M SBHs; Miralda-Escude & Gould 2000) to between 
AT.=3.2xl0 4 and/V.=7.4xl0 5 within r.=0.8 pc (forM. = 10M Q SBHs; Morris 1993). However, values 
of ./V, > 10 5 are inconsistent with the current dynamical measurements of the enclosed stellar mass on 
these scales (§3.2). Inside r. the SBHs are expected to be distributed in a r" 7 / 4 cusp (§3.1.3). 

Irrespective of the uncertainty in N m , this is an extreme degree of concentration. The inner parsec 
of the Galaxy may contain up to ~ 10~ 2 of all Galactic SBHs, in ~ 10~ 10 of the Galactic volume. Such 
a degree of mass segregation will significantly affect the run of the N+/M ratio with radius, which may 
fall by an order of magnitude from r~ 10 pc to r~0.01 pc (M. Freitag, private, comm.). 

It is possible to obtain a rough upper limit on the central concentration of SBHs from purely dy- 
namical considerations. The more concentrated the SBH cluster, the higher the rate at which two-body 
scattering will deflect SBHs into event horizon-crossing orbits. The "drain limit" is set by considering 
how many SBHs can be packed inside a radius r,, in steady state, so that the number scattered into the 
MBH over the age of the Galaxy equals the number enclosed (Alexander & Livio 2004; §6.1.1) 

1 dN. 1 1 

< — . (3.28) 



N. dt iog(2 ^J7 s )t r t H 
Equation (3.28) translates to the upper bound (solved numerically), 

21og(2v^7^) fm\ 2 P(r.) 
maxA r .(< r.)~ — — — f — — — . (3.29) 



31og(0.4maxAr.) \M.J t H 

This limit is particularly useful close to the MBH, where the enclosed stellar mass estimates become 
unreliable due to the uncertainties in the determination of the MBH mass and of N+/M. The drain limit 
for 1OM SBHs is ~ 10 4 inside r.=0.1 pc. 

The central cluster of massive compact remnants is not easy to detect. Three possible approaches are 
discussed in the literature: (1) direct detection of the compact objects; (2) detection of mass segregation 
by the central depletion of light "test particles"; (3) indirect detection of the compact objects by their 
effects on the environment. The evidence for mass segregation around the Galactic MBH is at present 
only circumstantial. 
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Direct dynamical detection of the smoothed potential of the extended mass by deviations from 
Keplerian motion requires very high astrometric precision over multiple orbital periods. The current 
upper limit is still ~30 higher than predicted by the drain limit (Mouawad et al 2005; §4.2.3). Orbital 
perturbations by occasional close interactions with individual compact remnants may be detectable with 
future very large telescopes (s> 10m) (Weinberg, Milosavljevic & Ghez 2005). The gravitational lensing 
properties of the MBH with respect to background stars will be modified by the SBHs orbiting it, in 
particular when the light rays of a lensed image pass near one of the SBHs (Gould & Loeb 1992; 
Miralda-Escude & Gould 2000; §5.2). However, the detection of background stars lensed by an MBH- 
SBH "binary" requires a photometric sensitivity ~ 100 times larger (K ~ 23 mag) than is presently 
available and even then the events are predicted to be rare (~0.06yr _1 ) (Chaname, Gould & Miralda- 
Escude 2001). Accretion of the dilute interstellar gas by fast moving compact remnants in the inner ~0. 1 
pc will probably not lead to significant X-ray emission (Haller et al 1996; Pessah & Melia 2003). On 
the larger, 1 pc scale there is an observed over-abundance of transient X-ray sources (4 sources observed 
compared to 0.2 expected on average). This is interpreted as evidence for a high central concentration of 
NS and SBHs, of which some are accreting from a binary companion acquired in the course of a 3-body 
exchange encounter with a stellar binary (Muno et al 2005). 

Chaname & Gould (2002) propose using the angular distribution of millisecond pulsars around 
the MBH for detecting mass segregation. Millisecond pulsars, being NSs, are lighter than the average 
stellar mass in the SBH-dominated center. Since they are older than the relaxation time (their ages are 
estimated from the spin-down timescale, P/2P~few x 10 9 yr), they should be depleted in the center due 
to mass segregation. However, the radio pulses are hard to detect because they are strongly scattered 
by the interstellar plasma in the GC. Lazio et al (2003) report finding only 10 candidate pulsars in 
the central ~ 200 pc of the GC. Deeper surveys at higher frequencies will be required to avoid the 
scattering and to significantly improve these statistics. Another class of potential light test particles are 
the long-lived, low mass MS progenitors (M* <2M Q ) of the He burning "horizontal branch / red clump" 
giants (Alexander et al., in preparation). While the low mass progenitors are themselves too faint to be 
observed by current instruments, their giant progeny are luminous enough (K ~ 16) and are readily 
detected. Observations show that these giants gradually disappear from the population with decreasing 
projected distance from the MBH (Genzel et al 2003b). One possible interpretation is that this is due 
to mass segregation operating on the progenitors. However detailed calculations are still lacking, and 
other explanations, such as collisional destruction (§3.4), have yet to be ruled out. 

Finally, the compact remnants could be revealed by their effect on the relatively few burning stars 
in the center. Collisional destruction by SBHs, perhaps accompanied by the creation of exotic collision 
products (e.g. Thorne-Zytkow objects) is a possibility that has not yet been explored in detail. Strong, 
non-destructive interactions between SBHs and MS stars of similar mass may be responsible for the 
capture of the young B-stars observed there (Alexander & Livio 2004; §7.2.4). 



3.4 Stellar Collisions 

With central densities above 10 8 M Q pc~ 3 , the inner cusp is the densest environment in the Galaxy, 
exceeding even the core density of massive globular clusters by up to a factor of 100. In this extreme 
density the mean time between physical stellar collisions, t co \\ (Eq. 3. 1 1) is shorter than the age of system 
which is also roughly the lifespan of a ~ 1M Q star. Since the probability to avoid a collision is 
exp(-f#/? co ii) <C 1, the dynamics and stellar population in the high density inner cusp are completely 
dominated by collisions. The GC is essentially a naturally occurring "stellar collider". 
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3.4.1 Collisional destruction of giants 

Images of the inner arcseconds of the GC clearly reveal a central cavity in the distribution of luminous 
late-type red giants (Fig. 2.2). The depletion is gradual (Fig. 3.2). Inside ~5" (0.2 pc projected distance) 
there are no K < 10 mag late-type giants, and inside ~ 1.5" (0.06 pc projected distance) there are even 
no K< 12 mag late-type giants. Apart for the luminous, short-lived He stars (§2.2), the maximal stellar 
luminosity gradually decreases in the inner 2". 

A possible explanation for this trend, selective obscuration of the center by a centrally concentrated 
distribution of dust, is not supported by an analysis of the measured extinction across the inner parsec 
(Alexander 1999). Alternatively, stellar collisions could be responsible for the depletion of the luminous 
late-type giants (Lacy, Townes & Hollenbachl982; Phinney 1989; Genzel et al. 1996) or the disappear- 
ance of CO absorption — a distinctive spectral signature of red giants — in the integrated light (Sellgren 
1990; Figer et al 2000). Luminous red giants have very large extended envelopes, and therefore a large 
cross-section for collisions with other stars. A strong collision may strip the giant's envelope leaving 
behind an almost bare burning core. This will drastically lower its IR luminosity because the IR spec- 
tral range lies in the Raleigh- Jeans part of the stellar blackbody spectrum (for effective temperature 
>4000K), and so the IR luminosity scales as L IR oc R^T+ whereas the total stellar luminosity scales 
as L* oc R\T^. Consider a collision that disperses the envelope of a ~ 100 Rq red supergiant and leaves 
a ~ IRq burning core. Studies of rapid mass transfer in close binaries (e.g. Podsiadlowski, Rappaport 
& Pfahl 2002) indicate that a stripped core retains its nuclear luminosity. The factor 10 4 decrease in sur- 
face area must then be compensated by a factor 10 increase in the effective temperature, which translate 
to a factor 10 3 (7.5 mag) decrease in the IR luminosity. 

Despite its low surface gravity, a red giant is quite robust against collisions. The dynamical timescale 
in the envelope is of the order of a year and the thermal timescale is of the order of several decades, 
so a collision in which a star "punches a hole" in the giant's envelope will probably not have a lasting 
effect. However, hydrodynamical studies of collisions between dwarf stars and giants indicate that the 
envelope can be completely disrupted when the impact parameter is significantly smaller than the stellar 
radius. This happens either by the ejection the giant's core from the extended envelope, which then 
rapidly disperses, or by the formation of a common envelope binary, which stirs the envelope and leads 
to its evaporation on a time scale much shorter than the giant phase (Livne & Tuchman 1988; Davies & 
Benz 1991; Rasio & Shapiro 1990,1991; Bailey & Davies 1999). 3-body collisions between a binary 
and a giant are even more efficient in disrupting the envelope (Davies et al. 1998). However, it is hard 
to estimate the overall contribution of 3-body collisions because the binary fraction near the MBH is 
uncertain. Note that collisional destruction is less likely to occur in early-type blue supergiants, which 
do not have a well-separated compact core / extended envelope structure. 

Figure (3.2) compares the observed distribution of stars in the plane of projected distance p and 
magnitude K with theoretical contours of the expected mean number of surviving giants in a collisionally 
dominated cusp. The stellar spectral type is indicated, where such information is available. The early- 
type stars, which are too young to be dynamically relaxed, clearly stand out in the p-K plane. They 
are too short-lived to be affected by stellar collisions. In contrast, the old late-type giants are relaxed 
and susceptible to collisions. The theoretical calculations estimate the collisions rate (Eq. 3.11) and 
survival probability by assuming solar mass projectiles, use cross-sections for collisional destruction 
calibrated by hydrodynamical simulations, and take into account the stellar DF (Eqs. 3.24, 3.22) and 
mass function (Eq. 2.1) of the GC and the evolution of the stellar radius and i^-band luminosity along 
the stellar tracks (Alexander 1999). The possible effects of mass segregation in enhancing the collision 
rate and efficiency are not taken into account. 

The results show that two distinct effects can cause the observed depletion. One is purely geomet- 
rical. The number of stars enclosed in an annulus of projected radius p and width dp scales as p dp 
in a n+ oc r~ a cusp. For cusps shallower than r~ 2 the number of stars decreases with p, and so does the 
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Figure 3.2: Evidence for the collisional destruction of the envelopes of bright red giants in a high den- 
sity stellar cusp around the Galactic MBH (updated version of figure 9 in Alexander 1999). The apparent 
stellar K-band magnitude of stars in the inner 5" is plotted against the projected angular distance from the 
black hole, p (VLT CONICA/NAOS data from Trippe, Ott et al. 2005, in prep.). The stars are sorted into 
spectral types (filled circles: late-type, crosses: early-type, open circles: unknown) according to their 
spectral features (when available) or according to their CO-band color m(CO)=m(2.29^m)-m(2.26^m) 
(from the Gemini science demonstration dataset), which is sensitive to the characteristic CO absorption 
feature in late-type giants (Genzel et al. 2003b). Most of the low luminosity stars without spectral 
identifications are probably also old late-type stars. The three dashed contour lines represent model 
predictions for the decrease in surface density of bright stars due to projection. The three solid lines 
include also the effects of collisional destruction in a high density n*=1.2xl0 6 (p/10T L4 stellar cusp 
(§3.2). The model predicts, on average, 1.5 (top contour), 1.0 (central contour), and 0.5 (bottom contour) 
dynamically relaxed (old) stars per 0.25 arcsecond bin that are brighter than the contour level. 



probability of observing the rarer, luminous giants. The second effect is dynamical. Once the stellar 
density rises above ~5x 10 7 M Q pc~ 3 , collisions become effective. The additional collisional depletion, 
above that due to projection, improves the fit of the model to the data in the inner ~ 1.5". Thus, the ob- 
served gradual depletion of bright late-type giants toward the center cannot, in itself, be taken as proof 
of a high density cusp. However, taken together with the observed rise in the surface density (§3.2), 
both properties indicate that a steep, high density cusp provides a self-consistent interpretation of the 
data. 
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3.4.2 Tidal spin-up 

The Keplerian velocity in the inner parsec around the MBH is a substantial fraction of the escape velocity 
from a typical star, and even exceeds it at r<r co n (§3.1.2). Stars moving rapidly in the dense stellar 
cusp near the MBH will suffer over their lifetimes numerous high-velocity, close tidal encounters with 
other stars. Although such encounters transfer some energy and angular momentum from the hyperbolic 
orbit (e < relative to the other star) to the colliding stars, they rarely remove enough energy for tidal 
capture. This is in marked contrast to the situation in the high density cores of globular clusters, where 
the colliding stars are on nearly zero-energy orbits and close collisions can lead to the formation of tight 
binaries. The effects of hyperbolic encounters on the stars are mostly transient. The dynamical and 
thermal relaxation timescales are short compared to the stellar lifespan, and thus apart from some mass- 
loss in very close collisions, the star is largely unaffected. It is however more difficult for the star to 
shed the excess angular momentum since magnetic breaking (the torque applied to a star when the stellar 
wind resists being swept by the rotating stellar magnetic field) typically operates on timescales of the 
order of the stellar lifespan (Gray 1992). Fast rotation is therefore the longest lasting dynamical after- 
effect of a close encounter. Over time, the contributions from many, randomly oriented tidal encounters 
will lead to a "random walk" buildup of the stellar spin (Alexander & Kumar 2001). 

The possibility that stars in MBH cusps are rapid rotators may have interesting implications for their 
evolution and the interpretation of their observed properties (see review by Maeder & Meynet 2000). 
For example, rotationally induced mixing may reveal itself in the spectral line ratios, and rotation may 
be directly observed in the spectral line profiles (but high rotation will be suppressed by expansion in the 
giant phase). Detection of such signatures in the spectra of the observed giants can provide additional 
evidence for the existence of an underlying dense cusp of low-mass MS stars, which at present cannot 
be directly observed. 

Consider a a star of mass M+ and radius that undergoes a single tidal encounter with a perturbing 
mass m (here a star or a compact object). When the tidal deformations in the star are small, the energy 
taken from the orbit and invested in raising the tides can be described by a linear multipole expansion 
in the periapse distance (distance of closest approach) r p =a(\-e) (Press & Teukolsky 1977), 



AE t fm\ 2 ^ T,(v,e) T 2 (3 _ 3Q) 



1=2 

where E+^GMl/R* is the stellar binding energy (up to a factor of order unity depending on the stellar 
structure), 7} are the dimensionless tidal coupling coefficients, which depend on the eccentricity e (e > 1 
for an unbound orbit, Eq. 3.15) and on the dimensionless periapse crossing timescale 77 = 7^/7-* where 

T p = Jr p /G{M*+m) is the periapse crossing timescale and = \/Rl/ GM± is the stellar dynamical 
timescale. The tidal coefficients 7} can be calculated numerically for any given stellar model and orbit 
(e.g. Alexander & Kumar 2001). For many applications it suffices to use the lowest order multipole, T 2 . 
The angular velocity at periapse is 

p p 

where is the relative velocity at infinity. Generally, T 2 peaks when Q, p matches the frequency of 
the fundamental mode of the star, Q p ~ 0(t~ 1 ). The higher the velocity of the encounter (the larger e, 
Eq. 3.15), the larger the value of r p (larger r]) for resonance with the fundamental mode. The distance 
dependence of AE, (to lowest order) can be understood by noting that the energy in the tidal oscillations 
is proportional to the square of the amplitude of the tidal elongation, AR+/R+ ~ (R^/rp) (m/M*). 
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The orbital energy, AE t , and angular momentum, A/ f , that are transferred to the star in an impulsive 
tidal encounter (such as a hyperbolic one) are related by (e.g. Kumar & Quataert 1998) 



For rigid body rotation 9 , the change in the angular velocity of the star due to the tidal interaction is 
AQ t = AJ t /I ,where / is the star's moment of inertia (assumed here constant), so that 



where Q/, = \/GM±/R\ is the centrifugal breakup angular velocity, where the star sheds mass from its 
equator. 

Equations (3.30) and (3.33) apply only for soft encounters (r p > 2R+), where the linear expansion 
holds. Hydrodynamic simulation can then be used to extend Eq. (3.33) to strong encounters. Simula- 
tions show that as r p decreases, AO f first grows faster then predicted by linear theory, but then as r p 
approaches R+, the formal divergence of AQ t is truncated when mass is ejected from the tidally disturbed 
star and carries away the excess angular momentum (Alexander & Kumar 2001). Figure (3.3) shows the 
predicted mean rotation of solar mass stars in the GC cusp after 10 Gyr of stochastic tidal spin-up, as- 
suming inefficient magnetic breaking. The stellar spin reaches an r.m.s value of fi/fij, ~0.3 at the center 
(rotational velocity of 0.3V* ~ 130kms _1 ), about 60 times higher than is observed in solar type stars in 
normal environments. The spin-up effect falls only slowly with distance from the MBH (£l/Qb ~0.1 at 
0.3 pc) because the higher tidal coupling in slower collisions compensates for the lower collision rate 
there. Thus atypically high spin could be a characteristic of long-lived, low mass stars over a substantial 
fraction of the MBH's sphere of influence. The spin-up, being a random walk process, grows with time 
as f 1 / 2 and become significant only for long-lived stars. It is not relevant for short-lived massive MS 
stars or for the relatively short post-MS giant phase. For example, a 5M Q star with V* ~ 600 km s -1 and 
a MS lifespan of 4~ 10 8 yr will be spun up close to the MBH to only 0.3V* \Jt^jtu ~ 20km s" 1 (scaling 
from the calculated spin-up for a 1M Q star and assuming a similar spin-up efficiency). This is much 
less than the typical rotational velocity observed in such stars in normal environments, whose origin is 
related to their formation and not to tidal spin-up. 

3.4.3 Stellar mergers and exotic objects 

A possible outcome of a collision between two stars is a merger, where a new object, possibly more 
massive than either of its building blocks, is created. The high rate of stellar collisions near the MBH 
raises the question whether it is possible to create massive stars from lower mass stars by successive 
mergers, or exotic stars by the capture of a compact object inside a normal star (Morris 1993). This is 
of interest in view of the puzzling presence of young massive stars near the MBH, where normal modes 
of star formation are thought to be impossible (§7). The hypothesis that the very massive emission line 
stars in the central 0.4 pc are merger products is made less compelling by the fact that such stars are 
also found in the young massive stellar clusters in the central ~50pc, where the central densities are 
much smaller (§2.1). In addition, the collision rate is much too low at the relatively large distance of 

9 The timescale for angular momentum re-distribution due to convective transport in a red giant is ~ 1 yr (Zahn 
1989) and so rigid rotation is achieved on a timescale similar to that of the collision itself. The timescale for 
angular momentum re-distribution in radiative MS stars is not well known, although it is likely to be shorter than 
the stellar lifespan. 



AE t = AJ t n p . 



(3.32) 




(3.33) 
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Figure 3.3: The average spin-up of a solar type star by star-star tidal interactions over 10 Gyr as function 
of distance from the Galactic MBH (Alexander & Kumar 2001). An a= 1.5 density cusp is assumed. 
The rotation grows over time in a random walk fashion by repeated close passages. The angular rotation 
frequency 5Q = 5^l/Qb = ^ corresponds to rotation at the centrifugal break-up velocity. In addition to the 
total spin-up, the separate contributions from collisions with main sequence stars (MS), white dwarfs 
(WD), neutron stars (NS) and stellar black holes (BH) are shown. (Reprinted with permission from the 
Astrophysical Journal). 

the emission line stars (Eq. 3.11). However, it is still relevant to consider the merger scenario for the 
S-stars, which are both much lighter and closer to the MBH. 

A successful chain of mergers requires that the mean time between mergers be shorter than the 
lifespan of the merger product (usually assumed to be similar to the MS lifespan of a normal star of the 
same mass, and therefore becoming shorter as the merger chain progresses). The time constraint does 
not appear to be an insurmountable obstacle, if glancing collisions can lead to a merger (Genzel et al. 
2003b). However, it is very unlikely that the merging efficiency can be so high for such high velocity 
collisions. 

The first requirement for an efficient merger is that the impactor star (the smaller of the two with 
mass Mj) can be stopped in the target star (the larger of the two with mass M t ). This translates to the 
requirement that the target mass contained in the cylinder "punched out" by the impactor equals or 
exceeds the impactor mass, that is approximately Mi<RfR t (M t /Rf), where Rj and R t are the radii of the 
impactor and target. Thus, the stopping condition can be formulated as the requirement that the column 
density of the impactor be less or equal to that of the target, 



M i< M j _ 
Rj ~ Rj ' 



(3.34) 



Effective stopping does not in itself imply an efficient merger, since when Vqo > V e , the specific 
kinetic energy of the impactor is enough to unbind a unit mass of the target star. The total amount of 
mass that can be ejected from the target depends then also on the mass ratio between the impactor and 
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the target. The relative velocity at impact is > V e even when Vqo = 0, which is of the order of the sound 
speed in the star. Therefore such collisions will be supersonic, and the kinetic energy can be efficiently 
converted to heat by the shocks and drive the mass loss. The outcome of a penetrating encounter depends 
therefore not only on the mass ratio and the relative velocity, but also on the assumed stellar M*-to-/?* 
relation, and on the impact parameter (since the column density is a strongly varying function of position 
in the star). 

The actual situation revealed by numerical simulations is more complex than implied by the simpli- 
fied discussion above. For example, an impactor can be destroyed by heating even when it is not stopped 
by the target, and mergers can occur in slow enough collisions even without initial physical contact be- 
tween the two stars if they are tidally captured and form a decaying binary. The simulations generally 
show that the merger efficiency decreases rapidly when Vqo > V e , but the exact details depend strongly on 
the assumed stellar structure model. Lai et al. (1993) find that mergers can occur even when the relative 
velocity is as high as Voo ~ 1.5V e , provided that the collision is very close to head-on. However, Freitag 
& Benz (2005), using more realistic stellar models, find that mergers are much less efficient. Freitag, 
Giirkan & Rasio (2004) find no mergers in Monte-Carlo simulations of stellar collisions in a galactic 
nucleus with a MBH. Thus, massive star buildup by mergers very near the MBH does not appear likely. 

A related process is the capture of a compact object (NS of SBH) inside a star, and the possible 
formation of an exotic, accretion powered "star" (Thorne-Zyfkow object). The typical relative velocity 
at impact between the compact object and the star is v~ 0(Gm/r), which is much smaller than the 
escape velocity from the compact object, V e ~ 0(c). The collisional cross-section of the compact object 
is then dominated by the gravitational focusing term (see Eq. 3.1.2), which is much larger than its 
physical size. Therefore, the effective size of the impactor is the Bondi accretion radius, /?; ~ 0(GM,-/v 2 ) 
(Bondi 1952), and the column density ratio between the compact impactor and the target star is ~ 
(R t /r) 2 m 2 / (MjM t ). This ratio can be of order unity if the collision is atypically slow or occurs at a 
large enough distance (e.g. r~0.1pc for m = 3.5x 10 6 M Q and a collision between a M t = 1.4M Q NS 
and a M t = 1OM , R t = 4.5/?* star). A star powered by accretion on a NS may appear as a red giant or 
supergiant with peculiar photospheric abundances (Eich et al. 1989; Biehle 1991,1994), whereas one 
accreting on a SBH may look like a WR blue giant (Morris 1993). The evolution of such configurations 
in uncertain. 

4 Probing the dark mass with stellar dynamics 

Stars near the MBH are effectively test particles since their mass is negligible compared to that of the 
MBH, M*/m~ 10~ 7 -10~ 5 . As long as their trajectories do not take them too close to the MBH, where 
dissipative processes can affect their motion (§6), stars directly probe the gravitational potential near the 
MBH. This is in contrast to gas dynamics, where the possible influence of non-gravitational forces due 
to thermal pressure, radiation pressure or magnetic fields can substantially complicate the interpretation 
of dynamical data and limit its usefulness. The analysis of stellar orbits near the MBH, well within the 
radius of influence, is further simplified by the fact that the orbits are Keplerian to a good approximation 
(Eq. 3.15). 

The type of information about the dark mass that can be extracted from the data, the level of analy- 
sis possible, and the methods used, all depend on the nature of the observations (stellar position, proper 
motion, radial velocity, magnitude, colors, spectral classification), on their quality, on the properties of 
stars observed (proximity to the MBH, velocity, luminosity, age), and on the total duration of the moni- 
toring. Two approaches can be used to weigh and locate the dark mass (§4.1). One is statistical, which 
typically requires less data for each star (and so was the first to be applied), but needs a meaningfully 
large sample and is limited by various necessary assumptions about the parent population. The second 
involves the derivation of partial or full orbital solutions for individual stars. In principle, a single well 
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sampled stellar orbit is all that is needed to obtain the mass and location of the dark mass (assuming 
the dark mass is fully enclosed in the orbit; §4.2). However, this is only feasible near the MBH, where 
the periods are short enough for accelerations to be measured. Farther away from the MBH (r>0.1pc, 
P> 1000 yr) there is no substitute for statistical mass estimators to measure the total enclosed mass (stars 
and MBH). 

Following initial attempts to probe the potential of the central dark mass by gas kinematics (Wollman 
et al. 1976; Lacy et al. 1980), the use of stellar kinematics progressed to include radial velocities 
(McGinn et al. 1989; Sellgren et al. 1990; Krabbe et al. 1995; Haller et al. 1996; Genzel et al. 
1997), stellar proper motions (Eckart & Genzel 1996; Ghez et al. 1998), and stellar accelerations (Ghez 
et al. 2000; Eckart et al. 2002). Each of these subsets of the full orbital data were used to place 
some constraints on the dark mass (§4.1.2). Finally, it is now possible, using high-precision infrared 
astrometry and spectroscopy, to follow the trajectories and radial velocity curves of individual stars as 
they orbit the MBH (e.g. Schodel et al. 2002; Ghez et al. 2003b; Eisenhauer et al. 2003,2005). The 
combined input from stellar orbits and from high precision radio measurements of the position and 
proper motion of Sgr A* have all but ruled out non-BH alternatives for the dark matter in the center of 
the GC (4.2). This has also made it possible to directly measure the distance to the GC (§4.1.3) and 
strongly constrain the presence of a secondary BH in the GC (§4.3). A very different class of stars that 
may provide information on the central dark mass are high-velocity stars, well on their way out of the 
Galaxy (§4.4). The discovery of such stars can prove that the dark mass must be extremely compact to 
eject them so energetically, and is therefore likely a MBH. 



4.1 Weighing and pinpointing the dark mass 

The primary parameters of the dark mass are its mass, size, location and velocity. An empirical determi- 
nation of their values has many implications. The mass and size fix the mean density, which determines 
whether the dark mass is indeed a MBH (§4.2). The formal density associated with a BH (mass over 
volume of a Euclidean sphere whose radius is the event horizon) is 

3 m _ 3(c 2 /2G) 3 _ / m \~ 2 f l.lxl0 3 gcnT 3 
PBH ~4^r 3_ 4Trm 2 _ V4xlO 6 M J * \ 1.7 x lO 25 M p C - 3 ' (-) 

The higher the lower bound on p/pbh, the more plausible is the assumption that it is a MBH, and the 
less likely are alternative non-singular configurations. The question of the position the dark mass can 
be decomposed into three separate issues. (1) Where is the dynamical center (the center of attraction) 
in the GC? (2) Where is the dynamical center relative to the center of mass of the GC? (3) Where is the 
dynamical center relative to the radio source Sgr A*? The three positions need not be the same. 

The dark mass will execute Brownian-like motion relative to the stellar cluster around it due to small 
gravitational "kicks" arising from Poisson fluctuations in the stellar number density (e.g. Chandrasekhar 
1944). The amplitude of the Brownian excursions away from the cluster's center of mass will depend 
on how massive the dark mass is (the more massive it is, the less it will move) and on the properties 
of the stellar system around it (§4.2.3). It should be emphasized that due to the long range nature of 
gravitational interactions, Brownian motion in a self-gravitating system is fundamentally different from 
that in a gas. The statistical mechanics of the gravitational Brownian motion of a massive object, which 
itself dominates the potential, are still a subject of active investigations. Neither is it obvious that Sgr A* 
must lie at the dynamical center. This is of course expected if the source is an accreting MBH, since 
the emission should arise inside ~ 10rs~ 10 cm. However, the extremely low luminosity of Sgr A*, 
now also detected in the X-ray (Baganoff et al. 2001) and IR (Genzel et al. 2003a), does not require a 
MBH. It could easily be powered by accretion on a stellar mass compact object. In that extreme case, 
Sgr A* would merely be a test particle in orbit around the dark mass. To explore this possibility, it is 
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Figure 4. 1 : Infrared spectra of the early-type (hot) S-cluster stars. Top: Infrared spectrum and orbit 
reconstruction of S2/S0-2 (Ghez et al. 2003b. Reprinted with permission from the Astrophysical Jour- 
nal). Bottom: Composite spectra of 5 bright stars (K < 15 mag, top spectra left and right) in the H 
and /T-bands and 6 faint stars (15 <K< 16 mag, bottom spectrum left) in the i^-band (Eisenhauer et al. 
2005. Reprinted with permission from the Astrophysical Journal). 



necessary first to pinpoint the location of the radio source Sgr A* in the IR frame, where the stars are 
observed and where the dynamical center is identified 10 . This was achieved by the discovery of 7 red 
giants and supergiants in the inner 15" whose extended mass-loss envelopes (7?~5/?^) emit radio maser 
radiation (coherent light emission analogous to laser light) from a molecular transition of SiO (Menten 
et al. 1997; Reid et al. 2003b). Since these giants are also bright IR sources, they can be observed in 
both radio and IR and used to align the two reference frames and correct for various distortions. This 
enabled the determination of the position of Sgr A* relative to the IR stellar observations to within 10 

10 The IR reference frame is defined by the average angular position (p) = E^p,/A^ of a statistically large 
sample of stars around SgrA*. While each star is moving, (p) is assumed to be nearly stationary and the relative 
motion between (p) and the central dark mass is assumed to be negligibly small (Eisenhauer et al. 2003). 
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mas (la) and established that SgrA* indeed coincides with the dynamical center as determined by the 
orbits (§4.1.2). The radio/IR alignment also played a crucial role in the detection of the low level IR 
accretion emission from SgrA* (Genzel et al. 2003a; Ghez et al. 2004). 

Finally, a precise measurement of the position of the MBH along the line of sight, the distance to 
the GC, is important for reconstructing the structure and dynamics of the Galaxy (§4.1.3). 



4.1.1 Statistical estimators 

Astrometric observations of stars can provide the 2D projected angular positions p(f,) at discrete times 
{tj}, which by differencing yield the angular velocity (proper motion) p(t) and angular acceleration 
(proper acceleration) p(t). Spectroscopic observations can measure the Doppler shift of spectral features 
and provide the line of sight velocity z(fj), and the line of sight acceleration z(t). The position z along 
the line of sight can not be measured directly. In practice, not all this information is available for all the 
stars. Spectroscopy requires more photons than photometry, and so faint stars may have proper motion 
information, but not radial velocities. In addition, a reliable measurement of the Doppler shift requires 
distinct spectral features, which are not always available in the IR range of the spectrum. Proper motion 
can be measured only when the astrometric errors are smaller than the proper displacement, and so slow 
stars far from the center may have radial velocity information, but not proper motion. The astrometry of 
stars near the dense center can be unreliable because of source confusion and because they can be "lost" 
among the other stars due to their high velocity and large displacement between consecutive observation 
times. 

The available stellar data can be described statistically in terms of the mean surface number den- 
sity £(p) and the projected velocity dispersion of the fth velocity component <rj;(p)= ([v;(p)— v,(p)] 2 ) 
(assuming a value of Ro for converting proper motion to physical velocity). The means are taken over 
binned projected areas. It is usually assumed, as is indicated by the observed E, that the density distri- 
bution is spherically symmetric. In that case the bins are annuli of projected radius p and width Ap. 

The Jeans equation is the most direct way of using the mean quantities to estimate the enclosed mass 
(e.g. Binney & Tremaine 1987 §4.2. Id). It is assumed here for simplicity that the velocity dispersion is 
isotropic (see analysis of the non-isotropic case by Genzel et al. 2000). The projected quantities can be 
inverted by the Abel integrals, 



The derived functions 11 n(r) and a 2 (r) are then used to estimate the total enclosed mass M(r) from 
the Jeans equation (Eq. 3.21). If there is a compact mass in the center, the enclosed mass curve will 
converge to m with decreasing radius, otherwise it should decrease to zero (Fig. 4.2). 



11 In practice, the discrete /?-bins that are used for estimating a 2 and E are often too coarse for a stable 
solution of Eq. (4.2). Instead, some simple parametrized models are assumed for n(r) and <r 2 (r), and the 
best fit parameters are obtained by the inverse equations Yi(p) = 2R^ ^ nr(r 2 - R\p 2 )~ x l 2 Ar and Y,(p)a 2 (p) = 




(4.2) 



2R 2 J™na 2 r(r 2 -R 2 P 2 )-i/ 2 dr. 
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Figure 4.2: Enclosed mass distribution in the GC (Schodel et al. 2003) for R = S kpc (Reid 1993). 
Filled circles: Orbital solutions. Filled triangles: Leonard-Merritt projected mass estimator (data from 
Schodel et al. 2003 and Ott et al. 2003). Open triangle: Bahcall-Tremaine mass estimator (Ghez et 
al. 1998). Filled rectangles: Jeans eq. solution for a parameterized anisotropic model (Genzel et al. 
2000). Open circles: Jeans eq. solution for a parameterized isotropic model using radial velocities of 
late-type stars (Genzel et al. 1996). Open rectangles: Non-parametric mass estimator for an isotropic 
model (Chakrabarty & Saha 2001). "G": Mass estimates from gas Doppler shifts (Genzel & Townes 
1987). Solid line: Best fit to all data for a model of a central point mass with a stellar cluster. Dash-dot 
line: Mass in the stellar cluster alone — this model clearly does not fit the data. Dashed line: Visible 
cluster with a Plummer model of a hypothetical very compact dark cluster — this model marginally fits 
the data. (Reprinted with permission from the Astrophysical Journal) 

When the data is too sparse to allow meaningful estimates of £ and a 2 , more elementary statistical 
estimators can be used that do not require binning or inversion from projected to 3D quantities. These 
can be obtained by taking moments of the Jeans equation, (r i M(<r)) for integer k, and averaging over 
the stellar density (the virial theorem is also derived this way, see Leonard & Merritt 1989). Bahcall & 
Tremaine (1981) show that the virial theorem mass estimator 

M YT =^(z 2 )/(l/R p) , (4.3) 
where (...) denotes the sample average, is very biased and has other undesirable statistical proper- 
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ties. They propose instead a better behaved "projected mass estimator" for spherically symmetric and 
isotropic systems of test particles orbiting a dominant point mass, which uses line of sight motions, 



Leonard & Merritt (1989) propose a related mass estimator for a spherically symmetric non-rotating 
system, which uses proper motions and is anisotropy-independent, 



where p» and p± are the proper motions parallel and transverse to p. This estimator is unbiased only 
when the sample covers the entire radial extent of the cluster and for cluster density profiles that have a 
finite mass. However, the stellar data in the GC do not satisfy either of these requirements (the observed 
power-law density, Eq. 3.24, can not be extrapolated to infinity as the total mass does not converge). 
The bias in Mlm then depends both on the anisotropy and the steepness of the cusp. This can be partially 
compensated for by correction factors (Genzel et al. 2000). 

Chakrabarty & Saha (2001) propose a non-parametric mass estimation algorithm for a spherical 
and isotropic system that can use both line of sight velocities and proper motions. The method uses 
maximum likelihood to find the energy-binned DF and radius-binned potential that best reproduce the 
observed data. 

A different approach for estimating the mass, the "Orbital Roulette" method, was proposed by 
Beloborodov & Levin (2004). The method assumes (as does the virial theorem) that at any given time 
the stars are at randomly distributed phases in their orbits. For any trial mass distribution, the orbits can 
be calculated and the phases tested for randomness. If the trial mass is too low, the stars will tend to have 
high velocities relative to the local escape velocity of the trial potential, and therefore will appear to be 
atypically close to periapse. Conversely, if the trial mass is too high, the stars will appear to be atypically 
close to apoapse (r a = a(l + e), the maximal distance along the orbit). The best-fit mass distribution and 
its confidence interval can be identified by the location and shape of the phase-randomness maximum. 
This method cannot be applied directly to the stellar data in the GC, since it requires all 6 phase space 
coordinates of the stars, whereas only 5 are available at best. However, in the special case of stars in the 
inner thin disk, z can be derived from the constraint on the angular momentum J. Initial results indicate 
that the Roulette method mass estimate is consistent with the higher values obtained by orbital mass 
estimates (§4.1.2) (A. Beloborodov and Y. Levin, private comm.). 

Table (4.1) compares the statistical mass estimates for the central dark mass that were obtained by 
different authors using various methods and different data sets. The results all lie in the range m ~ (2- 
3)xlO 6 M . 

4.1.2 Orbital estimators 

The first step toward orbital reconstruction is to detect accelerations (curvature) in the stellar motion. 
Measurements of acceleration can provide quite strong constraints on the properties of the dark mass 
even without full orbital information. The mean enclosed density of the dark mass can be directly 
obtained from projected quantities through Newton's law of gravity (for an assumed angular position of 
the MBH), 




(4.4) 




(4.5) 




3 m 3 f 3 p 



(4.6) 



4-7T r 3 AttG r AnG p 



where r and r are the true 3D position and acceleration. The last equality holds because r and r 
are parallel, and so are decreased in projection by the same factor. For example, the position and 
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Table 4.1: Statistical mass estimates" of the central dark mass (in 1O 6 M ) 



JE h NT h 



BT b LM b 



Distance 



Refs f 



2.5-3.2 z 1.7±0.6z 
2.5±0.4b 2.6±0.8p 



1.9±0.6z — 
2.6±0.8p — 
2.6±0.2p — 



<0.2pc 
<0.08 pc 

< 0.2 pc 

< 0.2pc 



[1] 
[2] 
[3] 
[4] 
[5] 
[5] 
[5] 



2.5±0.2p 
3.3 b 2.5-2.6 p 



3.1±0.3p 2.9±0.4p 



2.0±0.7 p 



< 0.004 pc 

< 0.005 pc 
<0.05 pc 




a Proper motion data is designated by (p), radial velocities by (z), both by (b). 
b JE: Jeans Equation, VT: Virial theorem, BT: Bahcall-Tremaine 
b LM: Leonard-Merritt, NP: Non-parametric 

c [1] Genzel et al. (1996), [2] Eckart & Genzel (1997), [3] Ghez et al. (1998) 
c [4] Genzel et al. (2000), [5] Chakrabarty & Saha (2001) 



acceleration of the star S2 in 1997/1998 (/? = 150mas, p = 5.4masyr~ 2 , Ghez et al. 2000; Eckart et 
al. 2002) translate to a lower bound on the density of /)2D = 2x 10 12 M Q pc" 3 or p/pbh = 3x 10~ 14 for 
m=4x 10 Mq. Furthermore, since the projected acceleration vectors point to the projected position of 
the center of acceleration, they can locate the MBH in the IR frame (Ghez et al. 2000; Eckart et al. 
2002) (Fig. 4.3). 

When the inclination of the orbit is not known, the acceleration provides a lower limit on the MBH 
mass (for an assumed value of Rq), 

m > R 3 pp 2 /G = mcos 3 i , (4.7) 

where i is the inclination angle (see footnote 12). The measured accelerations of a few stars in the inner 
~0.5" give a lower bound of m > 1O 6 M (Ghez et al 2000; Eckart et al. 2002; Schodel et al. 2003). 

There is a close similarity between the S -stars orbiting the MBH and a planetary system (for ex- 
ample, M+/m~MEanh/MQ and Pjupiter^S Ar^S-Ppiuto)- The potential of the MBH completely dominates 
the dynamics, and perturbations by other stars can be neglected over the typical time span of the moni- 
toring campaigns, ~O(l0yr). Classical astronomical techniques for deriving the parameters of binary 
star systems from astrometric and spectroscopic observations can be directly applied, with only a few 
modifications (Salim & Gould 1999). 

When considered separately, each of the stars forms a binary system with the MBH. A binary orbit 
is described by 14 parameters: 6 phase space parameters for the star and the MBH each, and the two 
masses 12 . The mass of the star and the reflex motion of the MBH can be neglected because of the very 
large mass ratio. However, any constant velocity the MBH may have must be taken into account. This 

12 It is customary to describe a binary orbit in the reduced mass frame by 6 parameters (apart for the reduced 
mass): 2 parameters that determine the shape of the Keplerian ellipse, the semi-major axis a and the eccentricity 
e (Eq. 3.15), 1 parameter that determines the phase of the orbit, the time of periapse passage to, 2 angles that 
determine the orientation of the orbital plane in space, the inclination angle ; (the angle between the line of 
sight and the normal to the orbital plane) and the longitude of the ascending node SI (the position angle of the 
intersection between the plane of the sky and the orbital plane, the line of nodes), and 1 angle that determines the 
orientation of the ellipse in the orbital plane, the argument of pericenter u> (the angle between the line of nodes 
and the semi-major axis). 
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Figure 4.3: Pinpointing the position of the MBH in the IR frame by projected acceleration vectors 
(Ghez et al. 2000). The la error wedges on the acceleration of 3 stars overlap in a region ~0.05" from 
the nominal position of Sgr A* (at 0,0) as was initially derived from the alignment of the radio and IR 
grids by maser giants (Menten et al. 1997) (Reprinted with permission from Nature). 

is in contrast with the case of visual stellar binaries (binaries where both stars are resolved), where the 
proper motion of the binary's center of mass is irrelevant since the relative position between the two stars 
is directly observed. Until recently the MBH was not detected in the IR, where the stars are observed, 
and so the stellar orbits were measured against the IR grid, and not relative to the position of the MBH 
itself. Now that Sgr A* is also observed in the IR 13 (Genzel et al. 2003a; Ghez et al. 2004), the problem 
of solving the orbit will become easier. The results discussed here are still for the case of an invisible 
MBH. 

In principle, one well sampled stellar orbit with precise astrometric measurements is enough to 
solve all the parameters of the system apart for two degeneracies: (1) Only the combination m/R^ is 
fixed by the solution, but not m and Rq separately (cf Eq. 3. 15). (2) There remains a twofold degeneracy 
in the sign of the radial velocity. These degeneracies can be resolved with radial velocity information, 
in which case the MBH mass and the distance to the Galactic center can be determined (§4.1.3). The 
visual binary method for deriving distances utilizes the fact that the star's radial velocity z is measured 
via the Doppler shift of the stellar spectral features in terms of an absolute velocity, whereas the proper 
motion p is measured in terms of an angular velocity. The two are tied together by the orbital solution, 
rsoi(0 = (^oPjZX thereby yielding the distance to the binary. 

At present, the available orbital data for any single star is not sufficiently accurate to fully constrain 

13 Sgr A* is today easily detected in its flaring state, less so in its quiescent state, where it is only visible in part 
of the exposures. It is quite likely that with some improvements in the photometric sensitivity, Sgr A* will always 
be observable. 
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Figure 4.4: The projected orbits (left) and the radial velocity curve (right) of the 6 S-stars used by 
Eisenhauer et al. (2005) to solve for the mass, position and distance of the MBH. The points with error 
bars are the observed astrometric and spectroscopic data and the solid curves are the best fit multi-orbit 
solutions. (Reprinted with permission from the Astrophysical Journal). 



the solution. Tighter constraints can be obtained by solving simultaneously several orbits (Salim & 
Gould 1999). Since all the stars are orbiting a common point, the MBH, the ratio of data points to free 
parameters is much improved and this can significantly reduce the errors on the derived parameters (Fig. 
4.4). The best fit orbital parameters and their errors are then obtained by a non-linear least squares fit to 
the data. 

Full 3D orbital solutions provide yet tighter constraints on the density of the central dark mass than 
Eq. (4.6). If the mass is fully contained in the orbit, as suggested by the fact that the orbit is consistent 
with Keplerian motion (§4.2), the mean density is 



5 m 

P3D=4^- (4-8) 

The smallest periapse discovered to date, r p =0.22mpc for the star SO- 16, raises the lower bound on the 
mean density to p iD = 8 x 1 1 6 M Q pc" 3 , or p/pbh = 1 0" 9 f or m = 4 x 1 6 M Q (Ghez et al. 2005) . However, 
the tightest constraints on the compactness of the dark mass, p > 7 x 10 21 M Q pc" 3 , or pj pbh = 10 -4 come 
not from stellar orbits, but from radio observations (Reid & Brunthaler 2004). These yield a lower 
limit on the mass associated with SgrA* (>4x 10 5 M Q ) based on the small peculiar velocity (< lkms -1 
perpendicular to the Galactic plane), and an upper limit on the extent of the dark mass (<10 13 cm) based 
on the assumption that it is contained within the millimeter radio emitting region (Doeleman et al. 2001; 
Bower et al. 2004). 

Table (4.2) compares the orbital estimates for the mass, position and velocity of the dark mass that 
were obtained by different authors. The mass estimates all lie in the range m ~ (3^1)xlO 6 M , which 
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Table 4.2: Orbital mass, position and velocity estimates for the central dark mass 



m a 


N* b 




Ap c 


P 


Fixed parameters 


Refs e 


(1O 6 M ) 




(kpc) 


(mas) 


(mas yr" 1 ) 


for MBH d 




3.7±1.5p 


1 








p=z'=0, p,R 


[1] 


3.6±0.4p 


6 




1.5 




p=z=0, R 


[2] 


4.1±0.6b 


1 




2.4 




p=z=0, R 


[3] 


3.59±0.59b 


1 


7.94±0.42 


1.7 




p=z=0 


[4] 


3.7±0.2p 


7 




1.3 


1.5±0.5 


z=0, R 


[5] 


3.61±0.32b 


6 


7.62±0.32 


1.7 




p=z=0 


[6] 



a Proper motion data is designated by (p), both proper motion and radial velocities by (b). 
b Number of stars in orbital fit. 

c The uncertainty in the position of the center of acceleration in the IR grid, Ap=(Ap%+Apy) 1 ' 2 . 
d When fixed, Rq is assumed to be 8 kpc. 

e [1] Schodel et al. (2002), [2] Ghez et al. (2003a), [3] Ghez et al. (2003b) 

e [4] Eisenhauer et al. (2003), [5] Ghez et al. (2005), [6] Eisenhauer et al. (2005) 



is systematically higher than those derived from the statistical estimators m ~ (2-3)xlO 6 M , (§4.1.1, 
table 4.1). The discrepancy with the orbital mass estimates is most marked for the low value from the 
non-parametric method (Chakrabarty & Saha 2001). This is puzzling, since non-parametric estimators 
are expected be robust, but so are orbital mass estimators, which unlike the statistical estimators do not 
require assumptions about the properties of the stellar cluster. Orbital mass estimates are expected, if 
anything, to be slightly lower than the statistical mass estimates because they probe the potential closest 
to the MBH (r<0.02 pc), where the mass fraction of any extended mass (e.g. stars, compact objects) 
should be negligible. The systematic differences cannot be attributed to uncertainties in Rq, since all 
methods assume R = 8 kpc. This discrepancy has not been studied in detail yet, and the reasons for it 
remain unclear. 

4.1.3 Distance to the Galactic Center 

The Sun-GC distance, Rq, is a key parameter in the reconstruction of Galactic structure and dynamics as 
observed from our vantage point at the Solar system. For example, the current uncertainties in the values 
Rq and the local Galactic rotation speed ©o preclude the derivation of useful limits on the spheroidal 
axes ratios of the dark matter halo around the Galaxy, which could provide information on the nature 
of the dark matter (Oiling & Merrifield 2001; see discussion in Weinberg, Milosavljevic & Ghez 2005). 
The distance to the GC is also an important (though not the only) lower rung in the extragalactic distance 
scale, useful for calibrating standard candles such as RR Lyrae stars, Cepheids and red clump giants. 
Until recently, the adopted distance to the GC, 7?o = 8-0±0.5kpc (Reid 1993) was based on the average 
of various primary (geometric), secondary (standard candles) and tertiary (theoretical models) distance 
indicators (see brief summary in Eisenhauer et al. 2003). Geometric distance indicators are in principle 
the most reliable, the assumptions about physical properties of the system used are well founded. 

An early attempt to use stellar observations in the GC to determine the distance to the GC was based 
on a general relation that holds in a spherical system and relates the line of sight velocity to the projected 
velocity through Rq (e.g. Leonard & Merritt 1989), 




(4.9) 
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The combined radial and proper motion data yield /?o = 8.2±0.9 kpc (Genzel et al. 2000). 

The visual binary method (Salim & Gould 1999; Jaroszynski 1999) uses precision measurements 
of proper motions and radial velocities of stars orbiting the MBH to provide a direct geometric determi- 
nation of the Sun-MBH distance (§4.1.2). The only assumption is that the orbit is Keplerian. By virtue 
of its simplicity, this method is essentially free of systematic uncertainties due to uncertainties in the 
astrophysical modeling. 

The first application of the visual binary method was based on the orbit of the star S2 (Eisenhauer 
et al. 2003). Since the data was not restrictive enough to constrain the 13 free parameters (§4.1.2), it 
was assumed that the 3 velocity components of the MBH can be taken as zero, in keeping with the low 
peculiar velocity measured in the radio (~25kms _1 , Backer & Sramek 1999; Reid et al. 1999, 2003a) 
and the small error in the estimated velocity of the Solar system toward the GC. This analysis yielded 
/? = 7.94±0.42 kpc, e = 220.7±12.7kms- 1 and ra = (3.6±0.6)x 1O 6 M . A subsequent analysis with 
more data (still assuming zero velocity for the MBH) yielded compatible values, /?o = 7.62±0.32 kpc 
and ra=(3.6±0.3)x 1O 6 M (Eisenhauer et al. 2005). 

4.2 Constraints on non-BH dark mass alternatives 

A black hole, by its nature, is an object whose existence is hard to prove. Other than demonstrating the 
absence of a hard surface and a diverging gravitational redshift (the signature of an event horizon as seen 
by a stationary distant observer), one has to make do with placing ever tighter constraints on the volume 
enclosing the dark mass. Such limits can then be used to exclude non-BH dark mass models (that is, 
mass distributions that extend beyond their event horizon), such as a cluster of stellar-mass dark objects 
(§4.2.1) or "balls" of exotic particles (§4.2.2). Even if the existence of a MBH can be established, it 
is still possible, indeed quite likely, that it is surrounded by an extended distribution of dark mass, for 
example a dense cluster of massive compact stellar remnants that have sunk to the bottom of the potential 
well by the dynamical process of mass segregation (§4.2.3, §7.2). The stellar orbits tracked around the 
central dark mass probe the gravitational potential there at distances as small as ~ 10 light hours (~ 10 15 
cm). The information provided by these orbits offers the strongest case yet for the existence of a MBH 
in a galactic nucleus. Gravitational lensing may also be used to probe the distribution of the dark mass 
(§5.2.2, §5.2.3). Figure (4.5) summarizes the constraints on the size and density of the dark mass that 
are discussed below. 

4.2.1 Cluster of dark stellar- mass objects 

The most conservative dark mass alternative to a MBH is a very dense cluster of dark, or low-luminosity 
stellar mass objects. Possible candidates include very low mass stars (M* > 0.09M Q ), brown dwarfs 
(0.01 Mq <M* <0.09M q , failed stars that are not massive enough to ignite hydrogen), planets (M* < 
0.01 Mq, where gravity is negligible compared to the electrostatic forces) or compact stellar remnants: 
white dwarfs (WDs) (M* < 1.4M ), NSs (1.4M <M* < 3M ) and stellar mass black holes (M* > 
3M ). It is not at all obvious that stellar evolution and stellar dynamics can form a cluster of dark 
objects that is massive and concentrated enough to play the role of the central dark mass in the GC. 
However, irrespective of the formation history of such a hypothetical cluster, it is possible to invoke 
stability arguments to obtain robust upper limits on the cluster lifespan against the dynamical processes 
of evaporation and run-away collisions (Maoz 1998; §3.1.2). The lower this upper limit is compared to 
the age of the GC, the less likely is the dark cluster hypothesis. 

Maoz (1998) obtained such limits for a sample of galaxies, including the Milky Way, where a dark 
mass m was detected inside a radius r m . To obtain the most conservative limit, it was assumed that (1) 
the hypothetical cluster has the least centrally concentrated observed density profile (a Plummer model, 
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Figure 4.5: Compilation of the constraints on the nature of the dark mass in the GC (Schodel et al. 
2003). Filled circles: Lower limits on the size and density of the dark mass. Grey area: Constraints 
on the size of the variable X-ray emitting region (Baganoff et al. 2001). Filled squares: Different dark 
mass candidates, including the visible star cluster and its central cusp, the heavy fermion ball (Tsiklauri 
& Viollier 1998), the boson star (Torres et al. 2000) and a 3x 1O 6 M MBH. Horizontal lines: Lifetimes 
of hypothetical dark clusters of astrophysical objects, neutron stars, white dwarfs, or stellar black holes 
(Maoz 1998). The tightest constraint on the mass density, p> O(0.1 gem -3 ) (not shown) come from 
the radio proper motion of SgrA* compared to the velocities in the surrounding star cluster (Reid & 
Brunthaler 2004) under the assumption that the dark mass is enclosed in the millimeter radio source 
(Doeleman et al. 2001). All configurations are excluded except a black hole and a boson star. (Reprinted 
with permission from the Astrophysical Journal) 



p=po/ l+(r/ro) ), so as to minimize the collision rate; (2) it is composed of equal-mass objects, 
since otherwise mass segregation would accelerate its dynamical evolution; (3) the objects are at zero 
temperature, so as to minimize their radius for a given mass and thus minimize the collision rate (Eq. 
3.11). A maximal lifespan t m was obtained for each galaxy in the sample, based on the measured values 
of m and r m , and by choosing the dark object that maximizes t m (0.6A/q WDs in the case of the GC). 
The dark mass in the GC was estimated to have t m < 10 8 yr <Cf#. This argues against the possibility that 
the dark mass is a dark cluster. 
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4.2.2 Exotic elementary particles 

One interesting MBH alternative is that the dark masses in the centers of galaxies are the high density 
peaks of the cold dark matter distribution that is dynamically observed to exist on galactic and inter- 
galactic scales (Tsiklauri & Viollier 1998; Munyaneza & Viollier 2002). This idea's appeal lies in its 
economical use of a single component to explain dark mass on all scales. If dark matter is composed 
of fermionic elementary particles, then a massive compact "Fermion Ball" could support itself against 
gravity by degeneracy pressure, without having a hard surface (thereby avoiding tell-tale emission from 
matter impacting on the surface) and without being relativistic. Dynamically, such a configuration would 
be indistinguishable from a MBH for any process that occurs outside the ball. Deviations from Keple- 
rian motion will become apparent only for orbits that cross through the ball. Also, a Fermi ball, unlike a 
low-mass MBH, will not necessarily induce tidal disruption since its size can exceed the tidal disruption 
radius (§6.1). 

The maximal Fermi ball mass that can be supported by the degeneracy pressure of a fermion of 
mass nif is set by the Oppenheimer-Volkoff limit (Oppenheimer & Volkoff 1939), the GR analog of 
the Chandrasekhar mass for WDs, maxm~ g^' 2 m\/rrv\, where gf is the spin degeneracy factor of the 

fermion and mp = y/Hc/G is the Planck mass. By choosing a small enough m,f it is then possible to 
explain the entire range of galactic central dark masses as Fermi balls. The maximal measured central 
dark mass (maxm < 3 x 10 9 M Q ) then implies an upper limit on the fermion mass of 



maxnif = 0.62 1/4 x -^—tru ~ (g f /2) 1/4 15keVc" 2 . (4.10) 

v ' V maxm v ' 

The orbit of SO-2 is consistent, up to the measurement errors, with a Keplerian orbit around a point 
mass m. This one orbit can be used to rule out the universal Fermi ball alternative to MBHs (Munyaneza, 
Tsiklauri & Viollier 1998; Schodel et al. 2002). The lack of deviations from a Kepler orbit near the 
orbital periapse implies that the Fermi ball must be wholly contained within the periapse r p . The velocity 
of S2 at periapse is < 0.02c, and so a Fermi ball that fills the enclosed volume can be treated in the non- 
relativistic limit. The polytrope mass-radius relation (Chandrasekhar 1939) for a cold, non-relativistic 
degenerate Fermi ball results in a lower limit on the fermion mass 

minm / = 1.76( g/ /2)- 1/4 (^j'^mp ~ (^-^VOkeVc" 2 , (4.11) 

where £p = ^HG/c 3, is the Planck length. Thus, there is no fermion mass that can simultaneously 
accommodate the orbital constraints of S2 on the extent of the dark mass in the GC and support the 
largest measured central dark masses. This rules out the possibility that the same component singly 
explains both the diffuse dark matter on galactic and inter-galactic scales and the compact dark matter 
in the centers of galaxies. 

Another proposed exotic particle alternative to a MBH, which cannot be ruled out by the stellar 
orbits, is the "Boson Star". Such an object is composed of a scalar field and is supported against 
collapse by the Heisenberg uncertainty principle (Torres, Capozziello & Lambiase 2000). A wide range 
of boson star masses can be realized, including MBH masses, depending on the assumptions about the 
boson particle masses and their self-interactions. A boson star does not extend much beyond its event 
horizon, and is therefore strongly relativistic, but it does not possess a singularity, an horizon or a hard 
surface. Because both are so compact, stellar orbits cannot empirically distinguish between a MBH and 
a boson star. However, there are theoretical arguments against the dark mass being a boson star, since 
it is not clear that it can avoid accreting some matter and collapsing to a MBH. A boson star, unlike a 
MBH, has no capture orbits and so high velocity matter, such as that released in tidal disruption events, 
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will not accumulate in the center. However, this may not be enough to suppress the accretion of compact 
remnants on tightly bound orbits (§4.2.3) or gas flowing from the inner edge of an accretion disk. 

4.2.3 Extended dark mass around the MBH 

Even if a MBH does exist in the GC, it is still possible that there are other dark components around it, 
a cluster of massive compact objects (§3.3), a cusp of exotic particles, or a combination of the two. It 
should be noted that irrespective of the presence of the MBH, the Galactic dark matter halo is expected 
to peak at the GC (e.g. Navarro, Frenk & White 1996). The possibility that a steep, high density 
P~ Ph(r/rh)~ a cusp of self-interacting dark matter particles could form under the dynamical influence of 
the MBH (§3. 1.3) was investigated in the hope that the enhanced particle annihilation rate (oc p 2 ) could 
produce a detectable signal (Gondolo & Silk 1999; but see opposing conclusions by Merritt 2004b). 
The detection of strong TeV emission from the vicinity of Sgr A* (Kosack et al. 2004; Aharonian et al. 
2004), further adds to the interest in this possibility. However, with a (very uncertain) normalization 
of only p/, ~ 100 M Q pc~ 3 (Gnedin & Primack 2004), the mass in the dark matter particles is too small 
relative to the stellar cusp (Eq. 3.24) to be probed by the stellar orbits. 

The smoothed potential of an extended dark mass component can be probed by deviations from an 
exact Keplerian orbit that will be displayed by a star that passes through the mass distribution. Such 
orbits will display a retrograde shift of the angle of periapse (a receding rosette-like orbit), because the 
stellar acceleration at distance r from the center is proportional to the enclosed mass M(< r) (spherical 
symmetry assumed), and it decreases as the star approaches periapse. This shift is in the opposite sense 
of the effect expected for a GR orbit around a point mass (an advancing rosette-like orbit, §5.1). The 
retrograde periapse shift per orbit Alo due to a small spherical potential perturbation 5<j) on top of the 
potential of the MBH, 0=(/> MB h+<5<^, is given by (e.g. Landau & Lifshitz 1969) 

Au(a,e) = §j(jJ Q r 2 5<j)&u^ . (4.12) 

For a powerlaw extended mass distribution, parametrized as M e (< r) = m(r /r m j i l 2 ~ p (cf Eq. 3.17), the 
periapse shift per orbit can be evaluated numerically by (see also Ivanov, Polnarev & Saha 2005) 



3/2-p 2(1 _ g 2 )3 /2-p , w 4{2 -p)e+ [4(2-;7)-(3-2p)(l- g 2 )] cosa; 

e(l-2p) J (l+ecosw) 7 / 2 ^ UJ ' 

Mouawad et al. (2005) show that the lack of detectable deviations from Keplerian motion in the 
orbit of S2 places an upper limit of M e (< r a )/m< 0.05 on the extended mass, where r a ~0.01pc is 
the apoapse of S2. This upper limit still exceeds by a factor of ~ 30 the dynamical upper limit on the 
population of SBHs that can exist there (§3.3). Occasional orbital perturbations due to close encounters 
with individual compact objects can distinguish a dark matter component consisting of compact objects 
from a smooth particle dark matter distribution. Such perturbations may be detected by future very large 
telescopes (3> 10m) (Weinberg, Milosavljevic & Ghez 2005). 

The orbits are not sensitive to the mass distribution inside the smallest periapse (provided it is 
spherical). An upper limit on the fraction of the extended mass on these small scales (<fewxmpc) 
can be obtained from the Brownian motion of the mass that is associated with the radio source Sgr A* 
due to gravitational interactions with stars around it. The amplitude of the Brownian motion is strongly 
constrained by the radio limits on the peculiar motion of Sgr A* (motion in excess of the reflex motion 
induced by the Sun's orbit around the Galaxy). Brownian motion in a gravitating system is different 
from that in a gas, where the inter-molecular forces are short range forces. Equipartition arguments can 
not be directly applied. Inside r/, the potential is dominated by the MBH and the stars are "carried away" 



Au>(a,e) = 
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with it as it moves; conservation of momentum rather than equipartition dictates the relation between the 
MBH velocity and mass and those of the stars. Neither is it possible to assign a well-defined temperature 
to the stellar system; the MBH itself affects the velocity field ("temperature") of the stars near it (Merritt 
2004a), and even outside r/, the stellar distribution is generally not isothermal (Chatterjee, Hernquist & 
Loeb 2002b, 2002a). Reid & Brunthaler (2004) use the radio observations to first show that if the mass 
associated with Sgr A* is some fraction of the total central mass, then it must lie very close to the center 
of the region delimited by the stellar orbits, otherwise it would have orbited rapidly around the dark 
mass, in contradiction with the upper bound on the peculiar motion of Sgr A* (~ lkms -1 perpendicular 
to the galactic disk). They then proceed to assume that Sgr A* does lie in the center, and use N-body 
simulations of the Brownian motion to translate the upper bound on the peculiar motion to a lower 
bound of >0.4x 1O 6 M on the mass associated with Sgr A* (~0.1 of the dynamical mass measured 
by the orbits). This limit is consistent with the results of a more detailed formal analysis of Brownian 
motion in gravitating systems (Laun & Merritt 2004). It is unlikely that the large mass fraction (~0.9) 
that remains unconstrained is in fact non-BH dark matter. However, an order of magnitude improvement 
in the velocity resolution will be needed to confirm that the MBH indeed dominates the dark mass. 



4.3 Limits on MBH binarity 

Galaxies as observed today are thought to have formed by a chain of successive mergers of smaller 
galaxies. High redshift quasars and the large number of AGN found at smaller redshifts suggest that 
these mergers must also involve the formation of a binary MBH in the center of the newly merged 
galaxy, after the two MBHs rapidly sink in due to dynamical friction (§3.3). The binary subsequently 
decays, initially by dynamical processes, and finally by the emission of GW, and the two MBHs merge 
(see review by Merritt & Milosavljevic 2004). In a "major merger" the two MBHs have comparable 
masses, whereas in a minor merger the mass ratio is large. The formation history of the Galactic MBH 
is not known, but there are several arguments that suggest that the Galactic MBH did not undergo a 
major merger, at least since it reached its present mass scale. Its small mass can be easily supplied, 
over time, by the stellar winds near it (§7.3) and so there is no compelling reason to assume mergers. 
Another argument against a major merger is the observed stellar distribution around the MBH (§3.2). 
Initially the binary's orbit decays dynamically, by ejecting stars that pass near it out of the galactic core 
by 3-body interactions. A total stellar mass of the same order as the binary mass is ejected during a 
major merger, leaving behind a galaxy with a shallow central density distribution (Milosavljevic et al. 
2002). The GC, however, has a steep inner cusp. Furthermore, a binary MBH is expected to induce a 
tangential anisotropy in the stellar orbits by ejecting or destroying stars on radial orbits that approach it 
(Gebhardt 2003). This is not observed in the GC (Schodel et al. 2003). 

Minor mergers, on the other hand, could well have occurred. The presence of a ~ 1O 3 M IBH 
is in fact suggested by the apparently self -bound stellar cluster IRS 13, ~0.15 pc (projected) from the 
MBH (Maillard et al. 2004; but see opposing view by Schodel et al. 2005). It is also suggested by the 
"sinking cluster" model, which seeks to explain the young population in the center as originating in a 
dense, star forming stellar cluster with an IBH, which sank to the center by dynamical friction (Hansen 
& Milosavljevic 2003; §7.2.3). 

The binary MBH hypothesis can be constrained directly by the stellar orbits. The good fit of S2's 
motion to a Keplerian orbit suggests that a hypothetical binary MBH must be contained in a region 
much smaller than the semi-major axis of S2, ~5 mpc (see Jaroszynski 2000 for a quantitative analysis). 
However, it is very unlikely that such a tight binary (of moderate mass ratio) exists now in the GC, given 
its very short lifespan against orbital decay and coalescence by the emission of GW radiation. The orbit 
of circular binary MBH (masses ni\ and mj and semi-major axis an) would shrink to a point by the 
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emission of GW radiation in time 



5 



(4.14) 



256 G 3 (mi +m2)m\ni2 ' 



(Peters 1964; Shapiro & Teukolsky 1983 Eq. 16.4.10). For example, if m\ = m 2 = m/2= 1.5 x 1O 6 M 
and ai2=0.1 mpc, then the lifespan of the binary is only ?gw= 10 4 yr. MBH Binaries with smaller mass 
ratios are, however, much longer lived, and cannot be ruled out by this argument. 

The MBH motion in the plane of the sky, v ~2masyr _1 =76(/?o/8kpc)kms -1 , is obtained from the 
full orbital solutions (Ghez et al. 2005). A second BH orbiting the MBH would induce a reflex motion 
of the MBH around their common center of mass. This can be translated to an upper bound on the mass 
and distance of such a BH (assuming it lies outside the solved stellar orbits) of 



assuming a circular orbit in the plane of the sky and m% <Cmi = 3 x 10 6 M Q (see Hansen & Milosavljevic 
2003 for a a more detailed analysis of the astrometric constraints on m 2 and a) 

Tighter upper bounds can be obtained from the limits on the peculiar motion of Sgr A* in the radio 
of 18±7kms _1 in the plane of the galaxy and 0.4±0.9kms _1 perpendicular to it (Reid & Brunthaler 
2004). However, these are not motions measured directly relative to the central cluster, but rather 
measured relative to extragalactic radio sources, and so include also the reflex motion, @o/Rq, which 
can not be subtracted exactly, and the possible Brownian motion of the central cluster relative to the 
Galactic center of mass, which is unknown. The upper limit on the perpendicular motion of ~ lkms -1 
translates through Eq. (4.15) to an upper limit of m 2 < 10 4 ^/r/0.005 pcM (for r> 0.005 pc). This 
upper bound is still consistent with a 1O 3 -1O 4 M IBH in IRS 13. 

4.4 High-velocity runaway stars 

A possible stellar signature of the existence of a MBH, or possibly a binary MBH, would be the discov- 
ery of very high velocity stars tens of kpc away from the GC moving with v~ 0(1000 kms -1 ) on nearly 
radial orbits relative to the GC (the local escape velocity from the Galaxy at the Sun is ~ 600 km s -1 , 
e.g. Wilkinson & Evans 1999). The stars can acquire such high velocities by exchange interactions deep 
in the potential of the MBH, either by the break-up of a binary on a nearly radial orbit (Hills 1988), by 
an encounter with another single star, or an encounter with a binary MBH (Yu & Tremaine 2003). A 
slow-moving star with velocity Vqo far from the MBH will acquire near periapse a very high velocity, 
v 2 ~ 2Gm jrp^v 1 ^. A small velocity perturbation by a gravitational interaction with a second mass, Sv, 
can add v5v3>v 2 x3 /2 to the orbital energy of the star, and take away the same amount from the other 
mass (a similar mechanism forcefully expels tidal ejecta, see §6.1.2). Yu & Tremaine (2003) find that 
stellar binary exchanges or encounters with a binary MBH can populate the Galaxy with high velocity 
stars at a substantial rate (> 10~ 5 yr _1 ), resulting in up to ~ 10 3 high velocity stars (v^ > lOOOkms ) 
inside the Solar circle. 

Intriguingly, Brown et al. (2005) report the discovery of the highest velocity star ever observed in 
the Galactic halo, moving with a radial velocity of 710kms _1 at a distance of ~55kpc. Its direction of 
motion could be consistent with radial motion away from the GC, if the (unknown) transverse velocity 
components are small. The stellar spectrum does not identify its spectral type unambiguously, but one 
of the two possibilities is a B9 MS star, such as those found in the S-cluster (the other possibility is a 
blue horizontal branch giant). 




(4.15) 
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5 Probing post-Newtonian gravity near the MBH 



General Relativity (GR) is the least tested of the theories of the four fundamental forces of nature. An 
important goal of the study of stars near the Galactic MBH is to detect post-Newtonian effects and 
probe GR in the weak and strong field limits near a super-massive object. Three classes of effects are 
discussed in this context. The first is relativistic celestial mechanics. Stars are observed moving with 
velocities of up to v~ fewxO.Olc at distances of < 10 3 rj from a mass of ~4x 10 6 M Q . This regime of 
GR dynamics is virtually unexplored 14 (§5.1). For comparison, the high precision confirmation of the 
predictions of GR in the Hulse-Taylor binary pulsar (PSR 1913+16) were for masses of ~ 1.4M Q with 
v ~ 0.003c at r p ~ 2 x 10 5 r 5 (Taylor & Weisberg 1989). The second class of weak field phenomena that 
may be relevant in the GC is gravitational lensing (§5.2). The lensed stars are not necessarily related to 
the stars near the MBH, since it is only their projected angular distance to the MBH that has to be small, 
but their true (3D) position can be anywhere behind the MBH. The third class of GR effects that may 
be relevant in the GC is the emission of GW from very low-mass stars that spiral into the MBH. This 
dissipative process is discussed separately in (§6.2.2). 



5.1 Relativistic orbital effects 

At present, all the available orbital data can be adequately modeled, within the measurement errors, in 
terms of Newtonian motion. With improved resolution, higher precision and a longer baseline, devi- 
ations from Keplerian orbits may be detectable. This will likely be a challenging task. The chances 
of success will be much improved if future observations detect fainter stars on even tighter and more 
eccentric orbits than the relatively bright ones observed today. This is because post-Newtonian effects 
are larger the tighter the orbit, while the uncertainty due to the dynamical influence of the unknown 
amount of extended mass around the MBH (§4.2.3) is minimized. However, it is not clear how many 
faint, low-mass stars exist in the inner few mpc, since mass segregation and possibly also collisional 
destruction by the massive remnants there are expected to suppress their numbers there. 

The highest stellar velocity recorded to date in the GC was that of the star SO- 16 during periapse 
passage (Ghezetal. 2005), which reached ~ 12000 km s" 1 (/? ~ 0.04, f3=v/c ) at a distance of r p ^600r s 
from the MBH. This corresponds to a relativistic parameter T = r$/r at periapse of Y(r / ,)~' 1.6 x 10~ 3 . 
For comparison, the relativistic parameter on the surface of a WD is a few x 10~ 4 and on a NS a few x 
10 _1 . The relativistic parameter changes along the orbit. When the orbit is only mildly relativistic, it 
obeys a nearly Keplerian energy equation (Eqs. 3.15, 3.16) and T can be expressed as 

T = /? 2 + ^/? 2 + r . (5.1) 
c L a 

When e S> and the star is near periapse, f3 2 3> To and T ~ 1 . The deviations from Newtonian mechan- 
ics become more pronounced as T increases. For T <C 1 (as is the case here) it is useful to expand the 
post-Newtonian deviations in /3, and classify each effect by the order of its /3-dependence (equivalently, 
by the order of its \frs/r p or \Jr^ja dependence). 



Gas processes probe the dynamics much closer to the MBH, but are also much harder to interpret both 
observationally and theoretically. Emission from accreting gas is observed from the strong field near the event 
horizon of MBHs. There is tentative evidence of gravitational redshift in the Fe Ka line profiles in AGN (e.g. 
Fabian et al. 1995, but see e.g. Page et al. 2004). The emission from Galactic SBHs and AGN sometimes display 
quasi-periodic oscillations that can be interpreted as coming from the last stable orbit around a spinning BH (see 
Genzel et al. 2003a; Aschenbach et al. 2004 for such evidence from the Galactic MBH). 
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5.1.1 Gravitational redshift 

The measurement of stellar proper motions and accelerations in the inner arcsecond of the GC (Eckart 
& Genzel 1996; Ghez et al. 1998) preceded by more than a decade the measurement of radial velocities 
for these stars (Ghez et al. 2003b; Eisenhauer et al. 2003). For this reason, most theoretical studies of 
the detection of post-Newtonian orbits in the GC focused exclusively or mainly on proper motions and 
effects related to them, such as periapse shift (Jaroszynski 1998b; Fragile & Mathews 2000; Rubilar & 
Eckart 2001; see also general treatise on relativistic celestial mechanics by Brumberg 1991). However, 
with the advent of adaptive optics assisted IR imaging spectroscopy, it now appears that it is the radial 
velocities, rather than the proper motions, that may provide the tightest constraints on deviations from 
Newtonian orbits. The current quality of spectroscopic observations allows the determination of radial 
velocities from the Doppler shift of stellar spectra to within <5v~25kms~ ! , or <5A/A~ 10~ 4 (Eisenhauer 
et al. 2005). This limit is partly instrumental and partly due to the fact that the stars very close to 
the MBH have hot ("early type") stellar spectra (§2.2, §7) with only a few strong IR lines, which are 
typically broad with (v rot sin i) ~ lSOkms" 1 (Ghez et al. 2003b; Eisenhauer et al. 2005; Fig. 4.1). 

The observed radial velocity [3 r can be expanded in terms of the true magnitude of the 3D stellar 
velocity f3 as 



Pr = AA/A = B +B l (5+B 2 p 2 + O(P 3 ). (5.2) 

Several Newtonian and post-Newtonian terms contribute to the three leading orders. The relativistic 
Doppler shift 

can be decomposed to the first order classical Doppler effect (i? is the angle between the velocity vector 
(3 and the line of sight) and to the second order transverse relativistic Doppler redshift. The gravitational 
redshift of a source at distance r from the MBH as measured by an observer at infinity (Eq. 5.1) is 

& = ^= (l-y)" 1/2 -l * \T(r) = \ [T (a) + (3 2 ] . (5.4) 

The gravitational redshift contributes a constant term, which depends on the orbital semi-major axis, 
or equivalently, the energy, and a second order term of the same magnitude as the transverse Doppler 
redshift. 

The Roemer (or light travel time) delay is a Newtonian effect caused by the fact that the light travel 
time from the star to the observer changes with the orbital phase when the orbital plane is not observed 
face-on (see also Loeb 2003). The phase-dependent mismatch between the signal's emission and arrival 
times leads to apparent deviations from the Newtonian orbit. The exact value of the effect depends on the 
orbital parameters. Its magnitude and /3-dependence can be estimated by considering a circular edge-on 
orbit, in which case the maximal delay is the light crossing time of the diameter of the orbit 5t~2a/c and 
the maximal apparent velocity shift is (3r ~ \(3 r \5t ^ (Gm / ca )8t = 20 1 . The Roemer time delay enters 
radial velocity measurements as a 0(P 2 ) correction, with a pre-factor Br that is typically larger than that 
of either the transverse Doppler shift or the gravitational redshift. Note that the Roemer delay expresses 
itself in the proper motions as a relative shift 5r/a of order 0((5), since Sr / a ^ (2a / a)(5t /P)~2/?. 

The observed radial velocity (Eq. 5.2) can therefore be written as 



/?© + & !g al + A,*+2 T 



+ [costf]/?+[l+fi R ]/3 2 + 0O3 3 ), (5.5) 
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where the constant coefficient Bo is expressed in terms of (1) the constant velocity shift (3 & due to 
the compound motion of the Sun and the Earth relative to the GC as well as the blueshift due to the 
potential wells of the Sun, Earth and planets ; (2) the constant velocity shift /3 Z)ga i due to redshift by 
the potential of the Galaxy; (3) the gravitational redshift /3 Z)i , due to the star's potential 15 , and (4) The 
constant contribution To/2 from the gravitational redshift of the star in the potential well of the MBH. 

Since stars are observed with j3 ~ few x 0.01 at periapse and since #2 = 1 +Br~0(1), it follows that 
second order effects are detectable with existing instruments, B2P 2 ~ 10~ 3 > 5X/X ~ 10" 4 . Third order 
GR effects, such as frame dragging (§5.1.3), Gravitational lensing and the Shapiro delay (relevant only 
for a nearly edge-on orbit and even then very small, Kopeikin & Ozernoy 1999), or GR periapse shift 
(§5.1.2), which is only a A/3~ Au;r / ,/r p ~ C(/3 3 ) (Eq. 5.6) effect in the radial velocities (although it is 
a 0(P 2 ) effect in the proper motions) will not be detectable with the available radial velocity data. It 
should be noted that second order effects involving proper motion, such as the relativistic periapse shift, 
will be very hard to detect with the present astrometric precision. 

In practice, the ability to measure the post-Newtonian deviations from a Newtonian orbit requires an 
accurate determination of the Keplerian parameters, and so reliable proper motion data are still needed. 
In particular, radial velocity data for a single star on a 0{[3 2 ) post-Newtonian orbit can always be ex- 
pressed as a Newtonian orbit (including the classical Roemer effect) with suitably modified orbital 
parameters (the degeneracy between To and /?© in Eq. 5.5 is an expression of this fact). Multiple 
orbits and proper motions are needed to break this degeneracy (Zucker & Alexander 2005, in prepara- 
tion). Post-Newtonian orbit solving simulations with proper motions and radial velocities, such as those 
conducted by Weinberg, Milosavljevic & Ghez (2005) for Newtonian orbits, are needed to assess the 
practical observational requirements for detecting these effects. 

A similar approach to detecting post-Newtonian effects with radio pulsars, which may exist around 
the MBH, was studied by Pfahl & Loeb (2004). The problems are that it is unknown whether such a 
population exists, and that radio observations at the required sensitivity will be extremely challenging 
(cf §3.3). 



5.1.2 Periapse shift 

The deviation of the GR potential from the Newtonian l/r point mass potential leads to rosette-like 
orbits with a prograde shift in the argument of pericenter lo (see footnote 12) of Am per orbit, which for 
a non-rotating black hole is (Weinberg 1972, Eq. 8.6.11) 

A"=7rV*- (5-6) 

The periapse advance is a 0(f3 2 ) effect. The effect of the periapse advance on the orbital shape is 
maximal at the apoapse. The change in the angular position of the apoapse, p a , per orbit (neglecting 
corrections for the orbital orientation) depends only on the orbital eccentricity, 

a(l +e) 3irrs 

A Pa ~ Aw- = — . (5.7) 

la R (l-e)R 



15 The constant terms /3 Q , /3 : . ga i and /3 zir are very small compared to B2P 2 (cB^li 1 ~ O(100 kms" 1 ) at periapse 
for S2/S0-2). The peculiar radial motion of the sun relative to the GC is estimated to be < lOkms" 1 (Gould & 
Popowski 1998) and the transverse Doppler redshift is 9o 2 /2c — 0.08 kms" 1 . The constant gravitational redshift 
due the the Galactic potential (in addition to that of the BH, which is taken into account by Eq. 5.4) is c/3 Zjga i~ 
Oo 2 /c~0.16kms"'. The gravitational blueshift due to the proximity of the Sun is only -v 2 B /c = -0.003kms- 1 , 
where v® ~ 30kms _1 is the Earth's orbital velocity. Likewise, the redshift due to the star's potential and the 
blueshift due to the Earth's potential are only e/3 z * ~ 1 kms" 1 and -GM§ /ci?® ~— 0. 0002 kms" 1 respectively. 
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For S2/S0-2, the shortest period star detected to date (P~ 15yr, a~4.5mpc, e~0.88 Eisenhauer et al. 
2005; Ghez et al. 2005), the shift is Ap a ~ 1 mas per period. A similar shift is predicted for the star 
S14/S0-16 with the most eccentric orbit solved to date (P~38yr, e~0.94). A more efficient use of 
the orbital data is to fit an orbital model including the periapse shift (i.e. with a variable argument of 
periastron to) to all the data points. Either way, unless a very eccentric, very short period star is discov- 
ered, the detection of the periapse advance will require very long monitoring (see detailed discussions of 
the observability of the effect by Jaroszynski 1998b; Fragile & Mathews 2000; Rubilar & Eckart 2001; 
Weinberg, Milosavljevic & Ghez 2005). 

The prospects of discovering the GR periapse advance are further complicated by the likely ex- 
istence of a distribution of faint stars and dark stellar remnants around the MBH, which introduces a 
retrogade periapse shift (§4.2.3, §3.3). Measurements of the periapse shifts in several stellar orbits will 
be needed to disentangle the two countervailing effects (Rubilar & Eckart 2001). 



5.1.3 Frame dragging 

There is some evidence from the variability of the IR and X-ray accretion flares that the Galactic MBH 
is spinning in at least half its maximal allowed rate (Genzel et al. 2003a; Aschenbach et al. 2004). A 
spinning mass with specific spin S. induces secular changes in the trajectory of a test particle orbiting 
it through gravitomagnetic GR effects known as "frame dragging" or Lense-Thirring precession 16 . The 
S, -induced change in the longitude of the ascending node Q (see footnote 12) corresponds to the pre- 
cession of the orbital angular momentum vector J around the MBH spin vector S,. The change in the 
argument of pericenter, to, corresponds to a periapse shift, in addition to the S. -independent shift (Eq. 
5.6). 

In the weak field limit, the change per orbital period in Q is (Thirring & Lense 1918) 

and the period of the precession is 

P Qlt = 4it^- — - , (5-9) 

c s \r s J 

where s = S,/(Gm/c) is the dimensionless spin parameter (0<s< 1). The change per orbital period in 
the argument of pericenter is 

Aw L T = -3ASl LT cos/, (5.10) 

where cos / = J • S,//5». The frame dragging shifts are 0(/3 3 ) effects. 

The predicted magnitude of the effect in the GC is small (Jaroszynski 1998b), and it is unlikely that 
it will be measured directly from the orbits even with the proposed Thirty Meter Telescope (Weinberg, 
Milosavljevic & Ghez 2005). However, over the lifespan of one of the young S-stars near the MBH, 
4 ~ 10 Myr, the accumulated precession can be substantial if the orbit is eccentric enough. Levin & 
Beloborodov (2003) note that for the eccentric, short-period star S2/S0-2 (P*= 15.2 yr, e=0.87, Schodel 
et al. 2002) the Lense-Thirring period is Pq lt ~ O(10 7 yr). They propose that S2 was initially a member 
of the thin star disk (§2.2), orbiting in the plane of the disk but on a particularly eccentric orbit, and 
that it was gradually dragged out of the plane by the Lense-Thirring effect, to the point where its orbit 
today appears unrelated to the disk. Such effect, if measured in two stars could be used to determine the 



16 A related but distinct gravitomagnetic effect, not discussed here, is the Lense-Thirring precession of the spin 
of the test particle itself. 
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spin axis. They further propose that this could be the origin of the other S-stars as well. However, the 
orbital solutions of several other S-stars do not support this idea, since for a large fraction of the S-stars 
Pn LJ ^>U (Eisenhauer et al. 2005). 

5.2 Gravitational lensing 

Gravitational lensing of stars by the Galactic MBH differs from the other stellar processes discussed 
here in that it is a process that affects the light emitted by the stars, rather than the stars themselves, 
and in that the lensed stars are typically well outside the dynamical region of influence of the MBH. 
However, the lensed images can masquerade as stars near the MBH or as accretion flares from the MBH 
and complicate the interpretation of the observations. On the other hand, gravitational lensing may be 
used to probe the dark mass and the stars around it, and to locate the MBH on the IR grid, where the 
stars are observed. 

There is yet no clear evidence for gravitational lensing by the MBH. This is consistent with various 
estimates that suggest that there are not enough luminous sources behind the MBH for gravitational 
lensing to be important for present-day observations, although future, deep observations may detect 
lensing events (Wardle & Yusef-Zadeh 1992; Jaroszynski 1998a; Alexander & Sternberg 1999; Alexan- 
der & Loeb 2001). Nevertheless, it is worthwhile to consider the possible roles of gravitational lensing 
in the observations and study of the Galactic Center. This is important not only in anticipation of future 
observations, but also because the estimates of the lensing probability are quite uncertain (they involve 
models of the unobserved far side of the Galaxy) and may under-estimate the true rate of such events. 

5.2.1 Point mass gravitational lens 

The MBH dominates the central potential and so gravitational lensing in the GC can be described to first 
approximation as lensing by a point mass. A light ray arriving from a background source to the observer 
is bent by an angle a L as it passes near a point mass lens (Fig. 5.1). In the small angle approximation, 
which is valid as long the impact parameter bi is much larger than the event horizon Qjl 3> r$), the 
bending angle is given by 



(see Schneider, Ehlers & Falco 1992 for a comprehensive treatment of gravitational lensing). The 
bending angle of a point mass gravitational lens diverges towards the center (&l=0) and decreases as the 
impact parameter grows. This is in the opposite sense to a spherical glass lens, where the bending angle 
is zero when the ray goes through the lens center and increases with the impact parameter. Thus, unlike 
a spherical lens, a point mass gravitational lens does not produce a faithful image of the lensed source, 
but rather breaks, warps and/or flips the image. A point mass lens creates two images of the source, one 
on either side of the lens. The two images are always in focus at the observer, regardless of the distance 
of the source behind the lens. By symmetry, the two images, the lens and the (unobserved) source all lie 
on one line (Fig. 5.2). The typical angular cross-section of the lens is given by the Einstein angle, 



where Dol is the observer-lens distance (Rq for the Galactic MBH), Dl$ is the lens-source distance, 
and Dos is the observer-source distance 17 . The physical size of the Einstein angle at the source plane is 

17 In a flat background spacetime, which is the relevant background for lensing on Galactic scales, Dos=Doi + 
Dis. In curved spacetime, which is the relevant background for lensing on cosmological scales, the distances are 



a L = 2r s /b L 



(5.11) 




(5.12) 



T. Alexander / Physics Reports 000 (2005 ) 000-000 5 1 




Figure 5.1: Schematics of gravitational lensing by a point mass in the small bending angle limit. An 
observer (0) sees light rays from a background source (S), which is located at a distance Dos away, 
deflected by angles an, olli in the potential of a massive lens (L), which is located at a distance Dot 
from the observer and Dls from the source. The lens mass and the distances define an Einstein ring of 
angular opening Be (Eq. 5.12), perpendicular to the optical axis (O-L). The source (S), which is not 
observed directly, is located at an angle yO^ relative to the optical axis. It appears as two lensed images 
(Ii, I2) at angular positions x\9e and *2#e» one inside the Einstein ring and one outside it. 

R e =Dos@e- At a fixed observer-lens distance, 9e reaches a maximal value for a source at infinity, 

0£ = 0oo\/P-, #oc=l/4^~ 2 " intheGC - ( 5 - 13 ) 

V D os V c 2 D 0L 

The relation between y, the angular position of the source relative to the observer-lens axis (the 
optical axis) and x\£, the angular positions of the images, can be derived from the geometry of the light 
paths in the small bending angle limit (Fig. 5.1), 

y = Xi,2-lAi )2 , or *i, 2 = ^ (y±V4+y 2 ) (5-14) 

where x l2 and y are measured in terms of 6 E and x 2 <0 by definition. The angular separation between 
the two images is O(0e)- Gravitational lensing conserves surface brightness, and so the magnification 
A 12 in the flux of each image relative to that of the unlensed source is proportional to the apparent 
change in the angular area of the source, 

Ai, 2 = |ay/ax li2 |- 1 = |l-x7 4 2 |" 1 . (5.15) 
the angular diameter distances, and this simple sum no longer holds. 
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The primary image at x\ is always magnified. The secondary image at x 2 can be demagnified to zero. 
The two magnifications obey the relations 



A x =A 2 + l>l, A=A { +A 2 = -. (5.16) 

yVy 2 + 4 

When y = the magnification formally diverges and the image appears as a ring of angular size 9e, the 
Einstein ring. This divergence is avoided in practice by the finite size of the source (e.g. a star). Finite 
sized sources are also sheared tangentially around the Einstein ring as the magnification increases. In the 
limit of high magnification, or small source angle (y<C 1), A~ 1 jy. It is customary to define A(y=l)= 1.34 
as the threshold of a microlensing event. 

Relative motion between the source, lens and observer will manifest itself as two moving images 
with time-variable magnifications (Fig. 5.2). The term "macrolensing" is used to describe the situation 
where the images can be resolved, and "microlensing" the situation where the images are not separately 
resolved, and only their joint variable magnified flux is observable. When the relative angular velocity 
y between the source and lens is constant, the microlensing light-curve is given by substituting y(t) = 
b ; o + 3 ;2 ( f-f o) 2 ] 1//2 hiEq. (5.16), where yo is the angular impact parameter of the source trajectory relative 
to the lens and to is the time of maximal magnification when y=yo (Fig. 5.2). Microlensing light-curves 
are distinguished from other astrophysical flares by their characteristic time-symmetric functional form 
and by their wavelength-independence, which is a consequence of the equivalence principle. 



5.2.2 Gravitational lensing by the Galactic MBH 

In order to plan the observational strategy for detecting gravitational lensing, or to estimate how likely 
it is that an observed flare is due to lensing, it is necessary to evaluate the detection probability. Two 
quantities are commonly used for this puipose, the optical depth and the lensing rate. The probability 
Pl = 1 -exp(-r^) for finding at any given time a star lensed by the Galactic MBH is defined in terms 
of an optical depth tl = n+iTR E dDos, which expresses the mean number of background sources 
within an angle 9e of the line of sight to the MBH. Rough estimates predict 7x^0(1) (Alexander 2001; 
Alexander & Loeb 2001) for lensing by the MBH (assuming no limits on the photometric sensitivity). 

The optical depth is a useful quantity in situations such as Galactic microlensing searches (tl ~ 
O(10~ 6 ), e.g. Afonso et al. 2003), when there are many possible lensed lines of sights, and ^r^ 1 ~ 10 6 
background stars must be monitored simultaneously to find the rare one that is lensed by an intervening 
star. The observational situation for gravitational lensing by the MBH in the Galactic Center is different 
because the position of the lens is known and fixed, and so there is only one line of sight to monitor. 
The optical depth does not take into account the relative motions of the observer, lens and source, which 
reshuffle their random alignment and introduces a timescale to the problem. A more relevant quantity 
in this situation is the lensing rate Tl, which indicates how long it is necessary to monitor the GC on 
average to detect a lensing event with peak flux above a flux detection threshold Fq, 



T L (> Fo) ~ / w -^dD 0S , for A > ; " , (5. 17) 

4nDi 



2Re 

I Dot A ^l^ u bs 

where v± is the velocity of the sources perpendicular to the line of sight relative to the optical axis, 
including the contribution from the reflex motion due to the observer, L+ is the source luminosity and 
A> 1 is assumed. For practical applications, equation (5. 17) has to be modified to take into account the 
range of stellar luminosities and velocities, dust extinction, the total duration of the observations T and 
the sampling rate AT (the mean duration of events amplified by more than A is 1l = ttRe /2Av± ; only 
events with AT <ti<T can be detected) (Alexander & Sternberg 1999). 
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Figure 5.2: A time sequence (left to right, top to bottom) of a gravitational lensing event of a point 
source by a point mass. The two resolved images and their trajectories (fuzzy circles, whose size and 
shade are proportional to the flux, with curved lines tracking the image trajectories) are superimposed 
on a plot of the microlensing light-curve (total flux F from both images as function of time t in arbitrary 
units). Note that the tangential shear of the images is not shown here. The background source, which is 
not observed directly, (open circle with straight line tracking the source trajectory) moves in projection 
from left to right behind the MBH (open circle at center) with an impact parameter of 0.1#£. The two 
images move in tandem clockwise about the Einstein ring (large dotted circle). The strongly magnified 
image (top) is always outside the Einstein ring and is always brighter than the source. The weakly 
magnified image (bottom) is always inside the Einstein ring and can be strongly demagnified. 

The observational limitations, Fq, T and AT and the telescope's angular resolution 5(f), place restric- 
tions on Dls and the impact parameter for which sources can be detected, and affect the appearance of 
the lensing event, the typical timescales and the peak magnifications that are likely to be observed. The 
angular distance between the two images close to peak magnification is ~2#£ (Eq. 5.14). A lensed star 
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will appear as a microlensing event if it is closer than D M behind the MBH, 




<5(/>/10mas 



) 



2 



{0oo/W) 2 -l 



0.2pc 



#oo/2" 



(5.18) 



Stars farther away will appear as two resolved images. The Einstein angle grows from zero to ^with 
source distance behind the MBH (Eq. 5.13), and so the lensing cross-section and the lensing duration 
(for a given max A and for constant vj_) also increase. In addition, high magnification events typically 
have longer time scales because the source trajectory must have a smaller impact parameter and so 
spends more time inside the Einstein radius. Therefore, observations with limited temporal sampling 
will tend to pick out high magnification events. 

A couple of early detections of transient flaring events very close to Sgr A* (Genzel et al. 1997; Ghez 
et al. 1998) were considered as lensing event candidates. For one of these a light curve was recorded, 
but as it was under-sampled only estimates of a timescale (~ 1 yr) and a typical magnification (~A > 5) 
could be derived from it. The a-posteriori probability of detecting a lensing event (for the relevant obser- 
vational limitations) was estimated to be only 0.5%, but on the other hand, the observed timescale and 
magnification were close to the expected median value (Alexander & Sternberg 1999). The interpreta- 
tion of this event remains inconclusive. The likely alternative to lensing is that these were faint stars 
near the MBH that were observable only intermittently when they were well separated from the more 
luminous stars nearby. It should be noted that the T^-band IR flares observed from Sgr A* have a typical 
timescale of 0(1 hr) (although year-long variability is observed in longer IR wavelengths, Clenet et al. 
2004a). They are too frequent to be microlensing of near background stars and furthermore, their irreg- 
ular light-curves and quasi-periodic oscillations are inconsistent with the time-symmetric microlensing 
light-curves (Genzel et al. 2003a). 

Macrolensed image pairs, if hidden among the many stars observed around the MBH, can be used 
to pinpoint the MBH on the IR grid directly from the astrometric data, independently of other methods 
(§4.1) and with no assumptions about m, Ro or the distribution of stars behind the MBH (Alexander 
2001). The MBH lies on the line connecting the two images of any background source it gravitationally 
lenses (Fig. 5.2), and so the intersection of these lines fixes its position. The point-mass lensing equa- 
tions (Eqs. 5.14, 5.16) imply simple relations between the fluxes of the two images, and between their 
angular positions and velocities and as measured relative to the projected position of the MBH, 



where 9 r j is the radial component of the projected angular velocity of image i (/= 1,2), 6 t j is the tan- 
gential component and F, the flux of the image. The simultaneous search for the position of the MBH 
and for lensed image pairs proceeds by enumerating over a grid of possible MBH positions. For each 
trial position, scores are assigned to all possible pairs of stars according to how well they satisfy Eq. 
(5.19). The MBH position is then chosen as the one maximizing the total lensing score, and the likeliest 
lensed image pair candidates can be recorded for spectroscopic validation (lensed image pairs should 
have exactly the same spectra, up to small deviations due to differences in the dust extinction along the 
two light paths). Such a search on an early astrometric data set produced a statistical detection of the 
MBH inside the error region on the center of acceleration (Fig. 4.3), with a chance detection probabil- 
ity of only 5 x 10~ 4 , where most of the signal came from the intersection of two pairs of lensed image 
candidates, tentatively identified as far background blue and red supergiants (simple models of the dis- 
tribution of stars and dust in the far side of the Galaxy predict none within 6e of the MBH). However, 
a subsequent recalibration of the data, and additional color information on the lensed image candidates 




(5.19) 
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(Clenet et al. 2004b), cast strong doubt on this result. This method has yet to be applied to the newer, 
more comprehensive and accurate astrometric data sets. 

The MBH can modify the observed stellar surface density inside ~ 9^, (Wardle & Yusef-Zadeh 
1992). A lens magnifies by enlarging the angular size of the unlensed sky behind it. Since surface 
brightness is conserved, the fluxes of sources are magnified by the same amount. If all the stars can 
be detected even without magnification, then the effect of lensing is to decrease the surface density of 
sources. However, if the magnification reveals faint stars, which could not be observed unlensed, then 
the observed surface density can be larger than the unlensed density if enough faint sources are mag- 
nified above the detection threshold to over-compensate for the decrease in surface density ("positive 
magnification bias"), or it could be smaller if the converse is true ("negative magnification bias"). The 
lensed luminosity function (stellar number surface density S per flux interval) is related to the unlensed 
one by 





= A" 2 
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In many cases the luminosity function is well approximated by a power-law, dS/dF oc F~ s (§2.1). For 
s=2 the decrease in the total surface density is exactly balanced by the magnification of faint stars above 
the detection threshold. The chances for the detection of this effect in the Galactic Center appear small. 
A statistically meaningful detection requires a very high surface density that probably exceeds even 
that around the MBH (Wardle & Yusef-Zadeh 1992), and furthermore, models of the stellar luminosity 
function in the inner Galactic Center suggest that s ~ 2 (Alexander & Sternberg 1999). 

Gravitational lensing will also affect the astrometric measurements of stars in the near vicinity of 
the MBH (Jaroszynski 1998a). The bright lensed image of a star on nearly edge-on orbits (within about 
±5°of the line-of-sight to the MBH, as measured in a coordinate system centered on the MBH) will 
exhibit a small astrometric deviation of 0(0.1 mas) in the Keplerian orbit as it passes behind the MBH 
(Nusser & Broadhurst 2004). The astrometric precision required to detect such a tiny effect is well 
beyond present-day capabilities. If measured by future telescopes, the form of the orbital deviation can 
distinguish between different dark mass distributions (§4.2) and be used to break the th/Rq degeneracy 
of astrometric Keplerian orbits (§4.1.2). 



5.2.3 Beyond the point lens approximation 

The point mass lens model and the small deflection angle assumption provide a very useful approxima- 
tion for describing gravitational lensing by the Galactic MBH. Nevertheless, there are several reasons 
to explore more complicated mass distributions and large angle lensing. First, it would be useful if 
gravitational lensing could be used to dispel any remaining doubts that the dark compact mass in the 
Galactic Center is indeed a MBH, and not some other extended distribution of matter (e.g. Capozziello 
& Iovane 1999; Nusser & Broadhurst 2004; §4.2). Unfortunately, it can be shown that he behavior of 
high-magnification light curves near peak magnification is universal and independent of the details of 
the lens (Schneider, Ehlers & Falco 1992, Eq. 11.21b). For spherically symmetric mass distributions 
this implies that the light curves differ only in the low magnification tails, which are much harder to 
observe. The astrometric effects of lensing may provide a more promising way to probe the dark matter 
distribution (§4.2). Second, large angle lensing probes GR in the strong field limit. GR collections to 
the small angle approximation result in a variety of interesting effects (Virbhadra & Ellis 2000; Petters 
2003; Bozza & Mancini 2004), such as "retrolensing", where the Galactic MBH can redirect by #~7r 
the light from a close luminous foreground source (for example the star S2) back at the observer (De 
Paolis et al. 2003). However, these effects are minute and will be extremely hard to observe even by 
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the next generation of instruments. Third, the MBH is surrounded by a massive stellar cluster. Because 
the stellar mass is not smoothly distributed but is composed of discrete point masses, its effect on the 
lensing properties of the MBH is much larger than one may naively estimate by adding the stellar mass 
to that of the MBH (Alexander & Loeb 2001). Enhanced lensing by the presence of discrete stellar mass 
objects may be used to discover SBHs that have accumulated near the MBH due to mass segregation 
(Miralda-Escude & Gould 2000; Chaname, Gould & Miralda-Escude 2001; §3.3). 

The effect of stars on lensing by the MBH is similar to that of planets on microlensing by a star, a 
phenomenon that was studied extensively for the puipose of planet detection (e.g. Gould & Loeb 1992). 
The lensing cross-section of an isolated star is Q\(M±) j Q\(M,) = M*/M, < 10~ 6 times smaller than that 
of the MBH (Eq. 5.12). However, when the star lies near 9 E (M,), the shear of the MBH distorts its 
lensing cross-section, which develops a complex topology, becomes radially elongated and is increased 
by up to an order of magnitude. As stars orbit the MBH, their elongated cross-sections scan the lens 
plane. If these happen to intersect one of the images of a background source that is lensed by the MBH, 
the image will be split into 2 or 4 sub-images whose angular separation will be of order 9e(M+), and 
so the sub-images will not be individually resolved. However, their combined flux will be significantly 
magnified. This will increase the probability of high magnification events over what is expected for 
lensing by the MBH alone. The light curves of such events will no longer be symmetric as they are 
for a point mass, but will exhibit a complex temporal structure, and their typical variability timescales 
will increase sharply for images that lie near 6e(M.) because of the increased stellar cross-section for 
lensing. Enhanced lensing by stars in the Galactic Center is estimated to increase the probability of 
A > 5 lensing events by ~2 and of A > 50 events by ~3 (Alexander & Loeb 2001). 

6 Strong star-MBH interactions 

Some of the most spectacular stellar phenomena that can occur in the GC involve strong interactions with 
the MBH. In this context "strong" is defined as any interaction that cannot be described by Newtonian 
gravity operating on point masses. This could be because the stars' internal degrees of freedom can no 
longer be ignored, or because dissipative interactions have to be included, or because GR effects must 
be taken into account. 

Early interest in such processes focused on the consumption of stars by the MBH, either directly, 
when a whole star crosses the event horizon, or indirectly, when the gaseous debris left after the tidal 
disruption of a star is accreted by the MBH (Hills 1975; Frank & Rees 1976; Lightman & Shapiro 1977; 
Frank 1978; Rees 1988). The interest was driven by the implications for the growth of MBHs, for feed- 
ing luminous accretion in AGN and for detecting MBHs in quiescent galaxies by tidal disruption flares. 
However, tidal disruption is only one of several possible modes of interaction with the MBH. A useful 
analogy can be made between the MBH and the stars around it and an atomic nucleus and its surround- 
ing electrons. In spite of the fundamental differences between a macroscopic classical system and a 
microscopic quantum one, the analogy to atomic processes suggests a classification scheme for the var- 
ious modes of strong star-MBH interactions. This not only provides the "bookkeeping" needed to keep 
track of all the permutations, but also provides some analytical tools for calculating cross-sections, rates 
and branching ratios. The following atomic-like processes, and their possible observational signatures, 
are discussed here. 

Annihilation Direct infall preceded by tidal disruption, where a substantial fraction of the stellar mass- 
energy is released (§6.1). 

Absorption Direct infall where the stellar mass-energy is absorbed by the MBH (§6.1). 



T. Alexander / Physics Reports 000 (2005) 000-000 



57 



Deep inelastic scattering Tidal scattering, where the star is strongly perturbed by the MBH tidal field 
during a close "fly by", but is not destroyed (§6.2.4). 

Metastable decay vs collisional ionization / de-excitation Gradual orbital inspiral due to a dissipa- 
tional force, such as gravitational wave emission or tidal heating, in the presence of orbital per- 
turbations that can throw the star into the MBH or deflect it to a wide orbit (§6.2.1). 

Charge exchange Three-body exchange where an incoming star knocks out one that is on a bound 
orbit to the MBH (the 3rd body) and replaces it (§7.2.4). 

Ionization Three-body encounters where an incoming star eject a star tightly bound to the MBH into a 
wide, unbound orbit (§7.2.4). 

6.1 Tidal disruption 

The tidal interaction between a star and a black hole is characterized by three length-scales: the stellar 
radius R+, the Schwarzschild radius rs and the tidal radius r t . The ratios of these length-scales can be 
expressed in terms of the ratios of the masses and escape velocities, 



where c, the speed of light, is the escape velocity from the BH and Vl-IGM^jR^ is the escape velocity 
from the stellar surface (for a solar type star V e ~600kms~ 1 and (c/V e ) 2 ~2.5x 10 5 ). 

The nature and outcome of a "grazing collision" between a black hole and a star (r p ~rs+R±) depends 
on the relative magnitude of the black hole's tidal field, the stellar self-gravity and its gas pressure and 
whether the disruption occurs in the Newtonian or GR limit (see discussion by Gomboc & Cadez 2005). 
By definition, the tidal energy extracted from the orbit and deposited in the star as it crosses the tidal 
disruption radius is of order of its binding energy, E+. It also follows from the definition of r t that the 
periapse crossing time at the tidal radius, r t , is of the order of the star's dynamical timescale, since 
r t ~ \J rf/Gm = \/Rl/ GM± = r* (r* ~ 1600 s for a solar type star). Therefore, complete tidal disruption 
occurs as fast or faster than the time it takes for the stellar gravity or gas pressure to react. These can 
therefore be neglected during the disruption and the particles composing the star approximately follow 
free-fall trajectories in the gravitational field of the black hole (Laguna et al. 1993; Kochanek 1994). 
This approximation holds until the stream of disrupted gas completes an orbit around the black hole 
and intersects itself. At that point strong shocks likely occur and the gas hydrodynamics have to be 
taken into account. When r p >r t , the tidal interaction is weaker and slower, lasting Tp~T f ft 3 / 2 , where 
b~ l =r t /r p is the "penetration parameter" . It may then become longer then the stellar dynamical time 
scale to the extent that the stellar gravity and hydrodynamics do have to be taken into account. The ratio 
r t jr$ determines whether the disruption occurs in the Newtonian limit (r t /rs> 10) or in the GR limit 
Or As < few). 

The different regimes of tidal disruption can be classified by the mass ratio m/M+ and Eq. (6.1). In 
order of increasing m/M+, these are 

1. m/M+ <C 1, which corresponds to rs <C r, <Ci?*. In this case the tidal interaction is overall weak 
and Newtonian (only a small part of the star overlaps with the tidal disruption zone), and the 
process is dominated by the star's gravity and pressure. The star is not disrupted. If the black 




(6.1) 
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hole passes through the star it may accrete a small fraction of its mass. A more significant effect 
may occur if the black hole is captured inside the star and sinks to the center. This tidal regime 
could apply to the formation exotic stars that cannot be formed in the course of normal stellar 
evolution, such as Thorne-Zytkow objects (Thorne & Zytkow 1975), which are giant envelopes 
with an accreting compact object in their center. 

2. m/M*~ 1, which corresponds to rs<^r t ^R+. In this case a strong Newtonian interaction will 
occur, which would result in significant mass loss and possible disruption, depending on the 
central concentration of the stellar density profile. The process will be accompanied by strong 
shocks. This tidal regime applies to close interaction between a SBH and a massive star (e.g. the 
3-body exchange scenario, §7.2.4). 

3. m/M* ~ (c/V e ) 2 , which corresponds to /?* ~ r$ < r t . In this case the disruption is complete and 
occurs in the Newtonian limit (the subsequent flow of the gaseous debris into the black hole 
is relativistic). This regime applies to disruption by intermediate mass black holes (m ~ 10 3 - 
10 4 M Q ), whose existence is still a matter of speculation (§1.3). 

4. (c/V e ) 2 <m/M*<(c/V e ) 3 , which corresponds to R-),<^rs<r t . Tidal disruption by the MBH in the 
GC falls in this regime (r t ~9r$ ~ 160/?* for a solar type star). The disruption can be treated as 
Newtonian to a good approximation. 

5. m/M* > (c/V e ) 3 , which corresponds to R*<^ r t <C r s . This is the case for a MBH with mass 
> 10 8 M Q (for solar type stars). The event horizon lies outside the (formal) tidal radius and the 
star falls into the MBH as a point particle on a General Relativistic trajectory without being 
significantly perturbed. 

The various processes and phenomena associated with the tidal disruption of stars are reviewed by Rees 
(1988). 

6.1.1 Tidal disruption rate 

A star is tidally disrupted if the MBH mass is small enough so that r t > rs and if its orbital angular 
momentum J is small enough so that r p < r t . Such orbits are called "loss-cone" orbits 18 and in a spherical 
potential have 



where the last approximate equality assumes that r t <C rj, (the potential is dominated by the MBH) and 
e<^(r f ) (the star is consumed from a weakly bound orbit). In terms of Keplerian orbits (Eq. 3.15), 
which are a valid approximation inside r/,, tidal disruption requires that r p - a(l — e) < r t , that is either 
a very tightly bound orbit (small a, high e) or a very eccentric one (high e, low J). The solid angle 
subtended by the loss cone is usually very small. The loss-cone's opening angle for parabolic orbits is 
easily derived from angular momentum conservation (Frank & Rees 1976) 



The term "loss-cone" originates in plasma physics. For motion on a straight line (no gravity), the volume 
in velocity space that corresponds to disruption orbits originating from a point is a cone. With gravity's effects 
included, it is actually a cylinder (Miralda-Escude & Gould 2000). 



J 2 < 4(e) = 2r 2 t ) - e) ^ 2Gmr t 
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where the typical radius is taken as r~a and where 



J 2 = Gma = (Gm) 2 j 2e 



(6.4) 



is the maximal (circular orbit) angular momentum for specific energy e. For example, #/ c ~ 10 for a 
solar mass star on a a = 1 pc orbit around a 3 x 10 6 M Q MBH. 

Loss-cone theory deals with the processes by which stars enter loss-cone orbits, the rate at which 
this occurs, and how these depend on the parameters of the system. This problem has been studied 
extensively over the years (Hills 1975; Frank & Rees 1976; Lightman & Shapiro 1977; Cohn & Kulsrud 
1978; Magorrian & Tremaine 1999; Syer & Ulmer 1999; Alexander & Hopman 2003; Wang & Merritt 
2004). While different approaches and assumptions lead to somewhat different quantitative results, the 
process is well understood qualitatively. 

Any stars that are initially on loss-cone orbits are promptly destroyed on an orbital timescale once 
they reach periapse. From that time on, the tidal disruption rate, F t , is set by the rate at which these 
orbits can be replenished by the relaxation processes that redistribute stars in phase space to high-e or 
low-/ orbits. Diffusion in e-space occurs on the relaxation timescale, which roughly corresponds to the 
timescale for e to change by order unity, 



The relaxation time is approximately independent of energy in the GC (Eq. 3.7), and is assumed here 
to be also independent of /. Since t r is also the timescale for J 2 to change by order J 2 (Eq. 6.4), the 
angular momentum relaxation timescale, which roughly corresponds to the timescale for / to change by 
order unity, is 



where the square root dependence of / on tj reflects the random walk (diffusive) nature of the process . 
Typically, J/ C <^J C . In principle, stars can enter the loss-cone, / <// c (e) either by a decrease in /, or by 
an increase in e (up to the last stable orbit). In practice, diffusion in /-space is much faster: although 
tj~t r when /~/ c , once / becomes smaller than J c by order unity, tj<^it r . Because /-diffusion is so 
much more efficient than e-diffusion, only /-diffusion need be considered. 

The way a star enters the loss-cone depends on A/(e), the r.m.s change in /(e) over an orbital 
(dynamical) time P(e), 



The ratio A/(e)/// c defines two dynamical regimes of loss-cone re-population (Lightman & Shapiro 
1977). In the "Diffusive regime" of stars with large e (tight short-period orbits), A/(e)<C// f and so the 
stars slowly diffuse in /-space. The loss-cone remains nearly empty at all times since any star inside 
it is promptly destroyed. At / 3> // c the DF is nearly unperturbed by the existence of the loss cone. 
In particular, it can remain very nearly spherically symmetric and depend only on e (§3.2), but it falls 
logarithmically to zero at / > // c . In the "full loss-cone regime" (sometimes also called the "pinhole" 
or "kick" regime) of stars with small e (wide long -period orbits), A/(e)3>// c and so the stars can enter 

19 In fact, it is the velocity component perpendicular to the velocity vector, Avj_, that executes a symmenic 
random walk, and not directly / itself. This introduces a drift term (oc t) to the evolution of /, which is small 
compared to the diffusive term (ex y/t) as long as t <tj and can be neglected (Hopman & Alexander 2005). 



t r ~e/e. 



(6.5) 



tj~J 2 /(j2) ~[J/J c (e)] 2 tl 



(6.6) 




(6.7) 
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and exit the loss-cone many times before reaching periapse. As a result, the DF is nearly isotropic at all 

J>Jlc 

The tidal disruption flux density (stars disrupted per dtde) in the diffusive regime, f^if, is derived 
from the orbital diffusion coefficients for two-body scattering via the Fokker-Planck equation, with the 
boundary conditions that the stellar density fall to zero for J^Ji c and approach the isotropic distribution 
«iso(e) for J— >J C (Lightman & Shapiro 1977), 

ln[/ c (e)// fc (e)]f f .(e) 

The diffusion flux is regulated mainly by the relaxation time and is only weakly dependent on the size 
of the loss-cone (typically 5 < ln(/ c /7/ c ) < 10) 20 . The logarithmic term reflects the fact that the diffusion 
occurs in a 2D space of solid angle (Lightman & Shapiro 1977). It can also be understood qualitatively 
by analogy to heat diffusion on the surface of a sphere from the equator ("7 C "), to a small hole at the 
pole ("7 /c ") (Frank & Rees 1976). 

In the full loss-cone regime far from the MBH, where the DF is nearly isotropic, the tidal disruption 
flux density is simply oc7n? 2 L. Since the loss-cone is refilled every orbital time, the flux density is 



Jf c (e) n iso (e) 
JKe) P{e) 



F V in(£;Jic) ~ 5; A Z\ ■ (6-9) 



The transition between the two regimes occurs at a critical energy e c (or equivalently a critical 
typical orbital radius r c ), where Fdif(e c )=Fpin(e c ), or approximately where A/(e)~// c . The total tidal 
disruption rate is then given by 



rsc roc 

r,~ / F pin (e;J k .)de+ F m {e;J lc )de . (6.10) 

J0 he 



An analysis of the behavior of the integrands in Eq. (6.10) for typical observed values of stellar cusp 
slopes, reveals that most of the contribution comes from the diffusive regime at values near e max ~ 
max(e f ,£/,), where £/, ~V( r /j)> or equivalently, near r max = min(r c ,r/j) (Lightman & Shapiro ; Syer & 
Ulmer 1999; for numeric examples see Magorrian & Tremaine 1999). In most galaxies r c < r/, (Wang 
& Merritt 2004). Making use of this property, the tidal disruption rate can be approximated to within an 
order of magnitude by 



r, N ^l h) . (6.ii) 

\a^r h /r t )t r (r h ) 

Taking for the GC r h ~2 pc, r, ~ 10 13 cm, JV*« r h )~m/ (M*) ~3.5 x 10 6 and t r (r h ) ~4x 10 9 yr (Eq. 3.8), 
the estimated rate is T t ~ 10~ 4 yr" 1 , in agreement with more detailed calculations (T t ~ 5 x 10~ 5 yr -1 , Syer 
& Ulmer 1999). 

It should be noted that there are dynamical mechanisms that may increase the tidal disruption rate 
beyond that predicted by two-body scattering in a spherical system. The possibilities include the effects 
of massive perturbers, for example giant molecular clouds (Zhao, Haehnelt & Rees 2002), which indeed 
exist in the GC on scales of 1-2 pc from the MBH; enhanced rates due to resonant scattering (factor of 
~2 rate increase; Rauch & Tremaine 1996; Rauch & Ingalls 1998); deviations from spherical symmetry 
(factor of ~2 rate increase; Magorrian & Tremaine 1999); or chaotic orbits in triaxial potentials (factor 
of 10-100 rate increase; Poon & Merritt 2002; Merritt & Poon 2004). 

20 The diffusion timescale is longer than ?,. by the logarithmic factor because the stellar phase-space density 
approaches zero near the loss-cone. 
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6.1.2 Tidal disruption and its aftermath 

The observations of tidal disruption in the GC are realistically limited to those after-effects that persist 
for at least ry 1 ~ 10 4 yr (the mean time between events). As described below, there are two main 
channels by which energy is released in a tidal disruption event: one is a luminous accretion flare lasting 
a few years, and the other, which is much less luminous but longer-lasting, is a shock wave, analogous 
to a supernova blast, that sweeps through the interstellar medium and can energize it substantially for 
~ 10 5 yr (e.g. Khokhlov & Melia 1996). Observationally, the two effects complement each other. While 
tidal disruption flares may be, and perhaps were already observed in cosmological galaxy surveys (e.g. 
Komossa 2002), the supernova-like shock wave may be observed in the GC. 

The hydrodynamics, radiative properties, early and late evolution of tidal disruption in the New- 
tonian and relativistic limits were studied extensively by numerous authors, both analytically (Rees 
1988; Phinney 1989; Cannizzo, Lee & Goodman 1990; Kochanek 1994; Loeb & Ulmer 1997; Ulmer, 
Paczyriski & Goodman 1998; Ulmer 1999; Menou & Quataert 2001; Ivanov & Novikov 2001; Ivanov, 
Chernyakova & Novikov 2002) and by simulations (Nolthenius & Katz 1982; Evans & Kochanek 1989; 
Khokhlov, Novikov & Pethick 1993; Frolov et al. 1994; Laguna et al. 1993; Diener et al. 1997; Ayal, 
Piran & Livio 2000; Bogdanovic et al. 2004; Gomboc & Cadez 2005). While there is considerable 
uncertainty about details, a general picture has emerged. 

It is useful to begin by considering the energy budget of a tidal disruption event. The star starts on 
a nearly parabolic orbit (Eorb — 0) and so the total energy of the system (ignoring for now its rest-mass 
energy) is the stellar binding energy, < 0. The outcome of tidal disruption is that some energy is 
extracted out of the orbit to unbind the star and accelerate the debris. Initially about Mj„ ^M+/2 of the 
stellar mass becomes tightly bound to the MBH with energy E in , (-E m S> while the remainder 
A^out ~^*/2 of the stellar mass is forcefully ejected with energy E out , (+E 0Ut S> so that E[ n +E out = 
£*. 

There are 3 distinct energy scales in the process, each associated with a different time scale. The 
lowest is the stellar binding energy By virtue of the virial theorem, the internal heat of the star is of 
order U ~- E+/2. During the disruption, on the periapse passage timescale r p ~ 0(1 hr), the inner parts 
of the star are exposed as it is stretched and the heat can be radiated away. This may lead to a short flare 
immediately after disruption (Gomboc & Cadez 2005). 

Next in scale is the kinetic energy of the ejecta. While the mean specific energy of the debris is 
—E+/M* ~ (V* is the circular velocity on the stellar surface), the spread about this value can be 
orders of magnitude larger. As the star is disrupted, tidal torquing spins it up in the direction of the 
orbit, imparting an excess velocity > V e = v2V* above the mean orbital velocity v p 3> V* to the stellar 
hemisphere on the far side from the MBH, and a comparable velocity deficit on the near side. The 
difference in velocities leads to a spread in the specific kinetic energies of about (v p +V e ) /2-(v p -V e ) /2= 
2 3 / 2 v p V* 3> V? (Rees 1988). The spread in energy can also be expressed in terms of the work done on 
the stellar debris by the tidal force (Lacy, Townes & Hollenbach 1982), 

GmM / m \ '/ 3 

where A0mbh(t>) is the MBH's gravitational potential difference across the star at periapse. In a typical 
disruption event in the GC, W > lOOZs* (for Z?~ 1) and thus the debris is released with specific energies 
in the range (E+±W)/M+. Note that while W it is still much smaller than the orbital kinetic 

energy of the star at disruption, -b~ x (m / 'M*) 2 / 3 'E*. The kinetic energy carried by the ejected debris can 
significantly exceed that released by a normal supernova (~ 10 51 erg) if the orbit is highly penetrating, 
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Unlike a spherical supernova explosion, the mass will be ejected mostly in one half of the orbital plane 
(Khokhlov & Melia 1996). 

The highest energy scale is that due to accretion onto the MBH. The bound debris will eventually cir- 
cularize around the MBH after the stream of orbiting gas will cross itself, collide and shock. In order to 
fall into the MBH, the gas has to lose its angular momentum. If a mechanism exists that can exert torque 
on the rotating gas (for example shear due to effective viscosity mediated by magneto-hydrodynamic 
turbulence), it will gradually spiral inward on a viscous dissipation timescale. The gravitational binding 
energy of the gas to the MBH will be converted to heat, until it reaches the last stable circular orbit 
around the MBH (3r$ for a Schwarzschild BH), where it will promptly fall into the event horizon on a 
dynamical timescale. If the gas can dissipate the heat by the emission of radiation, or acceleration of 
particles, then between r] = 0.057 (for a non-rotating BH) to ?] = 0.42 (for a maximally spinning BH) of 
its rest mass energy can be extracted by the time it reaches the last stable orbit (e.g. Shapiro & Teukol- 
sky 1983). If, on the other hand, energy dissipation is inefficient, the heat may be carried by the flow 
(advected) into the MBH and add to its rest mass (advection dominated accretion flow [ADAF] solu- 
tion; Rees et al. 1982, Abramowicz et al. 1995; Narayan & Yi 1994, 1995; Narayan, Yi & Mahadevan 
1995). Alternatively, the gas will not fall into the MBH but will escape as a wind (advection dominated 
inflow-outflow solutions [ADIOS], Blandford & Begelman 1999) or by convection outward (convection 
dominated accretion flow [CDAF], Quataert & Gruzinov 2000; Narayan, Igumenshchev & Abramowicz 
2000). In the case of high radiative efficiency, 

E acc ~ri(M*/2)c 2 ~rib 2 ( — -J f^j E 0Ut >E 0Ut . (6.14) 



Immediately after the disruption and acceleration at periapse, the debris' self gravity and gas pres- 
sure are dynamically insignificant compared to its kinetic energy, and the gas streams follow ballistic 
(free falling) Keplerian orbits. The orbits are very eccentric l-(e) = r p /a = 2b(M*/m) 2 / 3 with a large 
spread in orbital period and apoapse, extending from P — > oo (r a — > oo) of the barely bound gas, down to 
a minimal period P mm , (r„ ^^(m/M^) 1 / 3 ), which corresponds to that of the most tightly bound gas with 
specific binding energy b^im/M^^V 2 (Eqs. 6.12, 3.15; Rees 1988; Evans & Kochanek 1989). The 
minimal period is 

. m /M* \ 1/2 / ^ \ 3/2 / M* \ " 1/2 , 
*V*~0.1yr -rsr TT 1 ) b - ( 6 - 15 ) 




M* J \ 10 6 J \RqJ V M e 



After an initial delay of P m i n the debris will start returning to the point of disruption at a rate M = 
(dM/d/s) (dE/dP). Simulations show that the spread in energy of the bound debris is nearly uniform, 
AM/AE ~ const (e.g. Evans & Kochanek 1989; Ayal, Piran & Livio 2000), and so (Eq. 3.15) 

M+ ( t \" 5/3 

M in (t)~—^[ — ) , (?>P min ). (6.16) 

When r p > r$, GR correction have to be taken into account. These lead to a stronger tidal interaction, 
and consequently a higher spread in energies (Eq. 6.12), a shorter P m [ n and larger peak debris return rate 
(Ulmer 1999). 
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Initially, the mass return rate exceeds the Eddington-limited accretion rate 21 

In the GC, maxM~0.1M Q yr -1 . The return rate remains higher than the Eddington-limited rate for no 
more than a few years, 



-K^) 2/5 ^ /5 (") 2,5 fe) 2/ --(^) 3/5 ^(#)" 2/5 (^) 2/5 - 

(6.18) 

where TEdd = wj/^Edd is the e-folding timescale of Eddington accretion with unit efficiency. A similar 
estimate is obtained by considering the time for accretion of all the debris at the Eddington-limited 
rate if all the mass returned at once, fedd ~M,„/maxM~5yr. The high initial mass return rate does 
not necessarily imply an early phase of Eddington-limited accretion. In order for the disruption event 
to produce a luminous flare, both the time for circularization, which sets the flow pattern required for 
dissipation, and the timescale for viscous dissipation, which extracts the angular momentum, have to be 
short enough so as not to be an impediment to the accretion. This indeed appears to be the case. 

Circularization is expected to be efficient. The evolution of the bound debris is driven by shocks 
between streams of gas that converge to their common point of origin at periapse. Relativistic preces- 
sion (Eq. 5.6) will also lead to shocks between incoming and outgoing streams. These shocks will 
redistribute angular momentum and energy between the gas streams. Some fraction of the initially 
bound debris will become unbound after undergoing strong shocks and heating on the second periapse 
crossing, leaving perhaps ~M*/4 of the gas bound to the MBH (Ayal, Piran & Livio 2000). While 
the exact details are uncertain, it is expected that circularization will occur rapidly, on a timescale of 
? circ ~fewxP min (Ulmer 1999). 

The viscous accretion timescale at the tidal radius is t acc ~ 2ttt { /airh 2 = 3 x lO~ 4 (r^/r ltj 0)/a7r/i 2 yr, 
where a is the dimensionless viscosity parameter (the proportionality factor between the viscous stress 
tensor and the pressure) and h is the ratio of the thickness of the accretion disk to its radius (h ~ 1 for 
a thick disk). When M^Mgdd, the radiation pressure puffs up the disk (h— > 1) and decreases t acc . The 
value of a is uncertain, but 10~ 3 <a< 10 _1 seems plausible 22 (Ulmer 1999), so that probably ? acc <0.01 
yr. It then follows that t acc < ? c i rc < ?Edd> which means that for the first few years after a disruption 
event the accretion will probably proceed at a super-Eddington rate (L^L^dd, small r\, M>Mgdd) with 
a radiatively inefficient accretion flow, such as ADAF, ADIOS and/or CDAF (Rees 1988; Menou & 
Quataert 2001). 

Most of the work on the observational signature of tidal disruption concentrated on "tidal flares", 
which correspond to this initial luminous phase, lasting months to years. Tidal flares are expected to 
peak in the UV/X-ray range, as can be estimated by assuming that the Eddington luminosity is emitted 



21 The luminosity released by matter accreting on a mass m exerts an outward radiation pressure on the infalling 
mass. The radiation pressure balances gravity and halts the accretion when the luminosity reaches the Eddington 
luminosity, LEdd = 1-5 x 10 38 (ra/MQ)ergs -1 (for isotropic luminosity and pressure due to Thomson scattering). 
The Eddington accretion rate M-E&d = £Edd/c 2 = 3 x lO~ 9 (m/Mo)M0yr _1 is defined as the maximal possible for 
unit efficiency. The actual Eddington-limited accretion rate is maxM=MEdd/f?- 

22 Alternatively, if 10~ 5 <a< 10~ 3 and rj > few x 0.01, the configuration may be analogous to a Thorne-Zytkow 
object, in the form of a large (~ 100r f ) optically thick star-like object with 10 4 K, powered by accretion. Such 
an object will emit at the Eddington limit mainly in the optical (Loeb & Ulmer 1997), unlike the UV/X-ray flares 
discussed here. 
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as blackbody radiation from a surface whose typical size probably lies between the last stable orbit and 
the tidal radius. The blackbody temperature is then rbb = (LEdd/47rr 2 cJsB) 1 ^ 4 ) and 

^(r t )-3xl0 5 (^_y i2 K, or r bb (3r 5 ) ~ 7 x 10 5 (-^-) (6.19) 

(Assuming a solar-type star for r t ), almost independently of the MBH mass. This temperature corre- 
sponds to a blackbody peak at ~ 0. 1 keV. 

There are a few tentative detections of tidal flares in other galaxies (see review by Komossa 2002). 
However, the short duration of the flare and the mean time between tidal disruption events in the GC 
(§6.1.1) translate to a very low duty cycle of ~ 10 -4 . The probability of directly observing a tidal flare 
in the Galaxy is small. 

A few years after the tidal disruption event, when the return rate drops below the Eddington rate, 
the gas will cool efficiently and accrete via a radiatively efficient, geometrically thin accretion disk, at a 
rate set by the viscous timescale. The accretion rate in the thin disk is expected to fall off slowly with 
time (Cannizzo et al. 1990). However, this result is at odds with observations. The predicted average 
luminosity due to accretion of debris from a past tidal disruption event is 2-4 orders of magnitude higher 
than that observed in the GC, or in the centers of nearby galaxies M31 and M32 (Menou & Quataert 
2001; §6.2.5). This discrepancy has been interpreted as indicating that most of the debris was blown 
away during the early super-Eddington accretion phase. Following that, the luminous accretion phase 
shuts off early, either because the low-mass accretion disk spontaneously switches over to a low radiative 
efficiency solution when M< 10~ 2 MEdd, or else because increasingly larger parts of the cold outer disk 
become neutral, unable to sustain magneto-hydrodynamic turbulence and therefore effectively inviscid 
(Menou & Quataert 2001). 

6.1.3 Signatures of tidal disruption 

While there is as yet no clear evidence of tidal disruption in the GC, there are several ideas about the 
form such evidence could take. 

The ejection of tidal debris with energies that far exceed those released in normal supernovae could 
provide a long-term signature of tidal disruption, in the form of unusual supernova remnant-like struc- 
tures. This was suggested as the origin of the Sgr A East radio source (Khokhlov & Melia 1996). Sgr A 
East is a prominent non-thermal radio continuum source with an elongated shell-like structure extend- 
ing over 10-20 pc, offset from Sgr A*, but engulfing it at the edge of the structure (e.g. Yusef-Zadeh & 
Morris 1987). The shell structure and the non-thermal spectrum are reminiscent of a supernova remnant. 
Early estimates of the energy required to evacuate the radio-emitting cavity inside the dense surrounding 
molecular clouds yielded an extreme lower bound of E > 4 x 10 52 erg (Mezger et al. 1989). The high 
energy requirement, large size and non-spherical shape are in line with what is expected from a deep 
tidal disruption event. However, more recent estimates of the energy from X-ray observations (Maeda at 
al. 2002; Sakano et al. 2004) yield a much lower value, which together with the metal -rich abundances 
observed in Sgr A East, are consistent with the energy and elemental output of a single normal core 
collapse supernova (but see §6.1.4 about the possibility of nucleosynthesis during a tidal detonation of 
a low mass star). The single supernova interpretation is also supported by observations of molecular 
transitions, which allow an estimate of the kinetic energy in molecular gas that is being swept by the 
expansion of Sgr A East (Herrnstein & Ho 2005). 

A tidal disruption event is expected to produce a short luminous UV/X-ray accretion flare (§6.1.2), 
which will likely also have a hard X-ray component, as is seen from accreting MBHs in AGN. These pri- 
mary X-ray photons will encounter gas clouds as they propagate from the MBH, and will be Thomson- 
scattered into the line-of-sight or absorbed and re-emitted as fluorescence lines, primarily in the 6.4 keV 
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Ka line from neutral Iron ("X-ray reflection nebula"). The reflected photons arriving at the observer at 
any given time were scattered (or re-emitted) from a parabolic surface around the MBH. Reflected X-ray 
emission is not an exclusive signature of tidal disruption. There could well be other causes for episodic 
accretion flares or short high-state epochs. However, the spatial distribution and temporal changes in 
the scattered X-ray continuum photons and the emission lines can at least provide information on the 
position of the source and the time and duration of the flare, as the paraboloid of equal arrival time 
sweeps the volume around the MBH at the speed of light (Sunyaev & Churazov 1998). The X-ray re- 
flection effect could effectively lengthen by ~ 2 orders of magnitude the time available for detecting an 
accretion flare. However, even if the distribution of gas clouds extends far from the center, it is unlikely 
that the fluorescence could be detected as long as T7 1 ~ 10 yr years after the tidal disruption event (< 1 
kpc from the center) due the geometrical dilution of the primary X-ray flux with distance. 

The present observed level of diffuse X-ray emission from the central ~ 700 lyr of the GC was used 
to argue that Sgr A* has not radiated at the Eddington level for even a day in the last 400 yr (Sunyaev, 
Markevitch & Pavlinsky 1993). The detected 6.4 keV line emission from several sites near the GC 
(Koyama et al. 1996; Murakami et al. 2000, 2001) was found to be well reproduced by an X-ray 
reflection nebula model. The luminosity required to reproduce the 6.4 keV line exceeds by 3 orders 
of magnitude any observed X-ray point source in the region, and by 6 orders of magnitude the present 
average X-ray luminosity of SgrA* (Baganoff et al. 2001; Goldwurm et al. 2003) . This was taken as 
evidence that Sgr A* was active > 300 yr ago. The detection of high energy X-rays (< 200 keV) from 
a giant molecular cloud further supports this interpretation (Revnivtsev et al. 2004). The reconstructed 
X-ray luminosity and hard X-ray spectrum slope are characteristic of low-luminosity AGN, and suggest 
that Sgr A* was active for at least 10 yr some 300-400 years ago. However, recent X-ray observations 
with high spectral resolution (Predehl et al. 2003) do not support this interpretation, due to the absence 
of a strong 7.1 keV absorption edge that is expected to accompany the 6.4 keV emission line in X-ray 
reflection nebulae, and due to the lack of the expected inverse correlation between the 6.4 keV line 
surface brightness and distance from the center. 

Another possible long-term after-effect of a strong tidal encounter is tidal stripping, the partial 
disruption of evolved giants (Davies & King 2005). Evolved giant stars have an extended, cool, low 
density envelope and a hot, massive compact core. Because of its low surface gravity, the envelope 
can be stripped off the core even at large periapse distances where denser main-sequence stars are safe. 
The bare core, while no longer burning, will likely remain hot (T > 10 5 K) and luminous (L > 10 2 L Q ) 
for 10 3 -10 6 yr after the tidal stripping event. The hot core will appeal - as a "supersoft X-ray source", 
which can be detected in the GC and in the centers of other nearby galaxies (Di Stefano et al. 2001). 
The challenge lies in distinguishing tidally stripped cores from normal young hot white dwarfs and the 
He cores of massive stars that lost their envelopes by strong winds or by interactions with a binary 
companion. In addition, Collisional stripping (§3.4.1) is expected to produce a similar effect. 

6.1.4 Tidal detonation 

When a disrupted star crosses deep inside the tidal radius, b~ l >5, it experiences a very strong and rapid 
tidal compression during periapse passage. Figure (6.1) shows a qualitative demonstration of the effect 
by a simple (low-resolution Newtonian) Smoothed Particle Hydrodynamics (SPH) simulation 23 . The 

23 SPH is an algorithm for simulating the hydrodynamics of 3D self-gravitating fluids, which is commonly used 
in the study of stellar collisions (Lucy 1977; Gingold & Monaghan 1977; see review by Monaghan 1992). The 
star is represented by discrete spherical mass elements, each with an internal density distribution that peaks in the 
center and falls to zero at the edge. The total density at a point is the sum of densities in all overlapping spheres 
that include the point. The resulting density field is continuous and differentiable, and so its thermodynamic 
properties can be evaluated everywhere once an equation of state is specified. Every time step, the positions of 
the mass elements are updated according to the gravitational force and the pressure gradient, and the sphere sizes 
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compression occurs perpendicular to the orbital plane and it more than compensates for the stretching 
of the star along the orbit, resulting in a rapid increase in the central density. This happens despite the 
fact that the tidal force tensor is traceless (i.e. the force lines do not converge) — Both the inward- and 
outward-directed tidally induced velocities relative to the star's center, V t ~/?*/r p , are of the same order 
of magnitude. However, expansion along one dimension leads to a relative increase in that dimension 
of only order unity, whereas collapse in another dimension leads to an arbitrarily large relative decrease 
in that dimension, so the net effect is a vanishing volume with diverging central density (i.e., on the 
periapse passage timescale r p the volume shrinks roughly as Rl — ► R*(R i < + V t T p )(R+- V t r p ) — ► 0). An 
order of magnitude estimate indicates that the maximal central compression ratio can be as high as 
maxpr/ p c ~ b~ 3 , where p c is the initial stellar central density. This upper limit is derived by assuming 
the full conversion of the kinetic energy of the adiabatically collapsing stellar ellipsoids into internal 
heat (Carter & Luminet 1983). The actual compression ratio that can be achieved, and its ultimate 
consequences depend on the uncertain details of the response of the gas hydrodynamics and geometry 
to the sudden compression via increased nuclear energy production and shocks. Simulations indicate 
that the maximal compression ratio is significantly smaller than the kinematic upper limit (Luminet & 
Marck 1985; Khokhlov, Novikov & Pethick 1993), but still very high, max p c /p c b~ 2 (Bicknell 
& Gingold 1983; Laguna et al. 1993). 

It has been suggested that such a strong tidal compression could lead to a ignition of nuclear reac- 
tions and the sudden release of energy, thereby producing a supernova-like explosion and an energetic 
outflow of metal-enriched stellar debris (Carter & Luminet 1982; 1983; Carter 1992). The nuclear en- 
ergy produced during the compression is E mc ~ Qr p , where Q is the nuclear power. The nuclear power 
from a tidally compressed burning core scales very roughly as 2~L^ (p c /p^\ k (r c /r*) n , where L*, p* 
and T* are the initial equilibrium luminosity, central stellar density and temperature. For hydrogen burn- 
ing, k=\, n~4 for the p-p chain and «~ 16 for the CNO cycle; for helium burning (triple-a), k-2 and 
«~40. As a crude estimate, consider a simplified scenario where the compressed stellar core continues 
burning by the same nuclear process, but with scaled-up power. Assuming adiabatic contraction of an 
ideal gas (i.e. neglecting nuclear burning during the contraction), with Poc p 7 and T oc p 7_1 for 7=5/3, 
including the scaling of the periapse passage time, t p ~T*b 3 ' 2 , and parameterizing max p c /p c ~ b~ u 
(u=3/2 will be assumed below), the nuclear energy produced varies with b as 

E nuc ~ n L*b- s , ff = «[*+n(7-l)]-3/2. (6.20) 

For a solar type star (p-p burning, E± ~4x 10 48 erg and r^ = 1600 sec, s= 10), the energy produced 
is negligible, E BUC /E+~2xlQT Z?~ 10 <C 1, as long as b~ l < 13. The p-p chain is too slow to be an 
important source of energy during the brief tidal compression. However, for a 10M Q main sequence 
star, such as are observed near the Galactic MBH (CNO cycle burning, £ t ~8 x 10 49 erg, ~4800 sec 
and L± ~6x 10 3 L Q , s= 16), E/E* ~ 10~ 9 Z?~ 16 ~ 1 even for a moderate penetration parameter of b~ l ~3.6. 
CNO cycle burning accelerated by the tidal compression will cause the star to explode while it is being 
disrupted. A helium burning giant, where E^och -415 , will be even more susceptible to tidal detonation. 

Since the cross-section for a close encounter with the MBH with periapse < r p scales as r p (§6.2), 
tidal detonation events are expected to occur at an even smaller rate than that of tidal disruption events, 
r^~r^< 10~ 5 yr (§6.1.1), and so it is not likely that such an event will be observed in the GC in 
real time. Moreover, the tidal detonation is energetically negligible compared to the kinetic energy that 
is transferred to the debris by the tidal field from the stellar orbital energy (§6.1.2). However, tidal 
detonation can affect the aftermath of a tidal disruption event more subtly. By increasing the velocity at 
which debris is expelled relative to its center of mass along the orbital plane, it can modify the properties 



are readjusted to reflect the changes in the local space density of elements. 
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of the subsequent accretion on the MBH by making the bound matter more tightly bound and strongly 
relativistic. An additional effect of "pancake detonation" is to modify the chemical composition of the 
interstellar medium around the MBH. Pichon (1985) finds that nucleosynthesis by pancake detonation 
leads to the synthesis of proton-enriched isotopes, in contrast with core-collapse supernovae, which 
synthesize neutron-enriched isotopes. Unusual nucleosynthesis may also be associated with a related 
process, the tidal detonation of a WD caused by a combination of tidal compression and enhanced 
self-gravity due to GR terms. The compression can lead to pyconuclear reactions and a thermonuclear 
runaway under highly degenerate conditions, even when the WD mass is well below the Chandrasekhar 
mass (Wilson & Mathews 2004). 

6.2 Dissipative interactions with the MBH 

A star whose orbit passes just outside the tidal disruption radius, or the event horizon if r$ > r t , will 
survive the close encounter with the MBH. Such "near misses" often involve some orbital energy loss, 
AE, for example by the emission of GW (when rs < r p < fewxrs) or by the work invested in raising 
stellar tides and exciting stellar oscillations (when r t < r p <few xr t ). If the star maintains its nearly radial 
trajectory for many orbital periods, the dissipational energy losses will gradually accumulate and cause 
the orbit to shrink and the star to spiral in. This will continue until either the dissipation shuts off, or 
the star is destroyed by the MBH, or is disrupted by the dissipational energy. An initial close encounter 
with the MBH is in itself not enough to guarantee that the star will reach the final stages of inspiral. The 
inspiral timescale is usually orders of magnitude longer than the orbital timescale. During that time the 
star is vulnerable to the same scattering process that deflected it into the low-/ orbit in the first place. It 
can be scattered again to an orbit with a larger periapse, where dissipation is inefficient, or to a loss-cone 
orbit where it is promptly destroyed. Thus the probability for completing inspiral is small, and the rate 
of inspiral events is much smaller than that of direct infall events (§6.2.1). In the GC r t >r$, and most 
stars that undergo a non-disruptive close tidal encounter with the MBH are subsequently scattered to 
a wide orbit and survive after experiencing an extreme tidal distortion, spin-up, mixing and mass-loss 
that may affect their evolution and appearance. Such "tidally scattered" stars eventually comprise a few 
percents of the stellar population within the MBH radius of influence (§6.2.4). 

If the inspiraling star avoids being scattered to a wider orbit or falling directly into the MBH, its 
orbital period P will decrease every peri -passage and the mean dissipated power, ^AE/P, will increase. 
If the emitted power grows high enough, it may have observed in the form of GW emission (§6.2.2), 
or high tidal luminosity (§6.2.3). A different class of mechanisms for orbital energy dissipation (not 
relevant for the present-day GC) is by drag against a massive accretion disk, if one is present (§6.2.5). 
This may play a role in feeding stars to MBHs. 

6.2.1 Orbital inspiral into the MBH 

Inspiral can be treated as an extension of loss-cone theory. The inspiral problem was studied in the 
context of GW emission from compact objects orbiting a MBH (Hils & Bender 1995; Freitag 2001; 
Sigurdsson & Rees 1997; Sigurdsson 1997; Ivanov 2002; see review by Sigurdsson 2003; Hopman & 
Alexander 2005), and was generalized to arbitrary dissipation mechanisms and applied to tidal capture 
by Alexander & Hopman (2003) and Hopman, Portegies-Zwart & Alexander (2004). 

Inspiral processes are characterized by an inspiral timescale, to(eo,Jo), which is the time needed for a 
star starting at an orbit with initial energy so and angular momentum To to complete the inspiral (e— > oo, 
P — > 0). The dissipated specific energy per orbit, Ae, typically increases sharply as r p is decreased 
and correspondingly to decreases sharply with r p (Eq. 6.29). Over the range of interest the inspiral 
time satisfies P<Oo<0, (when P<to, r p is so small that the inspiral is indistinguishable from direct 
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Figure 6.1: A sequence of snapshots from an SPH simulation of a M* = \0M & , R*=4.5Rq star (rep- 
resented by an n = 1.5 ideal gas polytrope) being disrupted by a very close passage to a 3x 10 6 M Q 
MBH on a parabolic orbit with r p = 4R* = 0.06r, = \Ars- The simulation assumes Newtonian gravity 
and neglects nuclear burning, and so while the results are qualitatively correct, they are quantitatively 
inaccurate: GR effects analogous to periapse shift, omitted here, will actually cause the expanding de- 
bris in penetrating encounters (b~ l 3> 1) to fan out in a thin crescent-like shape centred on the MBH 
(Laguna et al. 1993). Top left: face-on view of the star as it orbits (counter clockwise) the MBH (dot at 
bottom right corner). Top right: edge-on view of the same orbit (the two left-most snapshots shown in 
the face-on view are omitted from the edge-on view, for clarity). Bottom: time evolution of the central 
density. Time is measured in terms of the star's dynamical time t* = 4840s. Note the strong and rapid 
tidal compression at periapse (a density increase by a factor of 1 1 on a timescale of ~ 100 s) . 



infall. Conversely, when t$<t r , r p is so large that dissipation is negligible). The condition to<^t r also 
implies that e-diffusion by scattering is negligible during the inspiral. The same is not generally true for 
/-diffusion. However, when the initial inspiral time is also much shorter than the angular momentum 
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relaxation timescale (Eq. 6.6), 



tQ(Eo,JoXtj(eo,Jo), (6.21) 

then the inspiral is completed with /~/o nearly constant. When Eq. (6.21) is satisfied the star decouples 
from the random perturbations that governed its dynamics in the initial scattering-dominated phase, and 
it enters the final dissipation-dominated phase where its orbit is deterministic and governed by the MBH 
potential and the dissipation mechanism only. 

It is not necessary to specify / to obtain a statistical description of t . The mean inspiral time t 
of stars that complete the inspiral is effectively a function of eo only because as long as tj <?o < t r , 
/-diffusion by scattering will continue until such time when / satisfies t <s^tj. Thus, the requirement 
of successful inspiral, together with the dynamics of random walk in /-space in the presence of dissi- 
pation, define a mean value of angular momentum in the inspiral phase, Jq{eq). The value of Jq and the 
distribution of Jq around it can be calculated once the specific properties of the stellar system and the 
dissipation mechanism are given (Hopman & Alexander 2005). 

In analogy to the case of direct disruption, Eq. (6.21) sets a lower limit, e c , on the initial energy 
(or an upper limit, r c , on the initial orbital radius) that is required for completing inspiral. Stars with 
so > s c (the diffusive regime) will typically complete inspiral, whereas stars with eo < £ c will be typically 
scattered to high-/ or loss-cone orbits and therefore have a vanishingly small chance of completing 
inspiral. Because t ^>P, the phase-space volume available for the diffusive regime of inspiral is much 
smaller than that of direct infall ( e J, ns P iral ^>ej, nfa11 , r J, ns P iral <^rj, nfa11 ) and consequently the inspiral rate is 
generally much smaller than that of direct infall. Inspiral can be likened to a race between scattering 
and dissipation, in which dissipation can win only if the track is short. 

It is sometimes the case that the energy dissipated per orbit, Ae, is constant over the course of the 
inspiral to a good approximation (for example when Ae depends mainly on the periapse distance, which 
remains nearly fixed when to<^tj and when the dissipated angular momentum A/ is small). In that case, 
the inspiral timescale in a Keplerian potential, to (calculated relative to a given time based on the values 
of Ae, eo and Pq at that time), is 

i r°° e P 1 

Ae Jeo Ae £ V2 Ae 

and the period and the semi-major axis evolve simply as 

P = Po(\-t/to)\ a = ao(\-t/to?. (6.23) 
The number of inspiral orbits by time t < to is then given by 

to 



Nmb= 2P 



t 

(i--r 2 -i 



to 

6.2.2 Gravitational wave emission 



(6.24) 



A compact enough object that can survive tidal disruption at a distance of r p ~fewxrs, will dissipate its 
orbital energy and angular momentum by the emission of GW radiation. As the orbit decays, both the 
emitted power and frequency of the GW increase. The specific orbital energy and angular momentum 
(per unit stellar mass) lost each orbit by the emission of GW are (r p > 3r$, m^M*) (Peters 1964) 

A £GW = -^/(e)^P) , f«)= f flf ■ (6-25) 
5V2 m \r s J (l+e) 1 ' 2 
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AJ GW = -—g(e) - , ^)=7T^2- (6-26) 

5 c \r s J (l+e) 1 

An orbit in the Schwarzschild metric ends in the MBH if its angular momentum is smaller than 
some threshold, which for bound orbits is almost independent of e. The effective loss cone is therefore 
defined as the threshold angular momentum for a zero-energy orbit (e.g. Shapiro & Teukolsky 1983) 

AGm 

Jlc = • (6.27) 

c 

Stars begin their inspiral on extremely eccentric orbits. The interplay between dissipation and scattering 
in the presence of a mass sink conserves the trend toward very high eccentricity until the very last stages 
of the inspiral (Hopman & Alexander 2005). For e~ 1 the periapse can be expressed as r p /r$ = 4(7/7/ c ) 
and 



AEGW = £| (_j , ^3^—. (6.28) 
The GW inspiral time increases rapidly with /, 

t ~2P ^-(^-) . (6.29) 

GW inspiral in the GC is of interest in anticipation of the planned Laser Interferometer Space 
Antenna (LISA) mission. LISA is expected to detect GW emitted by compact objects spiraling into 
MBHs 24 at cosmological distances (Barack & Cutler 2004; Gair et al 2004). The estimated rates are 
10" 9 -10" 8 WD yr 1 per galaxy (Hils & Bender 1995; Sigurdsson & Rees 1997; Freitag 2001 ; Ivanov 
2002; Hopman & Alexander 2005) and 10" 8 -10" 4 SBHyr -1 per galaxy (Miralda-Escude & Gould 2000; 
see review article on rate estimates by Sigurdsson 2003). For typical galactic density profiles, the rate 
depends only weakly on the relaxation time, because of the near cancellation of the ^-dependence 
of the scattering rate into the "inspiral cone" and that of the volume of the diffusive regime (§6.2.1). 
Furthermore, the rate depends only weakly on m because of the m/a relation (Eq. 1.1) (Hopman & 
Alexander 2005). Thus, the chances of detecting GW sources in any given galaxy are low, irrespective 
of its dynamical properties. 

The one exception is the Milky Way, where the proximity of the MBH makes it possible to detect 
even very low-mass objects with a weak GW signal (A£gw ocM 2 , Eq. 6.25) (Freitag 2003). Very low- 
mass MS stars dominate the population and their numbers far exceed those of compact objects (table 
2.1). The mean density of a star increases with decreasing stellar mass and reaches a maximum at 

<0. 1M Q , near the transition to brown dwarfs. Such stars can therefore survive the tidal field of a 



24 LISA is sensitive to GW in the frequency range 10 -4 < v < 10 1 Hz and is optimized for the range 10 3 < v < 
10~ 2 Hz, which corresponds to the frequency at the last stable circular orbit (LSCO) for MBHs of up to < 10 7 M Q 
(i'lsco ~ ^Gm/(3r s ) 3 /2ir = c 3 /(12V6nGm) ~ 2x 10- 3 (m/10 6 M Q )- 1 Hz). Extreme mass-ratio GW sources are 
of special interest because they are the cleanest probes of spacetime near a MBH, albeit with a weak signal (Eq. 
6.25; see review by Glampedakis 2005). This is to be contrasted with the yet intractable GR physics of 2 merging 
MBHs or NSs. The expected GW signal from inspiralling stellar mass objects at any given time is below the 
detector noise, but over the mission's lifetime of several years, ~ 10 orbital periods can be observed from an 
active source. It should then be possible to exttact the signal from the noise by utilizing the properties of the GW 
waveforms, which depend on the orbital parameters, in particular the eccentricity. 
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~ 10 Mq MBH and become GW sources. The mean number of active GW sources in the GC at any 
given time is Now ~lVo, where to is the inspiral time from an initial orbit with Po~ 10 4 s, the lowest 
frequency detectable by LISA. Freitag (2003) finds a capture rate of T, ~ 10~ 6 yr -1 for low mass MS stars 
in the GC (but see Alexander & Hopman 2003 who find a significantly smaller rate). This correspond to 
a few active GW sources at any given time in the GC, New ~ iVo > 1 for fo=2x 10 6 yr (Eqs. 6.28, 6.29 
with m=3 x 1O 6 M , M*=O.O5M , r p /r s = 10 and P = 10 4 s). It is not likely that LISA will detect GW 
from compact objects in the GC. The rates quoted above translate to Ngw^ 1 f° r the inspiral of SBHs 
and NSs and to N GW < 1 for WDs (Freitag 2003). 



6.2.3 Tidal heating, squeezars, and tidal capture 

Dissipative interactions of MS stars with the Galactic MBH proceed primarily via tidal heating, because 
the tidal disruption radius for a solar mass star lies well outside the event horizon, at r t ~ 10r^ . At that 
distance the tidal interaction dominates over the energy released by GW emission and the interaction 
can be treated to good approximation in the Newtonian limit. 

Inspiral by GW emission differs from inspiral by tidal heating in that the GW radiation is not local- 
ized in the star (the typical wavelength is Aqw ~ > fs ^>R*)- In contrast, the tidal energy extracted from 
the orbit is localized and thermalized inside the star. A star that is initially on a zero energy (parabolic) 
orbit has to dissipate many orders of magnitude its own binding energy in order to circularize, 

GmM* 1 / m \ 2/3 E± A 
E c = — — - ~ - — -± ~ 10 4 £* (6.30) 
2r p 2\Mj b 

where it is assumed that the circularization occurs at the initial periapse radius 25 . It is therefore unlikely 
that the star will survive tidal inspiral around a MBH. It will ultimately be destroyed, either by expanding 
beyond its tidally limited maximal size and breaking up (r t increases with R* until it overtakes r p ; Eq. 
3.2), or by exceeding its own Eddington luminosity and evaporating. Before that happens, the star will 
exist in a transient phase as a "squeezar", a star whose atypically high luminosity is powered by tidal 
interactions with the MBH rather than by its own nuclear burning (Alexander & Morris 2003). It should 
be noted that tidal capture is possible around an IBH because in that case E c /E± ~ O(10) (Eq. 6.30), 
and the tidal heat can be radiated away in O(10 4 ) yr (Hopman, Portegies Zwart & Alexander 2004). 

The energy deposition in the star in one periapse passage can be parametrized as (cf §3.4.2, Eq. 
3.30) 



AE t = T 2 (b 3 / 2 ) 
E± b 6 



**(0) 



(6.31) 



where possible stellar expansion is taken into account by the last term. As energy and angular momen- 
tum are exchanged between the orbit and the star, the star's mechanical and thermal properties change, 
and with them its response to the tidal force. Because of the feedback between the tidal interaction and 
the stellar properties, it is difficult to predict how the star and its orbit evolve in the course of tidal in- 
spiral. A similar problem is encountered in studies of the circularization of tidal capture stellar binaries 
in clusters. The results of such studies suggest two useful limits that likely bracket the true response: 
surface heating with radiative cooling (McMillan, McDermott & Taam 1987) in "hot squeezars", and 
bulk heating with adiabatic expansion (Podsdiadlowski 1996) in "cold squeezars". Hot squeezars dis- 
sipate the tidal heat in a very thin surface layer that expands moderately and radiates at a significantly 



25 If the impulsive relation holds between the tidal energy and angular momentum (Eq. 3.32), then r p = const 
and the initial r p is also the circularization radius. Conversely, if orbital angular momentum is conserved, then the 
circularization radius is 2r„. 



72 



T. Alexander / Physics Reports 000 (2005) 000-000 



increased effective temperature. This may apply to stars with large convective envelopes, where most 
of the energy in the oscillations is carried by eigenmodes that attain their maximal amplitude at the sur- 
face and dissipate there. Cold squeezars dissipate the tidal heat in their bulk, due to possible non-linear 
coupling between the high-amplitude, low-order oscillations that are directly excited by the tidal force, 
and an infinitude of high order modes which effectively dissipate instantaneously. Bulk heating leads to 
quasi-adiabatic, self-similar expansion at roughly constant effective temperature that doesn't rise much 
above its original value. 

With these and a few additional simplifying assumptions (motion in a Keplerian potential with 
constant r p ) it is possible to calculate the evolution of the stellar properties and the orbit. The orbital 
evolution is simple when the approximation Ais^const is valid (Eqs. 6.22, 6.23, 6.24). This is assumed 
for hot squeezars, but is not applicable to cold squeezars where AE t is not constant due to the adiabatic 
expansion (Eq. 6.31). In that case the evolution has to be calculated numerically. 

Numeric modeling of squeezar evolution (Alexander & Morris 2003) show that cold squeezars can 
brighten by up to AK ~ 3 mag. The temperature of hot squeezars can rise up to < 5 time the original 
effective temperature, and the star can brighten by up to AK^4 mag (Fig. 6.2). The luminosity and 
temperature increase could place low mass MS squeezars (M* >2M Q ) for the last < 10 4 yr of their 
life in the range of the S-stars. However, none of the S-stars for which an orbital solution is known 
approach the MBH close enough for tidal heating to be significant. Neither does the short squeezar 
cooling time allow for the possibility that these stars were tidally heated in the past and then scattered 
to wider orbits. Estimates of the inspiral rate for various tidal dissipation laws and their corresponding 
squeezar lifespans suggest that the mean number of squeezars in the GC at any given time is 0.1-1. 
Thus, the chances of detecting squeezars in the GC appear to be low. 

The end result of tidal inspiral (sometimes also called tidal capture) is disruption and the accretion 
of a sizable fraction of the stellar mass by the MBH. Tidal capture was hypothesized to be an important 
channel of stellar mass supply to a MBH, equal to, or exceeding by up to a factor of two the contribution 
of direct disruption (Frank & Rees 1976; Novikov, Pethick & Polnarev 1992; Magorrian & Tremaine 
1999). Simulations indicate that lower-mass MBHs (m< 10 7 M Q ) in low density galactic cores obtain 
most of their mass from direct stellar disruption (Murphy, Cohn & Durisen 1991; Freitag & Benz 2002). 
If tidal capture were indeed so efficient, this would have implied that such MBHs could be wholly 
constructed from the disruption of stars in the radius of influence, in seeming contradiction with the 
empirical m/a relation (Eq. 1.1), which holds on the much larger scale of the bulge (§1.1). The con- 
tradiction is resolved by noting that tidal capture, like any slow inspiral process, is strongly suppressed 
compared to direct infall. When the effects of scattering are properly taken into account (§6.2.1), it is 
found that tidal capture can contribute only an additional few percent above the mass supplied by direct 
disruption (Alexander & Hopman 2003). 

6.2.4 Tidal scattering 

Inspiral is very inefficient, and therefore most stars that are scattered to "near-miss" orbits pass by 
the MBH once and are then deflected to wider orbits. When the encounter is very close, the star will 
undergo an extreme non-disruptive tidal interaction and experience strong distortion, spin-up, mixing, 
and possibly some mass loss (Fig. 6.3). These effects may alter the star's subsequent appearance 
and evolution. While it is difficult to predict the long-term observational signature of such a "tidal 
scattering" event, plausible arguments suggest that the stellar luminosity will increase, photospheric 
abundances will be enriched by hydrogen-burning products (enrichment of 4 He, 14 N, 13 C and 26 Al and 
depletion of 12 C, 16 and 15 N) and the stellar colors will be bluer than expected for a normal star of that 
mass (Alexander & Livio 2001). It is interesting to note that precisely such abundance anomalies are 
observed in the M supergiant IRS 7, ~ 0.2pc from Sgr A* in projection (Carr, Sellgren & Balachandran 
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Figure 6.2: The evolution of the orbit (period and semi-major axis) and properties (temperature and 
K-band magnitude) a 1M Q hot squeezar in the GC, initially deflected into an orbit with r p = \.5r t {t =0) 
and Po= 1.4 x 10 4 yr, which corresponds to an inspiral time of ?o = 4.9x 10 5 yr. Time is shown relative 
to fdisrupt = 3.7x 10 5 , the time when r p = r t (^disrupt) an d the star expands to the point where it is tidally 
disrupted. At disruption the tidal luminosity is 642L Q , P=210yr and e= 1-2.3 x 10~ 4 . For n squeezars 
in the GC at any given time, the mean properties of the one in the most advanced stage of its inspiral can 
be read off the top axis. (Alexander & Morris 2003. Reprinted with permission from The Astrophysical 
Journal). 

2000). This suggests that IRS 7 has undergone mixing in excess of values predicted by standard models 
or observed in supergiants elsewhere, and has prompted Carr et al. (2000) to hypothesize that "extra 
mixing induced by rapid rotation may indeed be the fundamental difference between the evolution of 
massive stars in the GC and those elsewhere in the Galaxy." 

Whatever the exact nature of the after-effects of tidal scattering, their magnitude must increase with 
decreasing periapse. The definition of what constitutes a strong tidal interaction depends on the tidal 
effect that is of interest, and this in turn corresponds to some maximal periapse r p , beyond which the 
effect becomes insignificant. In spite of the smallness of the tidal radius, r t /rf, ~ 10~ 6 , the number of 
tidally scattered stars that accumulate inside the MBH radius of influence over the age of the Galaxy can 
be large (up to < 0.1 of the population). This can be shown qualitatively by noting that most of the stars 
that are deflected into low angular momentum orbits originate from c?~r/, orbits (Eq. 6.1 1) at the transi- 
tion between the diffusive regime and full-loss cone regimes. These orbits are approximately Keplerian 
and the DF is approximately isotropic, so that stars are deflected with equal probability to all directions. 
The cross-section for deflection into an eccentric orbit is then T,(<r p )ccr p , due to gravitational focusing 
(Hills 1975; Eq. 6.3). The rate at which stars are initially deflected into orbits with periapse smaller 
than r p but larger than r t is related to the direct disruption rate by Y s = (r p /r t — l)T t . This is also the 
tidal scattering rate to good approximation, because inspiral is a low probability process and most of 
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Figure 6.3: A sequence of snapshots (left to right, top to bottom) from an SPH simulation of a 10M Q , 
/?*=4.5/?0 star (represented by an n = 1.5 ideal gas polytrope) undergoing an extreme non-disruptive 
tidal interaction ("tidal scattering") as it passes near a 3 x 10 6 M Q MBH. Time is measured in units of 
the stellar dynamical time t*= 4840 s. The star passes near the black hole (located outside of the frame) 
on aparabolic orbit with r p = 100/?* = \.5r t =35Ars. After periapse passage (at t=2.2) the star appears to 
be on the verge of breaking in two (bottom left panel). However, by the end of the simulation, the two 
fragments coalesce, leaving a distorted, mixed and rapidly rotating bound object. 



these stars will be rescattered to a wider orbit (the star can also be scattered into the loss-cone, but the 
probability for that is smaller because loss-cone orbits occupy a smaller volume in phase-space than 
wide orbits). Another effect that will prevent once-scattered stars from interacting again with the MBH 
is its Brownian motion (§4.2.3). The apoapse of the scattered star's almost radial orbit will be r a ~2r/„ 
where the potential is no longer dominated by the MBH. The MBH executes Brownian motion relative 
to the stars inside 77, with an amplitude ~(r/,/ (/?*))((M*) /m) 5 / 6 r,3>r, and aperiod ~P(r/,) (Bahcall & 
Wolf 1976). Therefore, while the star is near apoapse the MBH position will shift relative to the stellar 
mas inside r/, without changing much the potential at r a , so that on its return to the center the star will 
miss the MBH. 

More accurate calculations (Alexander & Hopman 2003) show that the tidal scattering rate in a 
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power-law cusp n+ oc r a in the diffusive limit is, 
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(6.32) 



where P s ~0.9 is the survival probability against capture by the MBH and where typical values for the 
power-law index are (5=0.60 for a = 3/2, to <5=0.44 for a=7/4. For the full loss cone regime 5= 1, and 
so the actual effective power-law index is expected to lie in the range 0.5 < 5 < 1. 

A substantial fraction f of the total mass of a low-mass MBH (1O 6 -1O 7 M ) is supplied by tidally 
disrupted stars, which contribute each a fraction / m <0.5 of their mass to the MBH (§6.1.2). Simulations 
indicate that the total mass in tidally disrupted stars over the MBH lifetime lies in the range 0.15m < 
(ft/ f m )m< 0.65m, with the higher values typical of low-density nuclei where the population is weighted 
toward low mass stars (Murphy, Cohn & Durisen 1991; Freitag & Benz 2002). Assuming for example 
f t / f m = 0.25, 5 = 0.15 and a maximal periapse for strong tidal interaction of r p = l.5r t , it then follows 
from Eq. (6.32) that the mass fraction of tidally scattered stars within the MBH volume of influence is 
(ft/ fndiTs/Tf) ~ 0.07. The possibility that some of the unusual stars observed in the GC may be tidally 
scattered stars is intriguing. However, a better understanding of the long-term observable consequences 
of tidal scattering is needed to make progress on this issue. 

6.2.5 Star-disk interactions 

Drag against an accretion disk surrounding a MBH can dissipate orbital energy, as well as affect the disk 
itself (e.g. Ostriker 1983; Syer, Clarke & Rees 1991; Artymowicz, Lin & Wampler 1993; Vilkoviskij & 
Czerny 2002). There is no evidence today for an accretion disk in the GC. It is conceivable, however, 
that a dim, cold, low-viscosity "fossil" disk still exists around the MBH as a remnant of past periods 
of activity, perhaps associated with the creation and fragmentation of the gaseous disks that may have 
formed the observed star disks (§7.2.1), or that a small disk is continuously fed by cooling stellar winds 
from the young stars around the MBH (Cuadra et al 2005; §7.3). Stellar interaction with the disk can 
reveal its presence, or constrain its properties (extent, optical depth). If it is optically thick, the disk can 
eclipse stars behind it. Stars passing through the disk will likely emit X-ray and IR flares (Nayakshin 
& Sunyaev 2003; Nayakshin, Cuadra & Sunyaev 2004). The disk can reprocess the stellar UV light 
and re-emit it in the IR (Cuadra et al. 2003). However, observations do not support the presence of a 
disk. IR observations of the star S2 show no eclipses or IR flares and are consistent with the absence of 
a disk. Furthermore, the observed IR flares coincide to within a few milli-arcseconds with Sgr A* and 
display a ~ 17 min quasi-period (Genzel et al. 2003a), which can not be explained by passage through 
a disk. The temporal structure observed in the X-ray flares (2004) is likewise inconsistent with X-ray 
production by passage through a disk. 



The existence of young stars in the extreme environment so close to a MBH (§2.2) poses a challenge 
for theories of star formation and stellar dynamics (Genzel et al 2003b; Ghez et al 2003b). It has 
proved difficult to find a satisfactory explanation of how they could have formed so close to the MBH, 
or alternatively, of how they could have migrated inward from farther away in the course of their short 
lifespans. This is the so-called "paradox of youth". The question applies to any of the young stars in the 
inner parsec, but particularly so to the central cluster of the "S-stars", which exist a mere few hundredths 
of a parsec from the MBH. The problem of the young stars has become one of the major outstanding 
issues in GC research. 



7 The riddle of the young stars 
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The young population in the inner ~ 1 pc is often loosely described as the "OB-stars". This general 
designation can be misleading, as it fails to convey the significant systematic differences that exist in 
the population. An important open question is the nature of the connection, if any, between the S-stars 
inside ~*0.04 pc and the luminous emission line stars further out, on the 0.04-0.4 pc scale. While it is 
plausible to assume that these are different components of the same parent population, it should be noted 
the two groups have distinct locations, kinematics and stellar properties. The stars detected so far in the 
two young star disks at p~ 1"— 10" are luminous OB supergiants, giants and WR stars of various types 
(Genzel et al. 2003b; Paumard et al. 2001; Paumard et al. 2004). In contrast, orbital solutions (Ghez et 
al. 2005; Eisenhauer et al. 2005) show that the S-star orbits are randomly aligned and uncorrelated with 
the planes of the disks. Spectral identifications (Eisenhauer et al. 2005) reveal the S-stars to be B dwarfs 
(S2, the brightest, is a transitional O8V-B0V star; Ghez et al. 2003b). This distinction is important, 
because B-stars are more numerous, less massive (by about an order of magnitude) and longer lived (by 
1-2 orders of magnitude) than the luminous O-stars and WRs (Fig. 2.1). This significantly relaxes the 
constraints on models attempting to explain the central cluster. Thus there may actually be two separate 
issues: how to explain the presence of the OB stars on the > 0.04 scale, and an even more acute problem, 
how to explain the presence of the B-stars on the <0.04 pc scale. Some of the scenarios discussed below 
(§7.2) attempt to address both populations with a single solution, others focus only on one or the other. 

The stars in the central parsec can be described by stellar population synthesis models. Such models 
assume an IMF and a star formation history and use theoretically calculated stellar evolution tracks to 
follow the stellar population in time. Population synthesis models do not in themselves explain how 
the stars are born, or why the IMF and star formation history have a particular form. However, by 
fitting the models to the observed stellar population, it is possible to constrain these input parameters 
and estimate the total mass in the stars. Krabbe et al. (1995) find that the luminous early- and late-type 
stars can be modeled by a 7± 1 Myr old short burst of star formation with an IMF of dA^/dM* ocM~ 2 
between 1M Q and 120 M (that is, skewed to high masses as compared to the generic Salpeter IMF, 
dN+fdM* ocM~ 2 - 35 ), which formed ~3000 stars containing in total~ 1.5 x 1O 4 M . An additional, older 
star burst ~ 10 8 yr ago is indicated by the intermediate-mass asymptotic branch giants. This picture is 
consistent with the fact that the overall population in the central few parsecs is well represented by a 
continuous star formation model (Alexander & Sternberg 1999; §2.2). It appears that there is ongoing 
episodic star formation in the GC. This suggests that we are not observing the GC at a very special 
epoch, but rather that the presence of young stars there is a quasi-steady state situation. 

Any explanation of the luminous OB-stars on the >0.04 pc scale should account for their distinctive, 
disk-like distribution and kinematics and for the presence of the cluster-like structure IRS 13. Any 
explanation of the S-stars should account for three principal properties. 

1. The stellar properties. The stars appear to be entirely normal MS O8V/B0V to B9V stars, in 
terms of their (extinction corrected) luminosities, their absorption lines equivalent widths and 
line ratios (Eisenhauer et al. 2005). In particular, their rotational velocities are similar to those 
of nearby Galactic disk B-stars. The central cluster does not contain any star earlier than 08V, 
although such stars do exist outside the central arcsecond. 

2. The spatial concentration. The relative fraction of the young stars increases toward the center, to 
the near exclusion of any old giant stars stars in the inner ~0.02 pc (90% of all stars with K< 16 
mag in the central 0.02 pc are young main-sequence stars, Eisenhauer et al. 2005). 

3. The orbital properties. The orientations of the stellar orbits appear overall random, in marked 
contrast to the ordered planar rotation observed for the much more luminous emission line stars 
farther out. There are some statistically marginal trends observed in the orbital properties: higher 
than random eccentricity (Schodel et al. 2003) and a lower bound on the orbital apoapse of ~0.01 
pc. 
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The apparent normalcy of the B-stars is intriguing. This would seem to argue against any process that 
involves strong perturbations of the star, such as mergers (Genzel et al. 2003b), tidal heating (Alexander 
& Morris 2003) or stripping (Hansen & Milosavljevic 2003). However, this is by no means a conclusive 
argument, because it is not clear how efficiently a star relaxes after a major perturbation, and in particular 
how angular momentum is lost from the star or redistributed in it. 

A relevant comparison can be made with the properties of blue stragglers in clusters (MS stars that 
appear younger than the cluster's single age population), whose formation is thought to involve mergers, 
binary coalescence or 3-body interactions (Bailyn 1995). In spite of their violent birth, observations 
show that in many cases, late-B blue stragglers in intermediate age globular clusters, and early-B and 
O blue stragglers in open clusters do not seem to rotate faster than normal stars of the same mass in 
the field, in spite of the fact that, unlike lower-mass stars (later than F5V), they do not have efficient 
magnetic breaking (Leonard & Livio 1995 and references therein; §3.4.2). In this context it is interesting 
to note the claimed detection of a circumstellar disk around a blue straggler (de Marco et al. 2004), 
which may provide a breaking mechanism via magnetic anchoring. Neither is it clear whether collisional 
merger products should be mixed, and consequently, whether they should display unusual photospheric 
abundances. Lombardi, Rasio & Shapiro (1995,1996) find that the mixing by the collision itself is 
minimal (at least for lower mass stars near a globular cluster's turnoff-mass). Sandquist, Bolte and 
Hernquist (1997) find similar results. However, subsequent convection or meridional circulation can 
still induce mixing (Leonard & Livio 1995; Lombardi et al 1995), although this was not found to be 
the case in detailed numerical work by Sills et al. (1997) and Ouellette & Pritchet (1998). In summary, 
"absence of proof is not proof of absence"; if the analogy to blue stragglers is justified, then it appears 
that the lack of unusual spectral features in the S-stars does not place strong constraints on their origin. 



7.1 The difficulties of forming or importing stars near a MBH 

There is little evidence for star formation in the inner parsec at this time 26 . At present only low-density 
ionized gas is observed interior to the circum-nuclear molecular disk, which encircles the MBH at a 
radius of ~ 1.5 pc. Estimates of the molecular hydrogen density and the total mass in the disk vary 
from «h 2 ^5 10 6 cm~ 3 (tidally unstable, see Eq. 7.2 below) and M ~ fewxl0 4 M Q (based on molecular 
line ratios; Genzel et al. 1985; Jackson et al. 1993; Marr, Wright & Backer 1993; Marshall, Lasenby & 
Harris 1995) to «h 2 ~ 10 7 cm -3 (tidally stable) and M ~ 3x 10 5 M Q (assuming that the molecular clumps 
are gravitationally bound; Shukla, Yun & Scoville 2004). An additional < 1O 3 M may be contained 
in the gas and dust lanes (the "spiral") that extrude from circum-nuclear molecular disk inward (Liszt 
2003; Paumard et al. 2004; Paumard et al. 2004). Even if > 10 4 M Q of cold molecular gas existed in 
the past in the central parsec, as is implied by the mass in the observed blue and red giants (Krabbe 
et al. 1995), it is not clear how stars could condense from the gas so close to the MBH. The minimal 
proto-stellar cloud density that can resist the MBH tide is extremely high compared to values typically 
encountered in molecular clouds elsewhere in the Galaxy (n < 10 4 cm -3 ), where star formation normally 
takes place. For r < 1 pc, 



Some resolved infrared sources were initially explained as newly born stars still enshrouded in their dust 
cocoons (Ott, Eckart & Genzel 1999). However, more detailed observations later revealed these to be massive 
stars interacting with dust lanes extruding from the circum-nuclear molecular disk inward (Tanner et al 2002; 
Geballe et al. 2004; Paumard et al. 2004; Tanner et al. 2005). 




(7.1) 
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where m u is the atomic mass unit and where the gravitational potential of extended stellar mass is 
ignored. Forr> 1 pc, 

3(a-\)M fr\~ a 7 / M (lpc) \ f r \~ a _ 3 

mmn ~ 5 — ~4xl0 (a-l) —. cm , (7.2) 

A-nr\m u \r J \4x 1O 6 M J \ Ipc) 

where Mq is the total mass enclosed within ro and it is assumed that all the mass is distributed in a r~ a 
cusp. These high minimal densities exceed, or are just at the inferred densities in the circum-nuclear 
molecular disk. In addition, the molecular clouds in the GC have large turbulent velocities (~ 10km s -1 ), 
which may provide pressure support against fragmentation, and unusually strong magnetic fields (~ 1 
mG, e.g. Aifken, Moore & Roche 1998) are observed in the GC, which, if they permeate the clouds, 
would provide additional support against fragmentation and collapse (Moms 1993). The unavoidable 
conclusion is that if the young stars well inside the central parsec were formed locally, then they must 
have done so by a different mechanism than the collapse of self gravitating cold molecular gas clouds 
that occurs in normal star forming regions. 

The young stars are also too short-lived and too light to have formed far from the MBH, where the 
tidal field is weak, and then to have migrated in by dynamical friction. The "collection basin" for the 
migration of young MS stars of mass of M+ = 15M Q and lifespan ?* = 2x 10 7 yr (§2.2) is only maxrdf ~ 
0.2pc (Eq. 3.27), and for stars of mass M* = 3M and lifespan ?*=4x 10 8 it is only maxr^f ~0.5pc. 



7.2 Proposed solutions 

The solutions proposed so far for the riddle of the young stars (see reviews by Genzel et al. 2003b; Ghez 
et al. 2005) fall into three main categories: unusual modes of star formation near the MBH; rejuvenation 
of old stars from the local population; and dynamic migration or capture from farther out, where stars 
can form. While each has some attractive features, none is quite satisfactory. The paradox of youth 
remains unsolved at this time. 



7.2.1 Unusual in-situ star formation 

One class of in-situ star formation models invokes external pressure to trigger cloud collapse. Because 
the clouds are stabilized against collapse by their high turbulent velocities and magnetic fields, star 
formation in such clouds will be skewed toward massive stars, as are observed in the GC. Cloud-cloud 
collisions and the resulting shocks and cooling could in principle initiate fragmentation (Moms 1993). 
The near identical age of the two star disks could perhaps be explained by two colliding clouds, each 
fragmenting to form a disk (Genzel et al. 2003b). However, it is not clear that the required very high 
compression ratio (Eq. 7.1) can be achieved in such collisions. An additional objection to any scenario 
that postulates the existence of dense clouds near the MBH is that those clouds that are dense enough 
to resist the tidal forces at r < 2 pc will also be Jeans-unstable and will fragment and presumably form 
stars on that large scale before reaching the center (Vollmer & Duschl 2001). 

Alternatively, Moms, Ghez & Becklin (1999) propose a scenario of recurrent nuclear activity regu- 
lated by a limit cycle. At the current phase of the cycle, the circum-nuclear molecular gas disk, which 
has a central cavity inside ~ 1.5 pc, is prevented from filling the cavity by the radiation pressure from 
the hot massive stars in the center. However, once these stars die (in < 10 7 yr), the internal viscosity of 
the disk will set an inflow. The central cavity will be filled and the gas will reach the center, setting off 
a strong burst of luminous accretion. The radiation pressure will shock and compress the inner parts of 
the gas disk, presumably triggering an intense phase of star formation. Stellar winds and the combined 
radiation from the young massive stars and the accretion on the MBH will evacuate a central cavity, 
which will remain empty of gas as long as there are enough massive luminous stars inside, thereby 
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starting the limit cycle again. This scenario has yet to be studied in detail. Here again, the problem is 
that it is unclear whether radiation pressure can lead to the very high compression ratio that is required 
for fragmentation. 

Another class of models for in-situ star formation propose that the young stars were formed by the 
fragmentation of a gaseous disk around the MBH that gradually grew in mass (perhaps fed by tidally 
disrupted molecular clouds, Sanders 1998) to the point it became self gravitating (Levin & Beloborodov 
2003; Milosavljevic & Loeb 2004; Nayakshin & Cuadra 2005). The outer parts of thin massive disks 
are susceptible to fragmentation (Paczyhski 1978), so much so that it is, in fact, difficult to explain 
theoretically the observed existence of extended disks (Goodman 2003). While such a disk does not 
exist in the GC today, Milosavljevic & Loeb (2004) point out that masering disks observed around 
MBHs in other galaxies (disks emitting coherent radio emission due to population inversion by shocks) 
have properties similar to those needed to create the star disks in the GC. 

The criterion for disk fragmentation due to self-gravity (Paczynski 1978) can be expressed, up to 
factors of order unity, as a tidal limit, p='E/h>m/r 3 = p c , where £ is the disk's surface density. The 
disk's scale height h ~c v /J7 is set by the maximal height z an atom can reach when launched from 
the disk mid-plane with the sound speed c s against a vertical acceleration a z ~ £l 2 z (z <C r), where 
fi = y/Gm/r 3 is the local Keplerian frequency. At the critical density, Q equals the disk's gravitational 
free-fall rate, O = \J Gp c = t^ 1 . The Jeans (minimal) radius, Rj, and mass, Mj, of a collapsing fragment 
are set by the condition that the free-fall time be shorter than the sound crossing time (the time for 
the propagation of a gravity-resisting pressure adjustment), Rj > c s tg, and Mj > p c c 3 s /il 3 = c^/(G 2 £ c ). 
In addition, the gas must cool faster than the dynamical time to offset the increased pressure support 
from the contracting gas. For the mass and distance scales of the GC and for reasonable assumptions 
about the disk properties, M/~fewx 1M Q (Levin & Beloborodov 2003). The collapsing proto-stellar 
cloud then begins to accrete gas from the disk within its tidal radius, ~ (Mj /m) l l 3 r. The available 
mass there exceeds Mj by a few orders of magnitude (Levin 2003; Goodman & Tan 2004). Due to 
the differential rotation in the disk, the proto-stellar cloud acquires its own mini-accretion disk, which 
may itself become unstable and fragment to form groups of stars. Milosavljevic & Loeb (2004) propose 
that IRS 13 is such a gravitationally bound group of ~2500M Q (including low luminosity stars not yet 
observed), and that the S-stars were internally scattered from such groups by very hard encounters with 
binaries. 

Nayakshin & Cuadra (2005) derive limits on the initial total masses (stars and gas) of the two star 
disks observed today: a lower limit of ~ 10 M Q from the requirement of self-gravitation instability, and 
an upper limit of 10 5 M Q based on the magnitude of velocity dispersion in the outer disk that is caused 
by the potential of the inner disk. They are thus able to rule out the possibility that the young stars were 
formed by gas accretion on low mass stars that got captured in the disk. This would have required disk 
masses well in excess of the upper limit. 

The disk fragmentation model is a promising scenario for the formation of the star disks and the 
massive young stars on the >0.1 pc scale, but it is less clear whether it can account for the S-cluster as 
well. According to their stellar contents, the two star disks are contemporaneous to better than 1 Myr, 
and only ~ 5 Myr old (Krabbe et al. 1995 ; Genzel et al 2003; Paumard et al. 2005, in prep.). Any 
model that associates the S-stars with the star disks has to explain why it is only B-stars that are found 
in the S-cluster, and not any of the more massive stars that exist in the disks. Furthermore, if the B-stars 
are associated with the disks, then the relevant time constraint for redistributing the orbits is the age of 
the disks, and not the considerably longer MS lifespan of a B-star. 
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7.2.2 Rejuvenation 

The inner ~ 0.02-0.03 pc of the GC, the domain of the S-cluster, roughly coincides with the region 
where the stellar density is so high that stellar collisions occur more than once over the lifetime of a MS 
dwarf (Alexander 1999; §3.4). The strong tidal field of the MBH may also affect the stellar structure 
(§6.2). Rejuvenation models seek to explain the S-cluster in terms of low mass, long-lived stars that had 
enough time to migrate to the center from their formation sites far from the MBH (however, the central 
number density of such stars will be suppressed by mass segregation). 

All rejuvenation models proposed so far have serious problems in explaining the S-stars. The physi- 
cal processes were already discussed above in some detail, and are listed here briefly. Envelope stripping 
collisions (Hansen & Milosavljevic 2003) will reveal a hot core, but much fine-tuning is required for the 
bare core to masquerade as an apparently normal B-star (§3.4.1). Tidal heating (§6.2.3) requires a very 
small periapse, r p <2r t to be effective (Alexander & Morris 2003). Even the smallest measured periapse 
to date, that of star SO- 16 with r p ~ 600rj ~ 30r, (Ghez et al. 2005) is much too large for tidal heating 
to be of any relevance (tidal heating and scattering to a wider orbit can be ruled out as an explanation 
for its present luminosity since the scattering timescale is much longer than the stellar cooling time). 
Successive mergers of low-mass stars (Genzel et al. 2003b) are inefficient because the mass retention 
in a high velocity collision is typically very small (§3.4.3). Exotic objects formed by the capture of a 
compact remnant inside a normal star will probably not look like B-stars (§3.4.3). 

7.2.3 Dynamical migration 

Stellar dynamics are generally better understood than star formation, or the effects of far from equilib- 
rium conditions on stellar structure and properties. Therefore models that invoke dynamical processes 
to explain the S-stars tend to have more definite, falsiflable predictions. The orbits and the stellar spatial 
distribution thus appear to be a more promising way of discriminating between models. One class of 
dynamical models seeks to accelerate the dynamical friction (§7.1) by attaching the young stars to a 
massive "anchor". This allows the stars to form far from the MBH, where star formation is not strongly 
inhibited, and yet sink into the center within their short lifespans. 

A natural candidate for such an anchor is a dense, young, star-forming cluster (Gerhard 2001; §2.1) 
like the Arches and Quintuplet clusters, which lie within tens of parsecs from the center, and contain 
a young population closely resembling that near the MBH, in particular He stars (Figer et al. 1999). 
Such a cluster, being a compact massive object (M~ 10 4 M Q ), will undergo dynamical friction and can 
sink to the center in 10 7 yr from a distance of ~ 5 pc (Eq. 3.27). Much more massive clusters than are 
seen today in the GC, with M~ 10 6 M Q will be able to sink in from a distance of ~ 50 pc. The problem 
with this scenario is that the cluster will be stripped and dissolved by the tidal field of the GC before 
it reaches the central parsec (Portegies Zwart, McMillan & Gerhard 2003), unless it has an extremely 
high central density, p~ 10 8 M Q pc" 3 (Kim & Morris 2003). Alternatively, the center of the cluster may 
be stabilized against tidal stripping by a central IBH (Hansen & Milosavljevic 2003). 

The scenario of the dissolving cluster with an IBH proposes that a very dense stellar cluster con- 
taining an IBH, perhaps formed by runaway mergers in the cluster core (e.g. Portegies Zwart et al. 
2004), sinks rapidly to the center. The stars most tightly bound to the IBH avoid tidal stripping until 
the IBH reaches ~0.1 pc from the MBH, where they are deposited in a disk-like configuration. Some 
of them are subsequently scattered to tighter eccentric orbits by repeated interactions with the orbiting 
IBH. An analogous mechanism operates in the Solar System, where Jupiter scatters comets from the 
outer Solar System to tightly bound inner orbits. Subsequent, more detailed simulations (Levin, Wu & 
Thommes 2005) showed that a dissolving cluster with an IBH can deposit its last remaining stars in a 
fairly thin ring (if moving on a circular orbit), but that even repeated scattering of these stars by the IBH 
do not lead to the formation of a tightly bound cluster such as the S-cluster. Hansen & Milosavljevic 
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(2003) hypothesize that the S-stars are in fact O-stars from the dissolving cluster that were captured in 
very bound orbits by a close, nearly disruptive encounter with another star, which stripped their outer 
envelopes and whittled them down to B-star masses. The general concept of this model is supported by 
recent claims for an IBH embedded in the IRS 13 "cluster" (Maillard et al. 2004; but see dissenting view 
by Schodel et al. 2005). The weakness of this scenario for explaining the S-stars is that, lacking quan- 
titative calculations, it is unclear whether this method can actually form such a tightly bound cluster, 
and it is also unclear whether collisional stripping is efficient enough and consistent with the apparent 
normalcy of the B -stars. 

Another objection to this scenario is the large number of massive young stars that are implied by a 
dissolving M> 1O 5 M cluster (Kim, Figer & Morris 2004; Giirkan & Rasio 2004), and which are not 
observed in the GC. Giirkan & Rasio (2004) show that a very massive (M > 10 6 M Q ) and dense cluster 
starting at lOpc can sink in less than 10 6 yr to the central 0.5 pc, bringing with it 10 4 M Q of > 10M Q 
stars while undergoing core collapse and forming a 5000 M & IBH. However, if a lower mass cluster is 
assumed (and a correspondingly smaller initial radius, r& ~ 5pc) to mitigate the problem of the over- 
abundance of massive stars, then in order to carry stars into the inner parsec, the IBH must have a large 
mass relative to the cluster, m> 10 4 M Q and m/M> 0.1. This is about 2 orders of magnitude larger than 
can feasibly grow in runaway mergers, m/M ~ 10~ 3 (Giirkan, Freitag & Rasio 2004; Giirkan & Rasio 
2004). 

7.2.4 Exchange capture 

Another class of dynamical scenarios are those that invoke 3-body interactions to capture the B-star on 
a tightly bound orbit around the MBH. 

The massive binary exchange scenario (Gould & Quillen 2003) postulates that the B-star originally 
had a very massive binary companion (~ 100M Q ). The binary presumably originated in a radially 
infalling, disintegrating massive cluster. Its radial orbit brought it close to the MBH, where in the course 
of a 3-body exchange interaction the B-star switched partners and became bound to the MBH. Gould & 
Quillen (2003) find that the probability per binary for capture in an orbit like that of the star S2 is of the 
order of a few percent. The weakness of this scenario is the very low joint probability for (i) having a 
cluster on a radial orbit, (ii) containing a very massive star, (iii) paired in a binary with a much lighter 
secondary (iv) with the right orbital parameters for 3-body exchange with the MBH. It also does not 
provide an explanation why it is that only B-stars are captured, and not more massive stars, and why 
there is a lower bound on the apoapse. 

The exchange capture with SBHs scenario (Alexander & Livio 2004) proposes that the B-stars were 
originally formed far away from the MBH where normal star formation is possible (not necessarily in a 
massive cluster) and were then deflected into eccentric orbits that intersect the dense central concentra- 
tion of SBHs (with masses of ~ 7-10M Q ), which sank to the center from the inner ~5pc due to mass 
segregation over the lifetime of the GC (§3.3). The high concentration of the SBHs is also responsible 
for the collisional destruction of any tightly bound old red giants at the very center. Occasionally one of 
the B-stars passes very close to a SBH, knocks it out and replaces it in a bound orbit around the MBH (a 
3-body exchange involving a MBH-SBH "binary" and the single B-star). This mechanism can naturally 
explain several of the properties of the S-cluster. The central concentration of the S-cluster follows that 
of the highly segregated cluster of stellar BHs. The S-cluster is composed of B-stars because they are 
most closely matched in mass to the stellar BHs and so have the maximal probability for exchange. 
The lower bound on the apoapse, ~0.01 pc corresponds to the point where an exchange requires that 
the B-star pass so close to the SBH that it is tidally disrupted. The exchange mechanism also predicts 
a trend toward high eccentricities in the orbits of the captured stars. The weakness of this scenario is 
the high required central concentration of SBHs and the large required number of B-stars on eccentric 
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orbits. Detailed calculations of the exchange capture cross-section and rate indicate that the SBH con- 
centration must be close to the drain limit (Eq. 3.29) in order to sustain a steady state population of 
captured B-stars. These requirements may be relaxed if the source of B-stars is a dissolving cluster that 
approaches the SBH cluster (Giirkan & Rasio 2004), however, detailed calculations of this possibility 
have yet to be performed. 

The stripped AGB giant capture scenario (Davies & King 2005) suggests that the S-cluster stars are 
the stripped core of very luminous asymptotic giant branch stars, which are captured at typical S-star 
periapses and periods by the dissipation of orbital energy required for tidally stripping the extended 
giant envelopes. The criticism leveled at this scenario is that it requires an unrealistically large number 
of AGB giants to sustain a steady-state S-cluster population (total mass consumed in AGB stars larger 
than the stellar mass in r/,); that much fine tuning is required for the stripped cores to masquerade as 
B stars; and that the work done to strip the envelope will be taken out of the envelope's orbital energy 
rather than that of the core (Goodman & Paczyhski 2005). 

7.3 Feeding the MBH with stellar winds 

Given the yet unexplained fact that young massive stars do exist very close to the MBH, it is natural 
to consider the role that the mass ejected by their strong stellar winds may play in the MBH accretion. 
The massive He stars with their copious mass-loss winds can supply mass for accretion at a rate of 
M w > 10 M©yr . Estimates of the gas density and temperature from Chandra X-ray observations 
(Baganoff et al. 2003) imply a Bondi accretion rate (Bondi 1952; Melia 1992) of M B ~ 10" 5 M yr" 1 < 
M w , if the accretion is spherically symmetric. Accretion at this rate by a geometrically thin, optically 
thick accretion disk, which has a typical radiative efficiency of r/^O.l (Shakura & Sunyaev 1973; 
§6. 1 .2), can be ruled out since it would produce a luminosity of L B ~ r]M B c 2 ~ 10 41 erg s" 1 , about 10 5 
higher than observed. The implications are that (1) the He stars alone can easily supply all the mass 
required for accretion, and (2) Sgr A* is extremely under-luminous compared to available gas supply 
rate. It is believed that the dimness problem of Sgr A*, which is common to many MBHs in quiescent 
galactic nuclei, can be explained by radiatively inefficient accretion. The detection of linearly polarized 
sub-mm emission from Sgr A* (Aitken et al. 2000; Bower et al. 2003) as well as the observed spectral 
energy distribution of Sgr A* suggest that only M~ 10~ 8 MQyr _1 CM^ is accreted on the MBH, with 
a low efficiency of 77 ~ 10~ 3 . The rest of the mass is presumably carried away by a global outflow 
(see review by Quataert 2003). Note that Mtu <C m, so the present accretion rate cannot explain the 
mass of the MBH. Either the MBH was accreting at much higher rates (at the level of the Bondi rate 
or higher ) for most of its past, or else it grew by channels other than gas accretion, for example, a 
substantial fraction can be supplied by tidal disruption of stars, which can account in steady state for 
~O.5M r,? // ~2.5xlO 5 M ~O.lm (§6.1.1). 

A more realistic picture of the flow of the stellar winds into the MBH must take into account the 
discrete distribution of stars (Coker & Melia 1997; Rockefeller et al. 2004). Cuadra et al. (2005) note 
that stellar mass lost by slow winds (v<300kms _1 ) or in a faster wind, but in the direction opposite 
to the star's orbital motion, will not be shocked to very high temperatures and can cool on a dynami- 
cal timescale. They simulate the gas flow from the winds of an ensemble of orbiting stars in the inner 
~0.35pc (roughly the region where the 2 stellar rings lie, §2.2) and find that a 2-phase gas distribution 
develops, with hot, X-ray emitting gas in the inner 1" around the MBH and a cold, fragile mini-disk fur- 
ther out (containing ~ 10M Q ). Cold blobs occasionally detach from the disk and flow inward, leading, 
if they do not evaporate, to a variable accretion rate with a variability timescale of 10 2 — 10 3 yr. 

While the low level, quasi-steady state accretion of Sgr A* may be supplied by the shocked winds 
of massive stars in the inner parsec, a time dependent component may be contributed directly by stars 
on eccentric orbits that approach within ~ 10 3 rs of the MBH. Loeb (2004) suggest that a fraction of 



T. Alexander / Physics Reports 000 (2005) 000-000 



83 



the stellar wind in the direction oriented toward the MBH can fall ballistically into the MBH, thereby 
possibly introducing variability that correlates with the orbital phase. 

8 Outlook 

8.1 Progress report 

Table (8.1) summarizes the progress achieved to date (mid 2005) by the study of stars and stellar pro- 
cesses near the MBH in the GC, as well as the many remaining challenges. Several subjects are listed, 
broken down to sub-topics. For each topic, progress in three categories, theory, predictions and ob- 
servations, is graded on a qualitative scale, from "unknown" to "known with confidence". The theory 
category refers to the level of basic understanding of the particular object or process and the degree of 
confidence that it should exist or play a role in the GC. The predictions category refers to the ability to 
predict the observable consequences in the GC, either from a basic theory or from established empirical 
relations, assuming that the object or process is indeed relevant. The observations category refers to 
the extent at which observations were able to provide information on the matter (either supporting or 
refuting existence). 

For example, concerning the nature of the dark mass, the theoretical expectation is that the dark 
masses in galactic centers are MBHs (based, among others, on the close analogies with Galactic "micro- 
quasars" (SBHs accreting from a binary companion) which are not composite objects such as clusters 
and most probably are not made of exotic dark mass). Furthermore, a basic theory of BHs exists. The 
theoretical predictions of the effect of a MBH on stellar orbits are firm, and there are empirical expecta- 
tions for its mass based on the m/a relation. Even for the less- well defined exotic dark mass alternatives 
it is possible to predict that there should be deviations from Keplerian motion. The observations of stel- 
lar orbits provide a compelling case for a MBH. Consequently, theory, predictions and observations for 
this topic are all graded as "known with confidence". Another example, where the situation is less clear, 
concerns stellar collisions and mergers. While it is theoretically known with confidence that stellar col- 
lisions should occur in a high density environment, the predicted observational signature is uncertain, 
as it depends on the details of stellar physics. Nevertheless, there are promising, albeit still uncertain 
observational indications that giants are collisionally destroyed very near the MBH. 

It should be emphasized that it is often the case that there isn't one definitive theoretical prediction, 
either because of the complexity of the process, or because the predictions depends on various unknown 
parameters. The scores listed here reflect an admittedly subjective assessment of the overall robustness 
of the theoretical understanding and predictions. Likewise, the question whether or not available ob- 
servations can come to bear on a particular topic is tied to the theoretical progress. It is quite possible 
that relevant observations are already available, but remain "hidden in plain sight" due to the lack of a 
theoretical framework. 

8.2 Future directions 

The main achievements resulting to date from the study of stars near the Galactic MBH (table 8.1) in- 
clude the determination of the primary parameters of the central dark object (its nature, mass, location 
and distance), the detection of an extremely high density stellar cusp, and the observational characteri- 
zation of the puzzling stellar population near the MBH. There are still many challenges in making the 
connection between these discoveries and the underlying theories. This review concludes with a brief 
survey of some of the next research frontiers and the theoretical and observational advances that are 
required for making progress on these issues. 
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Table 8.1: Score card and progress report (mid 2005) 



Category 



Score " 
Theory/Predict./Obs. 



Remarks 



Primary parameters of the dark mass 



Nature of DM 
Position of DM 
BH mass 
BH spin b 
Binary BH 



V V 



?? 

V 
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MBH, from lower limit on density (§4.2). Dynamical 
upper limits on extended component (§4.2.3). 
Apparent position from orbital solutions (§4.1.2). 
Line-of-sight distance from radial velocities (§4.1.3). 
No basic theory, but mass consistent with m/a predic- 
tion (§4.1). 

! Quasi-period in IR accretion flares (Genzel et al. 

2003a; Aschenbach et al. 2004). 
? IRS 13? Bounds on secondary mass, distance (§4.3). 



Properties of the stellar system 

Density distribution J 
Dynamical properties ! 



Population composition 
Star formation mode 



?? 

99 



v 7 



? 

99 



v 7 
v 7 

v 7 

9 



High density stellar cusp observed (§3.2). 

Vel. dispersion, anisotropy measured on all scales 

(§4.1.1). 

Mixed old and young stellar populations (§2). 
Many ideas, but still no definite conclusion (§2.1). 



Stellar dynamical processes 

Relaxation 



Mass segregation 

Collisions/mergers 
Tidal spin-up 
Star disks 



v 7 

v 7 
! 

9 



? 
? 

99 



! 

99 

v 7 



Observed cusp slope and velocity dispersion of old 
population consistent with relaxation (§3.2). 
Strong segregation expected. Some evidence for over- 
abundance of compact objects (§3.3). 
Depletion of bright giants toward center (§3.4). 
Evidence of (rotational?) mixing (§3.4.2) 
Star formation in fragmenting gas disks? (§2.2). 



Stellar interactions with the MBH 

Tidal disruption 

Tidal scattering/heating ! 

Interactions with cold disk ?? 



? 

?? 

v 7 



Evidence inconclusive (§6.1). 

Process inevitable, signature uncertain (§6.2.3, 
§6.2.4). 

No evidence of cold disk (§6.2.5). 



Post Newtonian physics 

Gravitational redshift J 

Orbital precession ! 

Frame dragging ! 

Gravitational lensing ! 

Gravitational waves ! ! 



99 



Awaiting more orbital data (§5.1.1). 
Awaiting short period, better astrometry orbits (§5.1). 
Proposed detection depends on untestable assump- 
tions about the star S2 (§5.1). 
Expected lensing probability low (§5.2). 
Awaiting future detectors (§6.2.2). 



" ' ^J' known with confidence; '!' promising but uncertain; '?' uncertain; '??' speculative; ' — ' unknown or not applicable 
* Listed here for completeness. This estimate was derived from gas accretion properties, not stellar processes. 
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A yet unrealized goal is the detection of post-Newtonian effects in the stellar orbits (§5.1). The 
prospects of progress will be much improved if stars on even tighter and more eccentric orbits than 
those observed today are detected. Such stars, if they exist, must be less massive and less luminous 
than the faintest stars already detected near the MBH, late B-stars (M* > 3M®). The question whether 
such stars are present near the MBH is closely linked to the issue of mass segregation (§4.2.3). If, as 
is anticipated, most of the objects near the MBH are massive remnants, then the density of low mass 
stars there will be significantly suppressed. Furthermore, collisional destruction will limit the lifespan 
of tightly bound stars (§3.4.1). 

To detect and track stellar orbits in a field that becomes increasingly cluttered as more faint stars are 
revealed ("source confusion"), it will be necessary to increase the photometric, astrometric and spec- 
troscopic sensitivities. This can be achieved by much larger telescopes than the existing 8-10 m class 
telescopes, which will have a correspondingly smaller diffraction limit (e.g. the proposed Thirty Meter 
Telescope (TMT), for GC applications see Weinberg, Milosavljevic & Ghez 2005). Alternatively, opti- 
cal/IR interferometers such as the Very Large Telescope Interferometer (VLTI), the Keck interferometer 
and the Large Binocular Telescope (LBT) will achieve very high angular resolution and may also reach 
high photometric sensitivity with adaptive optics. Monitoring a short period star will require frequent 
sampling to avoid losing track of the star in the dense field, and to obtain good coverage of the orbit, 
especially near periapse, where the post-Newtonian deviations are largest (§5.1). 

A closely related problem is the detection and characterization of the extended dark mass around the 
MBH. Orbital monitoring can also detect deviations from Keplerian orbits due to a smooth distribution 
of dark mass or perturbations due to interactions with discrete objects (compact remnants). It may also 
reveal gravitational lensing events (§5.2) involving both the MBH and the cluster of compact remnants 
or faint stars around it (Alexander & Loeb 2001; Chaname, Gould & Miralda-Escude 2001). On a scale 
larger than probed by the orbits, X-ray surveys have already discovered an over-abundance of transient 
X-ray sources in the central parsec, which are thought to be NSs or SBHs accreting from a binary 
companion (Muno et al. 2005). Deep surveys at short radio wavelengths may probe the distribution of 
millisecond pulsars, which, as long-lived light test particles, are expected to be depleted in the center 
(Chaname & Gould 2002). The distribution of other long-lived light test particles, such as red clump 
giants, should also be sensitive to the degree of mass segregation. Theoretical predictions of the details 
of mass segregation around a MBH are only available for highly idealized situations (Bahcall & Wolf 
1977), and the application to real systems is unclear. Dynamical simulations of stellar systems around 
a MBH (e.g. Freitag & Benz 2002; Preto, Merritt & Spurzem 2004; Baumgardt, Makino & Ebisuzaki 
2004) that take stellar evolution and continuous star formation into account will be needed to confirm 
and refine these predictions. 

The detection of evidence of past strong interactions between stars and the MBH and between 
stars themselves (§6, §3.4) will require large area surveys, in particular spectroscopic surveys, to look 
for tell-tale peculiar stellar properties. There are however still no robust theoretical predictions of the 
appearance of stars that were subjected to major perturbations. Any progress on this complex problem 
will be very relevant for such searches, and conversely, the stellar population in the inner GC can be 
used test these ideas. 

The attempts to solve the riddle of the young stars (§7) raise many interesting possibilities about 
the contents and history of the inner GC. A deep spectroscopic survey of the inner ~ 15" will help 
disentangle the different stellar populations there and provide clues about the history and origin of the 
young stars and the properties of the star disks. Such a survey can establish whether or not apparent 
concentrations of stars, such as IRS 13, are indeed self-bound objects, perhaps harboring an IBH. 

The Galactic MBH provides a uniquely accessible laboratory for studying in detail the connections 
and interactions between a massive black hole and the stellar system in which it grows; for investigating 
the effects of extreme density, velocity and tidal fields on stars; and for using stars to probe the central 
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dark mass and post-Newtonian gravity in the weak- and strong-field limits. These issues are relevant 
for understanding the MBH phenomenon in general. The wealth of observed phenomena in the GC 
provides the impetus to study stellar processes in the extreme environment near a MBH. As shown in 
this review, such studies prove to be very fruitful, yielding valuable insights even in cases where they do 
not fully succeed to explain the observations. Many questions remain — it is quite likely that the most 
exciting science still lies ahead. 
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